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CA NNOUNCEMENT 


An agreement has been effected between the Johns Hopkins University 
and the American Mathematical Society for codperative publication of the 
AMERICAN JOURNAL OF MATHEMATICS. The Journal will be enlarged so as 
to help to meet the rapidly growing demands for the reasonably prompt publi- 
cation of memoirs. There is no intention of changing the policy. 

Substantial theories have, of course, the first place. As far as possible 
illustrations of such theories should be included, as they facilitate enormously 
the comprehension of the reader and thereby the spreading of the theories. 

In the second place the Journal will welcome the converse process of 
elucidating a specific problem either with the aid of available theories or as 
pointing to further theories. 

In the third place the Journal will gladly accept the applications of 
theories; i. e., the illumination of methods however general by specific and 
easily comprehended cases. 

It is hoped that Applied Mathematics may be a valued feature in the 
next lap of the Journal’s career. 

G. D. BrrxHorr, 
A. B. Cosiz, 
ABRAHAM COHEN, 


G. C. Evans, 
FRANK Mortey, Chairman. 
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Stability and the Equations of Dynamics.” 


By Gerorce D. BIRKHOFF. 


Introduction. 


For a long time attention has been directed towards the equations of 
dynamics, and more and more there has been the tendency to take these as 
being derivable from a special type of variational principle and as consequently 
possessing ‘ Hamiltonian’ or ‘canonical’ form. In consequence of the fact 
that the mathematical treatment of these equations involves difficulties due 
to the lack of flexibility of the underlying group of transformations, a very 
large formal literature has grown up. On the other hand, physicists, realizing 
the sufficiency of such equations for most applications and impressed with the 
purely mathematical developments, have yielded a position of the very first 
importance to them, just as have the mathematicians because of the physical 
applications. 

For these reasons there has been a disposition to attribute an almost 
mysterious significance to the Hamiltonian equations. 

Now there is one property of these equations, first explicitly and ade- 
quately recognized by Poincaré, which is of obvious significance. If there is 
given an equilibrium position or a periodic motion of a dynamical system of 
this type, which satisfies the simplest necessary conditions for stability, and 
if there are no relations of commensurability between the underlying periods, 
then there will be complete formal trigonometric stability. This fact is of 
course of great physical interest, and throws great light upon the inward 
nature of the stability existing in the solar system for instance. 

It is the main contention of this paper that essentially the only signi- 
ficance of the Hamiltonian (and thus the Lagrangian) equations is that they 
possess this characteristic property of complete’ formal stability. In other 
words, any set of n equations possessing complete formal stability at an equi- 
librium point, for example, may be given Hamiltonian form by appropriate 
formal change of variables. 

In the present paper attention is primarily directed to the case of a 


*See a preliminary note in the Comptes rendus for September 20, 1926. The 
present paper was read at the Philadelphia meeting of the American Mathematical 
Society on December 29, 1926. 
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system of equations of even order, although it is proved that the odd order 
case essentially reduces to that of even order. Furthermore, the treatment 
is restricted to the typical case of the equilibrium point, although I expect to 
treat the analogous important case of a periodic motion elsewhere, in so far 
as the generalization is not an obvious one. 

The precise nature of the Hamiltonian equations is not adequately ex- 
pressed in the mere statement that the property of stability is their one 
characteristic feature. Indeed, except from a purely formal point of view, 
these equations form a particular kind of completely stable systems, because 
of their convergence properties. The interrelation between the various kinds 
of stability of systems of ordinary differential equations, the nature of the 
trigonometric series employed and their degree of accuracy, the formal reduc- 
tion to Hamiltonian or Pfaffian form in the case of complete stability, must 
be studied in detail before the very interesting situation which exists can be 
rightly comprehended. It is my purpose in this paper to clarify this situation 
as far as possible. 

Many other problems are immediately suggested. I will only refer to one. 
In recent work on the theory of relativity there have been notable attempts 
made by Eddington, Einstein, Hilbert and Weyl to derive the general partial 
differential equations of physics from a ‘ Hamiltonian principle.’ There arises 
now the question: Is the significance of the principle the same here as in 
classical dynamics, namely to insure the complete trigonometric stability of 
equilibrium or periodic motion? 

In concluding these introductory remarks, I wish to express the hope 
that these important and interesting questions will receive the further atten- 
tion which they certainly deserve. In the present paper only a beginning 


is made. 


§1. Definition of Stability. 


We propose to consider a system of n ordinary differential equations of 


the first order 
(1) dx;/dt = Zn) 1, n), 


where the right-hand members are real analytic functions of 2, ° --, %n in 
the open connected continuum under consideration. In particular we shall 
study the properties of the solutions of (1) in the neighborhood of an ‘ equi- 
librium point’ of this continuum, i. e., a point for which all of the functions 
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X; vanish. It is especially questions concerning the stability of nearby 
motions that will occcupy our attention. 

There are two quite distinct properties of stability of equilibrium: the 
first is that the nearby motions must continue near to equilibrium for a cer- 
tain time; the second is that these nearby motions have a quasi-periodic char- 
acter. We propose to incorporate both of these conditions in our definition 
of stability. 


Definition of Stability. An equilibrium point of the differential system 
(1), taken at the origin in 2,,- * -, %» space, will be said to be stable of the 
mth order if it possesses the two following properties: (A) all points lying 
at a distance e from the origin at time f remain at a distance not exceeding 
Ke for an interval | t—t,| < Le” (K, L, definite positive constants) ; (B) 
for any fixed interval |¢—t,|<T the coordinates %,---, ¢ of any such 
motion and polynomials with lowest terms of degree s > 0 in the coordinates 
can be represented by trigonometric sums 


Ao +> (Ai cos pit + sin pit) 


of not more than NV +1 terms with an error not exceeding Qe™** (P, N, Q, 
definite positive constants) during the fixed interval. 

If the equilibrium point is stable of all orders it will be said to be com- 
pletely stable. 

It is obvious that stability of any order necessitates stability of all lower 
orders. 

The equations (1) are invariant in form under any point transformation 


vi = $i *, Zn) (t=1,---, n) 


where the functions ¢; are analytic with | 0$;/0%;| 540 at the transformed 
equilibrium point. The order of stability is clearly invariant under such a 
change of variables. In employing these transformations we shall always 
take the transformed point to be at the origin also, so that ¢; are given by 
power series in %,,° * +, Z, without constant terms. 

It will be one of our principal purposes to show the relation between these 
properties (A) and (B), and we shall find that ‘ perturbative stability’? (A) 
does not necessitate ‘trigonometric stability’ (B), though it turns out that 
the requirement of trigonometric stability suffices to insure the stability of 
perturbations. We find it convenient first to go as far as possible on the basis 
of property (A) alone, and then employ the more fundamental require- 
ment (B). 
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It is extremely instructive to consider a few characteristic examples in 
this connection. 
The simplest possible example of a point of equilibrium is furnished by 


the equation 
dz/dt = az, 


in which the point lies at 0. Here the exponential form of solution 
shows that there can not be first order stability even in the sense (A). 
A more instructive example is furnished by the equations 


with equilibrium point at (0,0). The general solution is easily determined 


to be 
To cos(6) + t — ty) To sin (0. + to) 


where 79, 6) are the polar coordinates of (z,y) at t). Thus r, is the distance 
e which enters into the above definition of stability. 

In this case the expressions for 2 and y remain of the order of e for 
|t—t) | S Le? (0< L <1), so that the equilibrium point is stable of the 
second order in the sense (A). But since the denominator in the expressions 
for x and y vanishes for t — tp = e€, it is obvious that the point is not stable 
of the third order. The approximate trigonometric form of the solutions with 
an error of order 3 in e for any fixed interval | t—t) |= 7’ shows that the 
situation is the same for stability in the sense (B). 

On the other hand, for the equations 


dx/dt = y + y(z? + y’), dy/dt = —z—2(2?+ y’), 
the general solution is readily found to be 


cosa(t — ty) —y sin a(t — to) (a=1+%*? + yw’), 
y=, sina(t—to) + y cosa(t—to) . 


Hence the origin is completely stable in this case in the full sense (A), (B). 

All of the preceding examples seem to indicate that trigonometric sta- 
bility is invariably associated with perturbative stability. The following 
example is of the simplest type which illustrates that this is not the case. 
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We shall take n = 4 with two pairs of variables 2, yi, V2, y2, and the equations 


da,/dt = — pry: (22? + yo"), dys/dt = pits + ays + y2"), 
dar,/dt = — + + yr"), = pot, + + yr”), 


with a, 8 of opposite signs. For convenience in dealing with the equations we 


shall introduce conjugate variables 


uw=atyV—1, (i=1, 2), 
and write the equations in the form 


du,/dt = + AUsV2)U;, dv,/dt = (— V—1 + atev2) 
duz/dt = (peV—1 + Buyv,) U2, dvo/dt = (— pe V—1 + Buir:) v2, 


where a, 8 are real constants with a8 < 0. It is to be noted that the equations 
themselves are of conjugate type: that is, when, for instance uw, and v,; are 
interchanged throughout and /—1 changed to —\V—1 in the first pair of 
equations, the equations are unaltered as a set. 

Now multiply the first and second equations of the set just written by 
V1, U, respectively and add; operate similarly with the second set. There is 


obtained the pair of equations 
dw,/dt = 2aw Ws, dw./dt = 2Bw,we, 


where w = U,V, and wz = UzV2 are positive or zero for all real solutions 2%, y;, 
Lo, Yo Of course. 

If w,, for instance, vanishes for any ¢, these equations show that it will 
vanish identically, with w, a constant.* This gives u,—v,—0, and the 
differential equations in Us, v2 enable us to write the corresponding solution 


Uy = 0, V1 0, Us eV-1 Pal t-to) Vo V2 e-V-1 pa(t-to) | 


Consequently for this type of solution the squared distance from the origin, 
measured by w, + we is constant, and the representation is trigonometric. 
But let us consider the general solution at distance « from the origin at 
t = t, so that 
w, w, w, 0, we 0, 


* This follows from the fact that the initial values of w,, w, at time t, determine 
the solution uniquely, in virtue of the elementary existence theorems, for w,=0, w, = 
const. is such a solution in the case at hand. 
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thus laying aside the possibility just considered. Introduce a parameter + 
for which dr = 2w,w.dt. The equations for w,, w2 show that we have 


W, =a(r—70) = B(r—t.) 


so that w, and w, vary monotonically and in opposite senses. 
Because of the integral relation 


Bw, — aw, = Bu, — aw.“ = const., 


in which B and a are of opposite sign by hypothesis, the positive quantities 
W1, W2, and consequently + also, vary in a finite interval as ¢ passes from 
—o to+o. 

Let us observe the relation between ¢ andr further. Taking fp = 7) —0 
for simplicity of notation, we find 


dr 
(ar + wi) (Br + 


Here the lower limit of integration lies between the two zeros of the denomi- 
nator of the integrand. Thus 7 ranges between these values 7’ and 7” only. 

The integral relation above demonstrates that w, and w. remain less than 
a fixed multiple of « for all time. Hence we have stability of every order m, 


in the sense (A). 
On the other hand we have from the original equations 


du,/dt = aw2)u,, dv,/dt = (— pp V—1 04, 
due/dt = (poV—1-+ Bui)u2, dv2/dt = (— poaV—1-+ Bu.) v2, 


whence we obtain by integration 


Uy — Uz erV-1 t+a fusdt, V4 = t+a 
Ue Uy ereV-1 t+B V2 e-P2V-1 t+B 


Consequently the coordinates u;, V1, U2, V2, have amplitudes which vary mono- 
tonically, but they are essentially trigonometric in the angle variables. 

It is obvious (see § 3) that 21, y:, %2, y2 and 4, V1, U2, V2 cannot be repre- 
sentable by trigonometric sums as demanded with an error of the second 
order in e. Consequently property (B) fails to hold for m= 1. 

Inasmuch as we are not following out very systematically the conditions 
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for stability of type (A), it may be remarked at this juncture that in general 
the property (A) will follow if there is an integral in which the initial terms 
give a definite quadratic form, just as they did in the above example. One is 
naturally led to conjecture that, in general, condition (A) can only obtain for 
some order m if there exists a polynomial P in the variables z,,° °°, %» with 
initial terms yielding a homogeneous form of degree & of definite sign, such 
that the series on the right in the equation 


dP/dt > X; 
j=1 


will begin with terms of degree m-+k& at lowest. At any rate the type of 
argument employed in the following paragraph shows readily that this con- 
dition is sufficient for (A).* 

Another example illustrating the possibility that (A) may be true for 
m = 1 without (B) holding for m —1 is furnished by the equations 


dx/dt = xy’, dy/dt = — x*y. 


Here x? + y? const. is an integral, so that (A) holds for any order m. 
The general solution is 


(1 + (1 + ya? 


which is not trigonometric in type. However, here the origin is not an iso- 
lated point of equilibrium, since all of the points of x0 or y=0 are of 
this type. On account of its degenerate character the example is not as illu- 
minating as the preceding one. 


§2. First Order Stability. 


The condition (A) for stability has reference to solutions of the given 
differential equations (1) lying for t= 7¢, at a distance e (less than some 
definite quantity) from the equilibrium point, which will always be taken to 
be at the origin. Now suppose that we write 7; = ey;, thus making a linear 


*In this connection the reader may be referred to Picard, Traité d’Analyse, Vol. 3, 
chaps. 8, 9, where the fundamental papers of Poincaré, Liapounoff and Levi-Civita are 
also cited. 
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transformation of the dependent variables. The equations (1) then take the 


form 
dyi/dt =1/e Xi ° yn) = ° Yn €) (1=1,° 


Since the X;’s vanish at the origin in 2, ° - +, @% space, Y; will be given as 
power series in 41, °° *, Yn, € without constant term. Moreover the domain 
of uniform convergence of these new series will certainly extend beyond a 
distance unity from the origin in Yn space if is sufficiently small. 

Corresponding to any solution of (1) at distance e from the origin to the 
21, °° *, % Space, we have a unique solution y; of the transformed equations 
at distance unity from the origin in y;,- + -, yn space. Let us write yi = ai 
at ¢=7¢, and regard these initial values a; as fixed, while e is a parameter 
which approaches 0. There will be defined a set of solutions y; varying ana- 
lytically with e«, according to the fundamental existence theorem for ordinary 
differential equations containing a parameter. For e=0 the set of functions 
will satisfy the corresponding equations of variation 


dy; n OX; 

@) 
obtained by setting «—0 in the differential equations for y;. These are n 
ordinary linear differential equations with constant coefficients. 

Suppose now that we make the hypothesis that the equilibrium point of 
(1) under consideration is stable in the sense of (A) for m=—1. It follows 
that = is uniformly bounded for | | of order Consequently, 
by allowing e to approach 0, we infer that the limiting functions y; are 
bounded for all ¢. Thus for stability of type (A), m1, it is proved that 
the equations of variation (2) have all such solutions y; bounded for all ¢. 
But the most general solution of (2) is derived from solutions of this type 
by multiplication by a constant factor, since the only restriction made was 
that the sum of the squares a@,,:-*, dn is 1. Hence it is a necessary condi- 
tion for stability of type (A), m =—1, that every solution of the equations of 
variations (2) be bounded for all f. 

The known nature of the solutions of (2) enables us immediately to 
determine the precise significance of this condition. 

Suppose first that the roots a of the characteristic equation 


0X; 


= 0 (8:3 = 0,447; 843 1) 


—adbij 
0 
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are distinct from one another. Then there are n particular solutions of the 
equations of variation which may be written in the form 


Yr = Yn = Cnje™* (j=1,° +9 


The most general solution is obtained by linear combination of these particular 
solutions. Here a;,- ° *, an may be real or complex of course, but correspond- 
ing real solutions may be obtained by combination from the pair of solutions 
corresponding to a complex root and its conjugate. It is apparent that such 
solutions are bounded only in case the expressions e“? are bounded. Hence 
all of the a;’s must be pure imaginary quantities or zero. In this way we 
see that in the general case of distinct roots of the characteristic equation, 
it is necessary that the roots occur in pure imaginary pairs except for possible 
zero roots. 

On the other hand in the more complicated case where the roots of the 
characteristic equation are not distinct, there may arise degenerate solutions 
not of the above type, but of the form 


corresponding to the case in which the elementary divisors corresponding to a 
multiple root aj are not simple. Evidently these are not bounded. In this 
way we are led to the following first conclusion: * 


For first order stability of type (A) at an equilibrium point it is necessary 
that all the roots of the characteristic equation belonging to the point be 
grouped in pure imaginary pairs except for possible zero roots, and that the 
elementary divisors of the corresponding a matrix be distinct. 


We are led at once to ask whether or not this necessary condition for (A) 
is sufficient. It is easily proved to be sufficient as follows. : 

Let us effect the linear transformation of (1) (in general complex) which 
reduces the linear terms in X; to the canonical fotm aiz;. This may always 
be done in the case of distinct elementary divisors. 

In accordance with the well-known theory of real linear transformations 
this may be accomplished in such a way that pairs of variables 2; and 2; are 
conjugate which correspond to conjugate roots a; and a; of the characteristic 
equation. We shall take pains to preserve this conjugacy at every subse- 
quent transformation. It should be observed that for such conjugate vari- 


* The main facts here are of course well known. 
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ables the squared ‘ distance’ from the equilibrium point is properly measured 


by the sum of the squares of the quantities | 2; |. 
The transformed equations take the form in the new variables 


dxi/dt = aia; + Xi (4=—1, n), 


where the power series X; in the variables 7,, °°, 2% begin with terms 
of the second degree, while the a;’s are either pure imaginaries or zero. The 
series X;, Xj corresponding to conjugate variables 2;, x; will be conjugate of 
course, while the series X; corresponding to zero roots of the characteristic 
equation will be real. If we make the further change of dependent variables 


we find the equations for z; to be 
dz; /dt X,(e%*z,,° -, 


Since the exponential factors occurring here have absolute value 1, we have 
clearly for some constant X 


| Xe] SX | |?) = Xe’, 


so long as 41,°* *, Yn OF 21,° * *, 2n are sufficiently small in absolute value. 
Hence we have 


d n d j 


the elementary inequality for real quantities b;,° - -, dn, 


4,j=1 


is employed in deducing this result, as well as the obvious relation 


| < day | sxw. 


Since u>£0 we infer 
—Vn Xdt S duf’? + Vn Xdt, 
whence by integration 


— Vn X(t—t) S 1/m—1/u S Vn X(t—t) (t > tr) 


with a like inequality for t < t. 


— 
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But at ¢==t, we have u—e of course. It appears then that wu will not 
increase to 2e or diminish to $e within an interval 


— <= 
2Vn X 


This is what we desired to prove. 


The stated necessary condition is sufficient for first order stability of 
type (A). 

A slight extension of the preceding proof leads us to the fact that if this 
condition (A) is satisfied, so that the general solution of the equations of 
variation is composed of trigonometric sums, then (B) will hold for m—1, 
i. e., the coordinates can be suitably approximated to by such sums. 

Within the same ¢ interval we have | X; | S 4Xup” so that from the differ- 
ential system in 2,° °°, Zn, the inequalities 


| 24 | = 4X Up" | t — to | 
result by integration, or, in the variables 2;, 
| — aj et (t-te) | = | E— |. 


In other words the n trigonometric sums which constitute the solution of the 
equations of variation with the same initial values as 2; at t, represent 2; 
with an error not exceeding 4Xe?|¢—¢,| during an interval |t—t|S 
«1/2\/nX. A like property must then hold for the solutions of (1) which 
are related to the solutions of the prepared equations (3,) by a linear trans- 
formation with constant coefficients. Hence we may state the following result: 


The stated necessary condition for first order stability is sufficient. If 


then we let 


be n linearly independent solutions of the equations of variation (2) so that 
a; are pure tmaginaries or zero and the determinant | ci; | 0, and tf we 
form the unique linear combination of these solutions taking on the same 
values as the solution of (1) at t = to, the real trigonometric sums so obtained 
are of the form 


Ao + (4; cos pit + B; sin pjt) 


j=1 
where A;, B; are of the first order in e«, and represent 2; with an error not 
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exceeding Ke? |t—to| during an interval |t—t,|S Le! (K, L, definite 


positive constants). Here pi,°° +, ps are the s absolute values of the 2s 


pure imaginary quantities ai, so that 2s Sn. 


§3. Supplementary Results for First Order Stability. 


It was stated in §1 that the condition (B) for stability necessitated the 
condition (A). We propose at first to establish the fact in the case of first 
order stability. In order to do so we prove a lemma concerning trigonometric 
sums of the type referred to in the statement of (B). 

Lemma Concerning Trigonometric Sums. If a sequence of trigono- 


metric sums 
N 
Ao + & (A; cos pjt + B; sin p;t) (|pi—p;|2=P>0) 
j=1 
where P, N are fixed positive constants, approaches a limiting function $(t) 
uniformly in the fixed interval | ¢—t,| <T then ¢ is itself a trigonometric 
sum of the same kind. 


Proof. The lemma may be proved by induction. Let us first establish 
it in the case V = 1, and let us take f) = 0. Here it is given that 


= lim [Ao + cos p, Mt + By™ sin ] 
1=00 


where the limit is approached uniformly by the sequence ¢:(t). 
Now the double integration of the term of trigonometric type here occur- 
ring simply divides it by — p,“*. Hence the sequence of functions 


di(t) dt? + 


consists of polynomials of the second degree in t. 

Suppose now that p,‘” has some limit p, (| =P since po” —0) 
other than © as / becomes infinite. Since the uniform limit of a polynomial 
of any degree can only be a polynomial of that degree, it is necessary that 


dt? + o(t)/p.? 


is a polynomial of the second degree in ¢. Consequently, by double differen- 


tiation, we derive 
d*p/dt? + = const., 


so that ¢ has the form 
= Ao + (A; cos pit + B; sin p,t) 2 P) 


| 
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in agreement with the lemma. 
If, however, p,‘” has only the limit o we conclude immediately by a 


similar method 
= Ap, 
in agreement with the lemma. 

We remark also that if additional polynomials of the second degree are 
added to the trigonometric sums, NV = 1, a similar proof shows that the limit 
¢(t) necessarily consists of such a sum augmented by a polynomial of the 
second degree. 

To prove the lemma for VN = 2 we need only isolate one of the two trigo- 
nometric terms with period 27/p,‘" in the approximating trigonometric 


sums. Proceeding as above we are led to consider 
Sot fot dt? + 


This expression lacks the trigonometric terms of the type isolated so that if 
po” has a finite limit po, 


foi p(t) dt? + 


can be approached uniformly by trigonometric sums of the above type of order 
N =1 augmented by polynomials of the second degree, and so is itself simi- 
larly constituted by the previous remark. By double differentiation, we con- 
clude that ¢ is a trigonometric sum of the second order and of the stated type 
in this case. When p.‘” has only the limit o the method can be carried to 
a conclusion as before. 

Evidently we have here the basis of a method of proof by means of induc- 
tion for N — 1, 2,:--. Hence the lemma must hold for all values of N. 

We may prove that the condition (B) implies (A) for m —1 by making 
use of this lemma. By hypothesis z;/e can be approximated to by a trigo- 
nometric sum of the prescribed type in any fixed interval |t—¢,| <T with 
an error of the first order at most in «. But by choosing 2; properly, zi/e 
can be made to approach any solution y; of the equations of variation having 
a> +:+++a,.7=1. By the lemma then since the functions y; can be 
approached uniformly by such trigonometric sums they must be such sums 
themselves. This means of course that the general solution of the equation of 
variation is trigonometric, which was seen to be a necessary and sufficient 
condition for stability of the first order. 


It follows that the condition (B) alone in the case m =1 implies con- 
dition (A), and thus stability of the first order. 


‘ 
& 
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The second result to which I wish to call attention is implicitly contained 
in the results for trigonometric sums specified in § 2 and has to do with certain 
-amplitudes. Let us formulate the simple notions involved. 

The amplitudes are defined to be the absolute values of the variables 2; 
under consideration, taken as a solution of the prepared equations (3,). If 
there are s pairs of conjugate imaginary roots a; and k& zero roots (n= 
2s +k) then there are only the amplitudes 


C1, Cs; Css15 C sak 


since conjugate imaginary 2;’s have the same absolute value. 
It is a fundamental property of the s+ % amplitudes that the sum 


serves to measure essentially the squared distance of a point in the original 
%1,° * *, &, space from the position of equilibrium. In fact the linear trans- 
formation relating (1) and (3,) evidently takes the above sum into a positive 
definite quadratic form in the values of 2, -°-, tm. Any such quadratic 
form serves to measure the squared distance from the origin as well as the 


sum of the squares used earlier. 
The fact which we desire to establish directly may now be formulated in 


the following way: 


In the case of first order stability the s + k amplitude variables C,,:--, 
Cs.x change only slightly in comparison with the initial total amplitude, 


e= (Cy? (t = to) 


during a time interval | t—t. | S L*e1 where L* is to be suitably chosen. 
For conjugate 2; and 2; the prepared form (3,) gives 


/dt = + xjXi, 
while for real 2; we have 
da;?/dt 24 


Hence if we write u? for the squared amplitude and recall that X; commence 
with terms of at least the second degree, we find, as in § 2, 


du/dt = Xu’, 


where X is an appropriately chosen positive constant. Now in §2 it was 
shown how from such an inequality we may infer that u will vary between 
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ju, and 2u, in a suitably taken interval |t— Le+. The same argu- 
ment permits one to infer that wu will vary between (1—7)u,) and (1 + 9)u, 
y arbitrarily small, if L is chosen small enough. In other words the sum of 
the squares of the individual amplitudes only varies slightly as compared with 
the total initial amplitude under these conditions. 


But the variables x,,- ++, 2% of the prepared equations (3,) are not 
uniquely determined, for we may replace 2%, respectively by 
Cn%m respectively, where Cn are arbitrary positive constants ex- 


cept that c; and c; associated with a conjugate pair x; and 2; are to be equal. 
In so doing the definition of the amplitudes is changed in that these are 
replaced by arbitrary constant multiples of themselves. Thus we conclude 
that, for any arbitrary set of s + & positive constants Ai, °° +, Asx, the sum 


changes only varies slightly as compared with its initial value for an interval 
of time of order e1. 

Now choose s + & linearly independent sets Ai, and an arbitrary 
small value of y >0. It will then be possible to choose for each correspond- 
ing sum and therefore for all, a constant Z* so small that the ratio of the 
corresponding sum of squares to the initial sum wp? at ¢ = ¢, will vary as little 
as may be desired for | t—t, | S L*e?. 

The theory of linear algebraic equations shows then at once that all of the 
amplitudes can vary only slightly as compared with the initial total ampli- 
tude, as first measured. This yields the result stated. 

It will be proved in the following section that when the amplitudes are 
more carefully defined in the case of first order stability, n even, the variation 
in these amplitudes remains slight for an interval of time, | t — fo | S L*e*. 


§ 4. Second Order Stability. 


In dealing with second and higher orders of stability, it is convenient to 
treat separately the cases when the order of the differential system (1) is even 
(n = 2s) and odd (n = 2s +1). Furthermore we have now to take account of 
commensurability relations between the characteristic frequencies p;/2z. 

For the case of even order we confine attention to the case where the 
inequality 

tipi tsp, 0 


a 
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holds for sets of integers i,,--*-, ts not all zero. This inequality excludes 
in particular any zero or equal roots of the characteristic equation. For the 
purpose of the discussion of second order stability here made, only these 
special restrictions are necessary. 

In the case of odd order there is always one zero root which we take 
to belong to the coordinate z,; we designate the 2s other pure imaginary or 
zero roots as in the even order case, and again make the assumption stated 
above, although only the case of further zero roots or of equal roots needs to be 
excluded as far as the argument of this section is concerned. 

It is with the intention of avoiding undue complexity that restriction to 
the ‘ general case ’ is thus made. 

We take up first the case n = 2s, and suppose that pairs of conjugate 


variables are called 


The prepared form of equations (3,) may then be written 


(3) dé,/dt = pi V—-1 + Bi, dyi/dt=— pi V—1 i + Vi 
++, 8) 


where #;, ¥; contain no terms of lower than the second degree. 
Let us attempt to proceed with the preparation of these second degree 
terms by means of a further transformation of the form 


8), 


where both of the functions F;, G; are conjugate homogeneous polynomials 
of the second degree in &,,- - -, 7s, so that if é;, 7; are conjugate imaginaries, 
8) will be, and conversely. The inverse transformation 


has the form 


where F;, G; are the same polynomials in &, 7: as Fi, Gi in €;, 7:1, and where 
higher order terms are only indicated. Such a transformation has already 
been observed not to affect the order of stability ($1). 

On making this change of variables in the prepared equations (3), we 
obtain the specific modified equations 


| 
| 
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= piV—1it+ + 
8 — OF; 
(4) 2 piv 


— 0G: 0G; ) 
piv 0; nj On; + 


where dashes over the €;’s and yi’s have been dropped, and where ;°, ¥;“ 
are the second degree components of ®;, ¥; respectively. The displayed parts 
contain the complete second degree term in the new equations. Hence this 
change of variables does not affect the character of the prepared form. 

If we let c; stand for the complete coefficient of a definite term of the 


second degree 


in ®;, while gq; is the analogous coefficient in F;, then the corresponding coeffi- 
cient in the second degree term of the first equation above is 


— [pi(li— mi) +° °° + (ik 
+ (ls—ms) ps] V—1 4. 


It is at this point that the hypothesis concerning the general case enters. 
The coefficient of gi cannot vanish unless 


Now the argument just stated has not made use of the fact that we are 
dealing with the simple case 1, + m,-+- +--+ m,—2 when there is evi- 
dently no possibility for the coefficient of gi to vanish. Hence we conclude 
that without any assumption of second order stability, we may further pre- 
pare the equations (3) so that no terms #;°, &;° remain. 

We may now prove the following result: 

In the general case of first order stability, n= 2s, there will also be 
second order stability. Let the combined linear and quadratic transformation 
from &i, ni to &:, ni be made which takes (1) into the prepared form (3) lack- 
ing second degree terms in ®, %;. If we substitute 


for i, mi in the formulas, the trigonometric sums obtained will represent 
2 


17 
| 
| | 
| 
| 
| 
| 
| 
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1,° * with an error not exceeding |t—t.| during an interval of 
time |t—t, |S Le®. 

The details of the proof are very much like those for the analogous result 
concerning first order stability (§2). To begin with, we note that in. the 
prepared form the non-linear terms on the right are of at least the third 
degree. 

Now, using exactly the argument of §2 except that we can assume 
| X; | S Xu’, we are able to establish the relation stated. 

Suppose now that we define the amplitudes C,,: - -, Cx. more precisely 
by the aid of this especially prepared equation (3). Following a natural 
modification of the argument made in the preceding section we are at once 
led to infer the following additional result: 


In the case of first, and thus second, order stability (n = 2s) the s indi- 
vidual amplitudes C,,- - -, Cs may be defined by means of prepared equations 
(3) in which ;, ¥; contain no second degree terms. With this definition, these 
amplitudes change only slightly in comparison with the total amplitude during 
a time interval | t—t) | S L*e? 


We turn next to the general case of odd order n = 2s + 1 for the system 
of equations (1), and propose to briefly indicate some analogous results. 
In this case, besides the 2s equations like (3), we have a last equation 


where X begins with second degree terms or higher. We attempt to prepare 
further these equations by a transformation 


&=&4+Fi, Gi, In +H 8), 


where Fi, Gi, H are homogeneous polynomials of the second degree. We find 
precisely as before that we can make #;°? = ¥;°?) —0 by a proper choice of 
F;, Gi, while the second degree terms in X have the form 


' If we take typical corresponding terms of the second on in X® and H, 
both with exponents ms, rin 7s, Un respectively, and with 
coefficients c and q, the new analogous coefficient is 


c—[pi(l,—m) +°* ++ ps(ls— me) ] V—14¢. 


| 
— 0H 0H : 
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In the case at hand the only possibility that the coefficient of q vanishes is 
that some 1; = mi =—1 while the other /;, m; must then vanish. Thus the 
terms é; and 2,” in X® cannot be removed by this method. 

In consequence of this failure in removing certain second degree terms 
it seems probable in advance that for the case n = 2s + 1, first order stability 
does not imply that of the second order. 

A simple example can be devised to show the truth of such a surmise. 
Take n = 3 and write the equations 


dé,/dt = /dt = (— V— + ¢x,")m, 
dx,/dt = (ce 0), 


containing one such non-removable term in the last equation. 
On writing u = é,m, v = 23 so that up + 1%” =e’, we obtain immediately 


du/dt = 2cv7u, dv/dt = eu, 
whence 
d*v/dt? = 2cv? dv/dt. 


This is immediately integrable, and yields 


Suppose now that we take vo = 0 and substitute 


ew 
t — typ = (3/2c) f 
0 


We find 
dw 


w* +1 

Taking the upper limit w—-+ o we conclude that the corresponding 
interval of time is of order e~*/*, since the integral converges. Thus v does 
not remain finite in an interval of order e~*/*, and property (A) does not hold 
for m2. Hence in this special example the presence of non-removable 
terms is associated with instability of the second order. 

We proceed now to prove that (B) can only hold if there are no such 
nor.-removable terms in the completely prepared form. Here we use once 
more the lemma on trigonometric sums of the preceding section. 

Let 


i 

| | 

| 

| 

| | 

| | 

| 

i 

| | 

| | 
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be that solution of the prepared equations for which 

n= In = € = fo). 
Here a; and 6; are conjugate imaginary quantities, and c“ real, with 


Also are solutions of the equations of variation 


é;* ePiV-1(t-to) , ni* e-PiV-1(t-to) = c, 


We consider now the second terms in these series. 
By differentiating the prepared equations twice as to e and setting «= 0, 


we find 


where the terms on the right in the last equation arise from the non-removable 
terms of X. Furthermore for t= fo, é;** = = 2,** = 0, since, for ex- 
ample, é;** is the second partial derivative of &; as to e, and &; is constantly 
equal to aie at t = bo. 

Thus we infer that &;**, 4:** are identically zero, while from the last 


equation 


But 2, can be approximated by a trigonometric sum of the stated type to the 
third order in e, by hypothesis. Thus (a, — ce) /e? can be similarly approxi- 
mated to the first order in e. Consequently the limit of this quotient, x,**, 
is such a sum by the lemma. This can only be a trigonometric sum for all 
choices of ai, 05, if =d=0 (t—1,- 8), which we wished to 
prove. 


In the general case n= 2s +1 a necessary and sufficient condition for 
second order stability is that a prepared form (3) exists with a supplementary 


equation 
dx,/dt X, 


where ®;, Wi, X begin with terms of at least the third degree. Thence 
follow results concerning the representation by trigonometric sums, and con- 
cerning the variation of the s+1 suitably defined amplitudes, which are 
entirely analogous to the results in the case n = 2s. 


| 
dé;** a 
a piV—1é:**, = — piV—17n:**, 
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In fact, we have reached a prepared form in the case n = 2s + 1 with the 
same characteristics as that obtained in the case n= 2s earlier, and can 
apply the same arguments with obvious modifications. 


§5. Third and Higher Order Stability. 


We propose now to assume stability of the third order, taking n = 2s 
first, and then, by making use of condition (B) only, to obtain a necessary 
prepared form of equations, which can readily be proved sufficient for stabil- 
ity of the third order. An entirely analogous treatment can be made for 
n= 2s-+1. 

By passing thus from the case of second to third order stability as well 
as from the first to the second, the general results can be formulated. 

Starting with the prepared form (3) we now attempt to remove as many 
of the third order terms as possible by writing 


8) 


where now Fi, Gi are homogeneous polynomials of the third degree in &, 

‘, 7s» By exactly the same formal work as in our treatment of second 
order stability, it turns out that all such terms may be removed except for 
terms of the form U;&;, Vini in the ith pair of equations, where U;, V; are 
linear combinations of the s products €sys. Thus we have 


dé;/dt = 
dyi/dt = (—~av—1 Vi) ni +: s), 


where the terms not explicitly written are at least of the fourth degree. From 
these equations we derive 


d (ini) 


(Us + Vaden 


where the omitted terms are at least of the fifth degree. 

Now we expand é;, i in power series in e as before. To determine the 
coefficients in this series we proceed precisely as in the case of second order 
stability and find series in e which we write fori—1,---,s, 


| 
j 
| 
{ 
| 
| 
Be 
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In virtue of the hypothesis of third order stability and the use of our lemma 
on trigonometric sums, we conclude that ¢i*, where 


are trigonometric sums. 


If now we substitute these different series in the equations above in qi, 


and compare terms in e* we find simply 
= [Ui (a,b, +5) + Vi (arb, +) 


Consequently unless U; + V; vanishes identically for i—1,---, s, the 
functions ¢;* will not be trigonometric sums but will be linear functions of ¢, 


which is not possible. 
This gives us the necessary form of prepared equations in the case of 


third order stability, n = 2s, 


(5) + Gj, 
dyi/dt =— Mini + Vi, (t=_=1,---, 8), 


where M; has a constant term p;\/—1 and linear terms in the s products 
€sys, and where begin with of at least fourth degree. 

Our first general conclusion, based on a step-by-step process to higher 
and higher orders, is thus the following. 


For n = 2s a necessary condition that (B) may hold for any odd order 
m of stability ts that the equations may be transformed to the form (5) in 
which the functions M; are polynomials in the s products &,m,-° °°, &sns of 
degree 4(m—1) at most, with initial terms V~—1 i, and having pure 
imaginary coefficients throughout,* and where ®;, ¥; commence with terms of 
degree at least m-+-1. In this case Bi, ¥; may be made to commence with 
terms of degree m + 2 at lowest without further assumption. 


To show that this form (5) induces (m+ 1)st order stability (m odd) 
we obtain first of course 


Sem) |= < Xyms, 


where = +° + ésys and thence 


< 


“These are pure imaginary quantities since the variables £;, 7; are self-conjugate, 
and M,é; must be conjugate to — M,7;. 


| 

| 

| 

| 
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Here X is a fixed constant. From this inequality we infer that w remains of 
the first order in ¢ for the interval | t—% |S Le“"*? by the same method 
already employed for m= 1. 
This gives us the property (A) as a consequence of the prepared form. 
Now we deduce also from the same inequality 


| — | S | to | 


in this same interval. On going back to the prepared form this yields at 


once 
| dé,/dt — | Kes | t— ty | + Ne? 


in which M; designates what 1; becomes when é;, 7: are replaced by &, 
Hence we infer 


| d (t-te) /dt | | t—t, | + Nem, 
and finally 


| — EMO | = to)? + Nem? | ty | 
in an interval of the order of e+), 

Hence we infer in passing that property (B) holds as well, and obtain 
additional properties of the approximate trigonometric representation. 


The necessary condition of the prepared form (5) insures stability in the 
complete sense that (A) and (B) hold, and so is sufficient. The explicit 
solution of the prepared form with terms ®;, %; omitted, and with initial 


values &, chosen to correspond to the initial values -,2n of 
t1,°* *, tn yields trigonometric sums by means of the composite transforma- 
tion which represent x; with an error not exceeding 
t — to)? 
Kems + Nem? | t— ty | 


during a time interval 
| t— | <= Lem) 


where « denotes the initial distance from the origin and where K, L, N are 
positive constants. 


If we define the amplitudes C,,--~-, Cs of a solution by means of the 
prepared form it is clear that we can establish the following result: 


The amplitudes C,,-+:*, Cs vary slightly in comparison with the total 


| 
| 
| 
| 
| 
| | 
| 
| 
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initial amplitude ¢ during an interval | t—t) |S L*e™, L* being a suit- 
ably chosen positive constant. 


There remains now to refer briefly to the case n odd. It is only the 


handling of the last equation in 2, that needs to be referred to. In passing 
from the second order case to the third it is seen at once that the third order 


terms in this equation may be entirely removed without any assumption con- 


cerning third order stability. The third order terms in the other equations 
can be treated as in the case n even. 


In the case of mth order stability, n = 2s +1, the prepared form is 
(6) dé; /dt = M;é; + i, dyni/dt =— Minit Vi, 
dx, /dt =X, 


where M; are functions of the products Ein; and of x, of degree m at most in 
°° Um, and X have no terms of degree lower than m +1. 
From this form, the properties (A) and (B), suitable trigonometric sums, and 
the slow variation of the amplitudes follow just as in the case n = 2s. 


§6. Formal Theory of Complete Stability. 


In this section we shall consider equations (1) of even order n = 2s, 
and for convenience we shall write the equations in the form 


(7) dé;/dt = 9; (é,,° dyi/dt = qe) 


where the conjugate power series ®;, Y; are without constant terms and where 
the first degree terms in ©; and ¥; are pj\V— 1 and — pi V— 1 respectively. 
We suppose furthermore that we are confronted by the general case as de- 


fined in § 3. 
In other words, of all possible equations (1) we are dealing with those 


typical of first order stability. 
It is obvious that any formal power series transformation 


(8) & + Fi, n= Hi + Gi 8), 


with F; and G; conjugate series in &;, 7, carries the system into another 


formal system of the same type, in which 4;, ¥; may not be convergent series, 
if Fi, G; are not. It is apparent that we have here a formal group leaving 
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(7) invariant, which will be perfectly consistent with itself in the sense that 
thc result of carrying out two successive transformations is the same as that 
derived by the composite transformation.* 

Now as already proved we can remove the second degree terms in %%, ¥; 
by using F;, Gi which are homogeneous quadratic polynomials. Following 
this with a transformation in which F;, G; are cubic, we can remove all terms 
of the third degree save those in ©; and ; of the respective types Uié; and 
Vié; where U;, V; are linear in the s variables &:y:,°- +, €sys. All of this 
goes just as in our discussion of stability in the earlier sections, and we 


arrive at the following first conclusion by an infinite succession of steps. 
Under the formal group the equations (7) may be given the form 

(9) dé;/dt =Uiéi,  dyi/dt = Vini (¢==1,°--, 8), 

where U;, Vi are power series in the s products é:m,° * *, Esns- 


This result is extremely interesting from the formal point of view for the 


following reason: 
If we write ui = Eini then we have the formal equations in u;,° * +, Us 


thus yielding an associated formal differential system of one half the order of 
the original system (7). 


We propose next to ask what is the most general transformation of the 
group (8) carrying the normal form (9) into the normal form. Let us 
consider terms of the second order in F;, G; which alone can modify those of 
the second order in U;, V;. According to the discussion carried out in the 
treatment of second order stability, these are uniquely determined by the 
initial and final forms of (9). But in neither are there terms of the second 
degree. Hence Fi, G; lack second degree terms. 

Pass next to third degree terms in U;, Vi which are only affected by the 
terms of equal or less degree in F;, Gi. Here we find as in dealing with 
third order stability that since U;é; and Vini have only terms of the third 
degree which involve a factor é; and 7; together with a product &;n; linearly, 
Fi, Gi contain only terms of the third degree of the same character. Thus 
we may continue indefinitely, and we infer 


*It is to be observed this group admits of a slight extension in that the initial 
terms é,,7; in the series for £,,7, might be replaced by ¢;£;, djn, respectively (c,d, being 
self-conjugate of course). 


| 
| 
| 
| 
| 
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The most general transformation (8) which carries the normal form into 


itself has the form 
(11) & = &Fi*, 8) 


where F;*, Gi* are series in the products &:91,° * *, &sns with constant terms 1. 


The form of this transformation indicates that the family of formal ‘ surfaces’ 
Q(u,° *, Us) = const. is invariant under the group. 


The general invariant theory of equations (7) is less simple than that 
of the special class of equations possessing complete formal stability. 


A necessary and sufficient condition for complete stability of a system of 
equations (7%) is that U; + Vi=0 (1—1,--: +, 8) so that the equations take 


a normal form 
(12) dé;/dt = Miéi, dni/dt = — Mini (+— 1, s) 


where M,,- ++, Ms are pure imaginary formal power series in the s products 


+, with initial terms +, psV—1 respectively. 


Eini, * 

For it is precisely a result of the earlier sections that stability of the 
mth order yields this normal form to the terms of order m. Also the new 
transformations employed at successive stages (m—1, 2, 3, --~°) do not 
affect terms of lesser degree. Hence we obtain a limiting formal transforma- 
tion which takes the equations into the given form. Conversely, if there is 
such a form, the actual transformation obtained from the formal transforma- 
tion by breaking off at the term of (m -+ 1)st degree will yield the same form 
of equations out to terms of the same degree so that the prepared form is 
actually obtained, and thus stability of any arbitrary order m is ensured. 

These series M; must be pure imaginary since the products iy; involved 
are self-conjugate. 

For this system (12) we have the integrals é:7; const. Hence in the 
most general transformation (11) which keeps the normal form (12) unal- 
tered we have &7: = const. also. On this account it is possible to determine 


the transformed equations at once. We find 


with like equations in 7. 


The most general transformation (11) leaves the normal equations (12) 


| 

dé; 1 d& 

| 
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unaffected except that in each M;, the products are replaced by 
Thus the subgroup for which G;* =1/F;* leaves the equations unaffected in 


every respect. 


The question naturally presents itself as to what are the invariants of the 
equations of completely stable type. The reply is obviously the following: 


Reduce the given equations with the property of complete stability to the 
normal form (12), thus obtaining s associated series 


The invariants of the completely stable type of equations will be the invariants 
of these forms under arbitrary transformations 


where the fi are arbitrary real power series except that they must have a 
constant term 1. 


In this normal form the equations can be at once integrated precisely 
because = give s integrals for i—1,-° -, s. 


The integrated form of (12) is 


When these expressions are substituted in the formal series giving 2%1,° °°, In 
in terms of &, °° +, ns the general formal trigonometric solution is obtained, 
involving &,- ++, 9s as the n arbitrary constants of integration. 

The products give real power series in the s variables 


evidently beginning with a quadratic form which is the product of two con- 
jugate linear factors. This result may be stated as follows: 


If the equations (1) have the property of complete stability there exist s 
first integrals = c; (t=1,-- -, 8), given by formal series beginning with 
a second degree term which ts the product of two conjugate linear factors, 
forming in all n = 2s linearly independent linear expressions. 


Of course an arbitrary Q(u:,- - -, us) = const. also yields a first integral, 
for any real power series Q. 

It remains to say a word about the case of odd order n= 2s+1. It is 
the normal form that is particularly to be noted. This may be written 
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(127) dé,/dt = dnyi/dt = Mini, dxn/dt = 0 
Here M; is a power series in the s products £m, °° *, €s and in @» only. 
Thus the equations separate into a system of the type obtained in the case of 
even order together with an equation yielding x, — const. for the remaining 


variable. 


iva) 


The normal form of equation for the completely stable case, n= 2s +1, 
consists of the equations (12’). These equations display the slight difference 
between the odd and even order case. In the odd order case there exists a 
real formal series integral ¢=const. (corresponding to an = const.) con- 
taining a first degree term, as well as s integrals $i = const. 


§ 7. Variational Principles and the Equations of Dynamics. 


Let Pi, +++, Pn, Q be n+1 real analytic functions (n= 2s) of a, 
‘+, Z, in a certain open connected continuum under consideration, and let 
C* be an are of a curve with continuously turning tangent lying within this 
continuum given by 
== 
for +, where z;* are continuous functions of ¢ with continuous 


derivatives. In this case the integral 


has a certain numerical value 7* along the curve C*. Now suppose that 


for i—1,- - :,m represents any arbitrary one-parameter family of such curves 
which gives C* for 40, and has the same end points as C* for all values 
of A. Suppose furthermore that the functions x;*(t,A) are restricted to have 
continuous partial derivatives as to t and A. A family 


will fulfill these conditions if the 8x; are continuous together with their 
derivatives and vanish at ¢ = t, and ¢ = ¢,, while A is taken sufficiently small. 
Under these circumstances J becomes a function J(A) when 2;*(t, A) are sub- 
stituted for 2; under the sign of integration, and this function 7() will have 
a continuous derivative as to A whose value for 40 we denote by 8J._ This 


= +A 1,-°-+, 2) | 
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variation 82 may vanish identically along a curve C*, i. e., 8J* 0, and it 
turns out by direct inspection of the expression for 87 that the condition is 
precisely that the coordinates 2,*, ---:, %* satisfy the set of n Pfaffin 


differential equations, 


j=1 0x; 02; dt 02; 


(15) == 

All this amounts merely to the statement of some elementary principles 
in the calculus of variations as applied to a particular and very important 
type of integral. 

On account of the fact that the vanishing of 82 along a curve in this 
sense is obviously invariant of the particular coordinate system, it follows that 
if we make a change of variables 


(16) = $i (Zi, En) ((=—1, n), 


where the correspondence between °°, and %,°* *, Zn space is one- 
to-one and analytic along the curve, the new equations can be obtained by 
determining P,, +--+, Pn, Q by substitution under the integral sign, and 
then writing down the new corresponding Pfaffian equations. But the formal 
relation so established can of course be stated and proved directly, without 
intervention of the principle of variation. 

Suppose now that we consider formal (i. e. convergent or divergent) 
series for P;, Q in the neighborhood of the origin in 2, - + +, Ym space, and 
also the corresponding formal equations (15), and suppose that we effect the 
formal transformation (16) in which ¢; are given by similar series without 
constant terms for which the formal Jacobian has a constant term not zero. 
The new equations can be obtained by the same rule of transformation as 
before, and in general the result of more than one such transformation is of 
course the same as that for the composite formal transformation. These 
statements can be proved by first taking the actual equations obtained by 
breaking off the series for P;, Q at terms of high degree m, while actual trans- 
formations are obtained by breaking off the series for ¢; in the same way. 
By considering the relations involved as m becomes infinite, the stated formal 
results follow at once. 

I state these facts at some length since it is important to realize that in 
this way the condition 60 yields a perfectly definite system of formal 
differential equations (15) subject to the same formal rules as in the case 
of actual convergence. 
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This point of view admits of a slight extension which will come into play 
in our first theorem. Suppose that Pi, Q are generalized to the respective 
forms P;/D, Q/D, where Pi, Q, D are formal power series, but where D need 
not contain a constant term. However, nevertheless we obtain an extended 


Pfaffian system 


i] 
{ 
aQ aD 
Vi OX; 


The equations obtained will be independent of any formal factor » intro- 
duced in all of the numerators and denominators of Pi/D, Q/D. Thus the 
formal significance of the condition J —0 is maintained even in this case. 


If the equations (1) are of completely stable type at the origin, and if 
&i, ni, Mi 8) are the sertes in %4,°* *, which appear in the 
normal form (12), then the equations (1) may be expressed in the extended 


Pfaffian form 


According to what has been said above, it is only necessary to verify that 
(17) leads to the equations (12) when é,,- - -, ys are taken as the dependent 
variables. This may be immediately verified if the conditions dM; —0 
(t=1,--°-, 8) are equivalent to dé; —0, as we shall take them to be. 
A still more formal attack is to change variables in (12), writing ui = Mi, 
vi = log 8). The condition (17) becomes 


(17”) Surat usar } 


' giving the desired equations 
dui/dt = 0, dvi/dt = 2ui, (+=—1, s) 


if we grant the legitimacy of the formal processes involved. 
The equations given by (17’) are of Hamiltonian form with u;, v; as 


8 
conjugate variables and with principal function H = > u,;”. Or, better still, 
j=1 


we can remain within the non-extended type of Hamiltonian equations in case 
M;, aside from its constant term which does not effectively enter in (17), 
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is divisible by i. In this case (Mj; — pj V—1) /€jn; is an admissible power 
series, and (17) may be modified to the equivalent form 


th 


j=1 
giving equations of the Hamiltonian form, if we take 


EimV—1 


as conjugate variables and 13 M;?/ V—1as principal function. Also, this 
j=1 

choice of variables will lie in the original formal group in case M; contains a 

term ci€ini, ci ~0, of the first degree. Of course it can contain no other 


such terms in this case. These results may be summarized as follows. 


The equations (1) of completely stable type can be made to assume the 
Hamiltonian form, at least if we admit the legitimacy of certain extended for- 
mal transformations. In any case if Mi—piV—1 is divisible by ini for 
t=1,-:-, 8, with linear term cigini (ci #0) then the equations (1) can be 
made to take Hamiltonian form by means of (18), without any use of such 
extended transformations. 


The last part of this result is obvious from the normal form (12) which 
becomes of the very special type 


dés/dt = (pi + Emi) V—1 = (— pi — Ems) V— 1 0 
(t=—1,-°-, 8) 
if use be made of the special variables (18). 
Consider now a Hamiltonian system (1) of general type with equilibrium 
point at the origin, which can be written 


du: 0H = dy 


in which the lowest degree terms in H are 


where the numbers p; arise in the same way as for any equilibrium point with 
first order stability. The variables u;, v; are conjugate imaginary variables 
of course so that a preliminary linear transformation of the real variables is 
presupposed. 
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It is well-known that the transformation for arbitrary K, 


=—_ 


preserves the Hamiltonian form, i. e. that the transformed equations are 
obtained by replacing ui, vi by their series in @, 0; in H to obtain H. We 
shall assume that 


where the terms not explicitly written are of higher than the second degree. 
If K has only the first term, the transformation becomes the identity 


Vi=Vi (t—=1,- 8) 


These transformations form a group. 

We propose first to take K of the third degree and thus get rid of the 
third degree terms in H; next to take K as the sum of the first quadratic 
term and a homogeneous fourth degree term and thus eliminate all the fourth 


degree terms of H except those quadratic in the s products w,1, °° *, UsUs; 
and so on. Thus by an indefinite sequence of such transformations we pro- 
pose to make a function of the s products °° Usvs. 


Now the first of these transformations is 


in which the variables on the right of each equation are @;, vi, and K* is a 
homogeneous cubic in these variables. If we express wi, v; in terms of i, Di 


they take the form : 


ak* 
i i + Ov; + i i Oil; 
where K* is K* with the variables %;, v; replaced by a, Ti, and the complete 
first and second degree terms are explicitly written. 
Consequently H‘* as obtained by the usual expansion in Taylor’s for- 


mula is 


8 
piv 1 (w du; Vj dv; 


where H“ designates the third degree terms in H, and the dashes over the 


letters are omitted. 
Evidently then the homogeneous cubic polynomial K* can be so deter- 
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Ov; 0; 


ree. 


ete 


he 


BirkHorF: Stability and the Equations of Dynamics. 33 


mined, and uniquely, that all of the third degree terms disappear in H. In 


fact a term in K* 


leads to a similar new term in H with coefficient 


j= 
where the term in brackets is not zero for 


This is in accordance with our statement since there can be no third degree 
term in the products * +, UsVs. 

Another transformation will lead to a modified fourth degree term in H 
of the same type, 

du; 

where K* is now the arbitrary homogeneous polynomial of the fourth degree 
at our disposal. By the aid of it we can eliminate all the fourth degree 
terms in H“ except those for which 1; = mi, i. e. the very terms we do not 
attempt to remove. 

By continuing in this way, we are led to the conclusion: 


For the general Hamiltonian equations (1) with first order stability, H 


can be taken in the form H(uiv1,°--, Usvs). Inasmuch as the Hamiltonian 
equations then take the form 
dui = Ui dvi 0H ( UiVi) 


of the completely prepared type (12), it is clear that first order stability 
necessitates complete stability in the case of the Hamiltonian equations. 


Of course it is well known that the Hamiltonian equations admit a com- 
plete formal trigonometric solution in the case of first order stability. 

We can now readily attack the question as to the possibility of reducing 
a completely stable system (1) to Hamiltonian form by means of a trans- 


formation belonging to the restricted formal group. 

It was proved in the preceding section that the most general transforma- 
tion of the formal group which kept preserved the completely prepared form 
(12) was given by (11). It was furthermore seen that M,,---, M, are not 
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altered, i.e. that M; is obtained directly from M; by substitution. But such 
a substitution changes the variables in Ms; from to 
-, 7%, in accordance with the formulas 


i*Gi*. 


Ti, ° 


‘ Imagine now the particular transformation performed which reduces the 
completely stable prepared equations (1) to Hamiltonian form, if that be 
possible. Since M,,--:-, Ms; are not altered, there must obtain 

M, = 0H /07i, . 
or, more briefly, 


dH = 


4 
A necessary condition that such a transformation to Hamiltonian form 


be possible is therefore that 
8 
2 M;(m, ms) dm; 
j=1 


becomes an exact differential for suitable real power series 
mi = mi fi (m1, 78) 8) 
where the constant terms c; in f; are not zero. 
In particular this requires that the first degree terms in 

8 

= Cj M 

j=1 
yield an exact differential. More precisely if this term in JM; is 

+ 

we must have 


Cj 
— = 1 q =] S 
(i, j=1,°°-,8) 


which is not in general true for s > 2. 
The special case first pointed out when transformation to Hamiltonian 


form is possible, namely when M;— piV¥—1 is divisible by ini, falls under 
the type treated with dr; dM;i/V—1. 

The necessary condition stated is also sufficient for transformation to 
Hamiltonian form. 

In fact choose F;*, Gi* in (11) so that 7 —&iqji, which is obviously 
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possible in an infinite variety of ways. Then M,,---, M; are unaltered, and 
also 
M, dz, + + M dH. 
It follows that 
0H 
which is what we desire to prove. 

It is apparent in this way that while completely stable systems are of 
equal generality with the Hamiltonian type in a very wide formal sense, they 
are more general if only transformations of restricted type are admitted. 

We turn now to the consideration of the general Pfaffian type of equation 
derived from the non-specialized variation problem. 

At the equilibrium point the equations of variation of a Pfaffian system 
(15) are evidently of the form 

dx; 

= ai j + ) =0 
where aij is skew-symmetric (i. e. aij sa origi and Bi; is symmetric (i. e. 
Bij = Bji). These are linear differential equations with constant coefficients, 
the nature of whose solutions depends on the roots of the determinant equation 
in » 


| + Bis | =0. 


In the case of stability the roots must be zero or pure imaginary, of course, 
while the elementary divisions of this A matrix must be simple. We shall 
assume that there are no linear homogeneous combinations of these roots with 
integral coefficients which vanish. This is the general case. 

In order to obtain a preliminary prepared form for the Pfaffian equa- 
tions, we note first that the constant terms in P; and Q can be omitted since 
they do not enter into the differential equations. Moreover, if the origin is 
an equilibrium point, there are no linear terms in Q. According to the prin- 
ciple of variation, we may directly transform the part of the integral I con- 
taining the linear terms so as to obtain the linear terms with the new vari- 
ables. But it is known that by linear transformation and omission of an 
exact differential, the linear terms in P; and P; become 2; and —2;, where 
x; and z; are conjugate variables in the sense of the Hamiltonian equations, 
while simultaneously Q takes the form 


Q= p(t? +277) 


| 

n 

ly 


36 BirkHorF: Stability and the Equations of Dynamvics. 


where only the lowest degree terms are written. Now, if we note that — x;dz; 
differs from z;dx; by a perfect differential we can write the lowest degree terms 


under the integral sign as 
2 > pj (xi? + 2;?) dt 
where the index 1 runs over s values while 7 takes the conjugate index values. 
But these are essentially the first order terms in the Hamiltonian integral. 
At this stage we may introduce the conjugate variables in; by setting 


where we find the linear terms above take essentially the form 


while 


It is this preparation of the Pfaffian system in its lowest degree terms 
that we shall assume to be made, and we shall seek to reduce it further to 
complete Hamiltonian form by a succession of changes of variables 


1,°° +, 8) 


where Fi, Gi are conjugate series in £,, °° -, js beginning with terms of at 
least the second degree. 
Let us write the prepared linear differential form as 


R,dé, ++ + ++ + Sidm +: Sedns + Qadt, 


and let us seek to eliminate the terms of second degree in f,,---, Rs by a 
change of variables of the above type in which F;, G; are arbitrary homo- 


geneous quadratic polynomials. 
We find immediately that the above linear form is preserved in character 


while we have 


8 0G; 
j=1 


if Ri, Ri denote the second degree terms in R; and R; respectively. Evi- 
dently, by the subtraction of an exact differential, we can replace the sums 
on the right by the more abbreviated expressions — Gj. 
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Hence there exists a set of G; which removes the remaining terms in 
R;, and it is clear at the same time that in this way higher and higher 
degree terms in fj, S; can be disposed of by taking Fi, G; as homogeneous 
polynomials of corresponding higher and higher degree. Thus a formal limit- 
ing composite transformation is set up which reduces R,,- - -, Rs to zero, while 
taking +, Ss into é,,° +, respectively, as desired. 


In the general case of first order stability the Pfaffian equations can be 
transformed formally to Hamiltonian form. 


§8. Questions of Convergence. 


The methods which we have been following are based on the use of certain 
series. Their convergence or divergence is immaterial for the properties which 
we have so far discussed. Nevertheless the questions of convergence are of 
great importance for the theoretic treatment of a given dynamical problem. 

For example, the Hamiltonian and Pfaffian types of equations (1) admit 
of the convergent series integral H — const. and P =const. respectively, 
whence in the case of first order stability we can immediately infer the per- 
manent complete stability in the sense (A). Or again, there may exist a 
convergent series P for which 


0/dx, +: -+0/day (PXn) =0 


in which case f Pdv is an invariant volume integral (dv, element of volume) ; 
from the existence of this invariant volume integral, various important prop- 
erties of recurrence can be deduced. 

In the most favorable circumstances the normal form of the equations 
(12) can be obtained by ordinary transformations defined by convergent 
series. Here the series M; and the trigonometric series solutions of the ori- 
ginal equations will converge. The precise character of the motion can then 
be specified for all time. Concerning this integrable case we prove the follow- 
ing simple result: 


For integrability of the system (1) tt suffices that the formal integrals 
oi = = const. (t=1,-°-, 8) 
are given by convergent series. 


In fact suppose that these real series do converge. Since they are form- 
ally factorable into conjugate factors &i:, the separate components are given 
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by convergent series of course. This follows from Weierstrass’s preparation 


theorem. Hence 
1 dé; 


— = 1/§ > 


are series given as a sum of formal quotients of two convergent series and so 


themselves are convergent. This demonstrates the complete integrability of 


(1) in case the series converge. 
The existence of ‘small divisors’ in the series §;, 7; shows that in general 


there will be no such convergent series ¢;, and the series M; appear to be 
divergent for the same reason. 

The possibilities of classification of completely stable systems (1) on 
the basis of diverse convergence properties are limitless. Every system of 
series or of differential equations which may be attached to (1) so as to have 
invariantive significance under the general formal group evidently affords a 
means of classification of systems (1) in which convergence or divergence 
plays a role. 

It is with these few remarks that we terminate our reference to converg- 
ence questions. While on the purely formal side any completely stable system 
of equations (1) is to be regarded as equivalent to the extended Pfaffian or 
Hamiltonian type derived from a variational principle, it is clear that these 
more special types possess at least one ordinary integral and other convergence 
properties. In this non-formal direction the complete characterization of the 
completely stable equations (1) which may be reduced to the classical types 


still remains to be made. 


Ry? 


Quaternary Quadratic Forms Representing all 
Integers. 


By L. E. Dickson. 


1. For the case in which only squares of the four variables occur, the 
forms representing all positive integers were determined by Ramanujan.* It 
is a much more difficult problem to determine the forms which involve also 
products of the variables whose coefficients may be odd. This paper develops 
a method of solving this problem completely. 

The fact that one of these new forms represents all positive integers is 
equivalent to the fact that a related form involving only squares shall repre- 
sent all positive multiples of a fixed integer, such as 4, 8, 12, 28, 44. A very 
few of the latter facts were known, and their importance has been emphasized 
in the literature. 

The first part of the paper proves the necessary auxiliary theorems on 
ternary forms. They are useful also in other problems in the theory of 
numbers. 


Part 1. REPRESENTATION BY TERNARY QUADRATIC ForRMS. 


2. Dirichlet ¢ gave the simplest proof of the fact that z*+ y* + 2? 
represents all positive integers except those of the form 4*(8n+ 7), but 
represents no integer of that form. This will be expressed briefly as 


(1) +y*?+2? represents all ~4*(8n + 7). 
From (1) we readily deduce { 
(2) + 2y?-+ 22? represents all 4*(8n + 7) ; 


(3) + 2y? + represents all 4*(16n-+ 14). 


* Proc. Cambridge Phil. Soc., 19, 1916-9, pp. 11-21. Proofs of the needed theorems 
on ternary forms are given in the writer’s paper cited in § 2. 

t Journal fiir Mathematik, Vol. 40 (1850), p. 228; Werke, II, p. 91. 

¢ Dickson, Bulletin of the American Mathematical Society, Vol. 33 (1927). 


389 


| 
L 
> 
| | 


40 Dickson: Quaternary Quadratic Forms Representing all Integers. 


The proof is similar to that now given for 


(4) 


If f is even, x= 2X, f= 2F, F + y*? + 422. Applying (3), we con- 
clude that f represents all even integers except those of the form 4*(8n + 7). 
If f is odd, x? = 1, 2y?=0 or 2 (mod 8), whence f= 1 or 3 (mod 8). But 
when m =1 or 3 (mod 8), (2) shows that m is represented by x? + 2y? + 227. 
Then 2? =1, 2(y? + z?) =0 or 2 (mod 8), whence y and z are not both odd. 
Permuting them if necessary, we may assume that z—2Z, whence m= 
x? 2y? + 

The writer proved (loc. cit.) 


f=2? + 2y?-+ 82? represents all ~ 8n + 5 and £4*(8n-+ 7). 


(5) 32? represents all ~9*(9n + 6). 


If g represents 3m, x and y are multiples of 3, = 3X, y= 3Y, and m 
is represented by 3X? -++ 3Y?-+ 2?. The converse is true since the product of 
the latter by 3 is of the form g. 


(6) 32? + 3y?+ 2 represents all + 9*(3n + 2). 


3. We shall need the following results: 

(7) =2+ay+y*?+2? represents all ~9*(9n + 6) ; 
(8) 32? represents all ~ 9*(3n + 2); 
(9) + + 227 + zz represents all 49*(49n + %e), e = 3, 5,6; 
(10) + 2y’ + — represents all 49*(7n + €), e =3, 5,6. 

By the Hessian H of a quadratic form whose cross products have even 
coefficients we shall mean the determinant of the halves of its second partial 
derivatives. By Hisenstein’s* table of improperly primitive reduced positive 
ternary quadratic forms of Hessian = 100, 2A, 2p, 2p, 20 are the only reduced 


forms with all coefficients even of Hessians 6, 18, 14, 98, respectively. 
We employ the form, lacking a term in zy, 


(11) + 2By? + + %Ryz + 


* Journal fiir Mathematik, Vol. 41 (1851), pp. 186-9. 
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Its Hessian will be 2D if 
(12) B= AA—D, A= 4BC — R’. 


Lemmal. If p=3 or’ and A is any positwe integer not divisible by p, 
2A is represented by the unique reduced form 2d or 2p of Hessian 2p with all 
coefficients even. 


Here D—p. We take A—4d, d=2pt-+1, and choose ¢ so that 
B=4dA — p= 8Apt + 4A —p isa prime. This choice is possible in infi- 
nitely many ways, since the coefficients of the linear function of ¢ are rela- 
tively prime. Since B=— p= 1 (mod 4), 


(—A/B) = (d/B) = (B/d) = (—p/d) = (d/p) = (1/p) =1. 


Hence the last equation (12) is satisfied by integers C and R=—2r. Thus 
(11) is equivalent to 2A or 2p. Hence the latter also represent 24. 


Lemma 2. If p=3 or 7% and A is any positive quadratic residue of p, 
2A is represented by the unique reduced form 2p or 20 of Hessian 2p? with 
all coefficients even. 


Here D=p*. Take A= 8pl, where / is prime to 2p. By (12), B= pb, 
b=81A—p. Then bd is a prime for infinitely many choices of 7. Since 
b=— p= 1 (mod 4), 


(1/b) = (b/l) = (—p/l) = (1/p), (p/b) = (b/p) = (214/p), 
(—2/b) = (2/b) = (2/p), (A/p) =1, 
(—A/b) = (—2/b) (p/b) (1/b) =1. 


To prove (7), note that, by Lemma 1, A represents every positive integer 
prime to 3. Since A= (x— y)?+ 2 (mod 3), A is divisible by 3 if and 
only if and are divisible by 3. Write y—=2-+ 3w, z=3Z. Then 

= 3L, L=2* + 3rw + 3w* + 32°. Since L=0 or 1 (mod 3), A repre- 
sents no number of the form 3(3n-+ 2). In Z replace x by X —Y and w 
by Y. We get » written in capital letters. By Lemma 2, » represents every 
positive quadratic residue 3n +1 of 3. Thus A represents every 9n + 3 and 
hence its product by 9%. Finally, if L=0, then e=0 (mod 8), and g, y, z 
are all multiples of 3, whence 2 is the product of a like form by 9. This 
proves (7). 
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Next, (8) follows from (7), since we proved that, when 3m is represented 
by A, m is represented by L and hence by p. The converse is true since * 


+ ay+y’?) X—c—y, dy. 


To prove (9), note that, by Lemma 1, p represents every positive integer 
prime to %. Since p= (x— 3z)? + y” (mod 7), pis divisible by 7 if and only 
if y=7Y. Then p=7R, R= 2224+ + 7X? 7Y?. 
Hence p represents none of the integers excluded in (9). For z= Z—2X, 
R becomes r= 2Z?— ZX + X?+ VY’, which is of type o in (10). By 
Lemma 2, o and hence r represents all positive quadratic residues 7n + 1, 2, 4 
of 7. Hence p represents their products by 7:49%. If R==0, then z=az 
= y= 0 (mod 7), and p is the product of a like form by 49. 

Finally, (10) follows from (9), since we proved that, when 7m is repre- 
sented by p, m is represented by # and hence by o. The converse is true 


since + 


(a? + 2y? — ay) =X? 4 2Y?—XY, X=x—4y, Y=2%e—y. 


The following results are not required in this paper. Lemma 1 holds 
also if p=—5 or 13. In the proof take d= 2pt-+yv, where v is odd and a 
quadratic non-residue of p. Lemma 2 holds also if p=: 5 and A is any non- 
residue of 5. As before, we prove that 


(13) xv? y? + 22? -+ zy +a2-+ yz represents all  25*(25n + 5); 
(14) 2x? + 2y? + 22? yz represents all 25*(5n +1). 


These forms are the only reduced ones whose doubles have the Hessian 10 


and 50 respectively. 


4. Forms of Hessian 6. The reduced forms are 


w = 2 + 2y? + 327, 2d. 


+ 62%, 


*For if w is an imaginary cube root of unity, the product of the factor «— wy 
of #? + wy + y? by the corresponding factor 1—w of 3 is X—wY. Also, (8) follows 
directly from Lemmas 1 and 2. If u is divisible by 3, e—=y-+ 3X; write y= Y—¥X. 
Then » becomes the triple of a form 4. 

7 For, if «? +2—e=0, the product of the factor w—ey by the corresponding 
factor 1—2e of 7=1-+2 -2?—2 is X—eY. Also, (10) follows from Lemmas 1 
and 2. If o and hence (#+ 3y)? is divisible by 7, y=27+7Z and o=7S, 
S =a? + + + 2%. For x= X—3Z, z=Y, S becomes p written in capital 
letters. 
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Suppose that »=10 (mod 16). Then »=2 (mod 4). Thus 
«= 2’ -+ 32? is even. Whether x and z are both even or both odd, «=0 
(mod 4). Hence y is odd and 10 =e + 2,e«=8 (mod 16). Ifzandz 
are odd, e==4 (mod 8). Hence r= 2X, z=—2Z, X?+3Z=2 mod 4), 
contrary to the above. Thus w represents * no integer of the form 2(8n + 5). 
This is used in proving 


(15) 7=2°+ y?-+ 62? represents all 9*(9n-+ 3). 


If represents 2m, then y= 2X, .—y=—2Y, m=X?+ Y? + 382, 
Conversely, the latter implies 2m = (X + Y)*-+ (X—Y)?+ 62%. Hence 
by (5) the even integers represented by 7 are the doubles of the integers repre- 
sented by g and include all ~9*(9m-+ 3), where m=2n+1. We next 
prove that + represents every positive odd integer a not divisible by 3. We 
employ 

Kk == an? + by? + cz? + 2ryz + 2zz, 


whose Hessian is H = aA—b, where A—bc—r*. Take H=6, A—16k, 
k=1 (mod 6). Then = 16ka—6 = 28, where B = 8ka—=3 is a prime 
for infinitely many choices of &. Thus 6 is not represented by w. Since b 
and the odd a are represented by x, « is equivalent to neither w nor 2A and 
hence is equivalent to 7. Since 


(—A/B) = (—k/B) = (B/k) = (—3/k) = +1, 


A = bc —7r’ is satisfied by integers c and r. 
If is divisible by 3, 3X, y= 3Y, r= 3T, T = 3X? + 3Y? + 22’, 
whence +r represents no 3(3n + 1). 


It remains only to prove that + represents a= 3A, where A is a positive 
odd integer 3n—1. We employ 


F = az? + 2By? + cz? + 4Ryz + 2szz, 
s=1+4A, B=8kA—3, k=l1 (mod6), 


with & chosen so that B is a prime. The Hessian of F is 6 if 


B+2AB+ +4 3—3ABe + 6AR?=0. 


* It was proved (Bulletin, 1. c.) that a positive integer is represented by w if and 
only if it is not of the form 4*(16n + 10). 
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Replace the first term B by its value and cancel A; we get 


Sb 6h*—BP—0, P=—3e—A—2. 


Hence 8k(1+ P) =0, P=—1 (mod 3). Hence ¢ is an integer. The 


congruence 48% + (62)*=0 (mod B) holds if and only if 


(—48k/B) = (—3k/B) = +1. 


The latter follows from 


(k/B) = (B/k) = (—8/k) = (k/3), 
(—3/B) = (B/3) = (8kA/3) = (k/3). 


5. Forms of Hessian 12. We shall prove that 
(16) § = 2a? + 2y? + 32? represents all ~ 8n-+ 1 and ~9*(9n-+ 6). 


If 8 is even, it is the double of a form (15). Henceforth, let z be odd. 
Evidently 5541 (mod 8). Let a be a positive integer such that a=3 
(mod 4), a4 9*(9n-+ 6). Then a is represented by g in (5). If z is even, 
g=0, 1 or 2 (mod 4), contrary to a= 3. Hence z is odd, 2° +y?= 
(mod 4), ~=2X, y=2Y. Hence a is represented by 2(X + Y)?+ 
2(X — Y)*-++ 32? and therefore by 8. 

It remains only to prove that 8 represents every positive integer 6 such 
that b==5 (mod 8), bA9*(9n +6). Then is represented by g in (5). 
First, let z be odd. Then x2? + y*?==2 (mod 4), whence z and y are odd. 
Thus y—2X, r—y=—2Y, and g=—2X*+ 32’, so that D is 
represented by 8. Second, let z—2f Thus r-+-y is odd. Since we may ] 
permute x and y, we may take ct = 2é, y odd. Then 5=g=4° +1+ 122 


(mod 8), and ¢ is odd. Thus 
(17) X= 4}(e—32), Z=i(e+2) 
are odd integers and 
(18) X? 4 3Z? = 2? + 32%. 


Hence we have a representation of b by g in which the unknowns are now all LE 
odd, as in the first case. F 
We shall need the fact that 


(19) e= 2’ + 3y? + 42? represents all ~4n + 2 and ~ 9*(9n + 6). 
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If e« is even, x= y + 2X, «= 4), where d is of the form (7), and con- 
versely. Thus (19) holds for even integers. Next, let a be a positive odd 
integer 9*(9n Then a is represented by gin (5). Thenz+y-+z 
is odd. First, let two of x, y, z be even. Permuting 2 and y if necessary, we 
may take x even. Then g —4X?-+ y? + 32’, whence a is represented by «. 
Next, let x, y, z be all odd. Changing the sign of x or z if necessary, we may 
take a==1, z==—1 (mod 4). Then X and Y in (17) are even integers for 
which (18) holds. Hence the second case reduces to the first. 


6. Corollaries. We require the results 


20 a = 32? + + 122? represents all ~4n-+1,4n + 2,9*(12n +11), 
Pp 


9*(12n + 8); 
(21) B=2? + 3y? + 122? represents all ~4n+ 2, 9*(38n+ 2). 


If m is represented by 8, 3m is represented by 32? -+ Y*-+ 4Z? with 
Y= 3y, Z = 3z, and conversely, whence (19) implies (21). 

If a is even, r= 2X, a is the quadruple of a form (6). Next, let a 
and hence x be odd. Then a=3 (mod 4). Any integer represented by a is 
represented by @ with z even. Conversely, if B==3 (mod 4), y is odd and z 
is even, whence B is of the form a. This proves (20). Finally, 


(22) y = 2? + 12y? + 122? 
represents all 4 4n + 2, 4n + 3, 9*(12n + 5), 9*(12n + 8) ; 


(23) = 3x7 + 4y? + 42? 
represents all 4n + 1, 4n + 2, 9*(36n + 15), 9*(36n + 24). 


If v is even, v is the quadruple of a form of type (5). Next, let v be odd. 
Then v=3 (mod 4). Write Y=y+z, Z=y—vz. Then v becomes 
§ = 32? + 2Y? + 2Z?. - Conversely, if 5==3 (mod 4), z is odd and Y+Z7 
is even=2y. Write Y—Z—2z. Then & becomes »v. Applying (16), we 
get (23). 

If y represents m, then 3m is represented by v with y and z multiples of 3. 
Conversely, the latter follows if v= 3m, whence ym. Thus (23) implies 
(22). 


Part 2. QUATERNARY QUADRATIC ForMs. 


%. Let + be a positive ternary quadratic form such that f=wu?++7r 
represents all positive integers. Let m be the minimum positive integer repre- 
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sented by r. If m= 3, f would not represent 2. Postponing the case m = 2, 
let m1. Let H be the Hessian of the form 2r all of whose coefficients are 
even. Since the minimum of 27 is 2, a known theorem* states that 2r is 
equivalent to a form 27 such that 47 — X? is a reduced form C = Ly? + 
2Myz-+ Nz* of Hessian 2H, where X = 2x-+ ty-+ sz, while the values of 
¢ and s are both chosen from 0,1. Further, we may assume + that M = 0 and 
omit the case L—=4, M—2,s=—0. Since the coefficients of X? + C=4T 
are all multiples of 4, we have 


(24) L+#?=N-+s?=0 (mod 4), M-+ts=0 (mod 2). 


The Hessian 2H — LN — M? of C is therefore a multiple of 4, whence H is 
even, H —2h. By hypothesis, f—u?+ 7 represents all positive integers. 
Hence 

F=4f=—4v7+ X¥°+0C 


represents all positive multiples of 4. Since C is reduced, its minimum > 0 
is ZL. If L>12, F=12 therefore requires that C—0. Since 12 is not a 
sum of two squares, this proves that L=12. Since t—0 or 1, (24) gives 
LI =0 or —1 (mod 4). Hence the only possible values of LZ are 4, 8, 12, 
3, 7,11. The least positive multiple of 4 not represented by 4u* + X? + Ly? 
is 1 = 28, 56, 24, 24, 84, 88 for the corresponding value 4, ---, 11 of L. 
Since F' represents 7, we have C </ for certain integers y, z, 240. Then 


(25) LC = (Ly + Mz)? + 4hz2? 
implies that 
(26) 4h = Li. 


Hence there is a limited number of values of f and hence of reduced forms C. 
The number of forms 7 is therefore finite. 


8. The simplest case is L—4. Since C is reduced, 2M = L, M=—0O, 
1,2. By (24), ¢=0, M=0O (mod 2), and N + s? is divisible by 4. 


First, let s=0. If M = 2, our ternary form is known to be equivalent 
to one with s=1. MHence there remains here the case M—0. By (26), 


* Landau, Handbuch... Verteilung der Primzahlen, 1909, p. 545. 
{ Eisenstein, Journal fiir Mathematik, 41, 1851, 141-90. 
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hX28. Since is a multiple 4m of 4,1 mS7%, and f=w+a?+ 
y? + mz® represents every positive integer (Ramanujan). 


Second, let s=1. We may write N—4m—1. There are two cases. 
(i) M=0. Thenh=JN, h=28, 
(271) 4f = 4u? + 4y?+ N2*, X (mod 2). 


This is of the form F = U? + X?+ Y?-+ N2*. Conversely, if F is a mul- 
tiple of 4, either z is even and U, X, Y are all even, or z is odd and just one 
of them is odd. Permuting them, we may take X odd, whence U and Y are 
even. In both cases, U = 2u, Y = 2y, X =z (mod 2), and we have (27). 
To show that F represents all multiples of 4, we apply (1). Let 
a= 4"*(8n +7) be a multiple of 4, whence 21. For z= 2*1,a— N2?= 
4'-1P, where P = 32n-+ 28—N=1 (mod 4). Thus P is positive and a 


sum of three squares. 
(ii) M=2. Thnh=—=N—1, ADH28, 


(272) 4f = U? + X? + W?+ (4m—2)2*, U =2u, 
A=—2 +2, W=—2y+2. 


Conversely, if this form is a multiple of 4, either z, U, X, W are all even, 
or z is odd, and just one of U, X, W iseven. Permuting them, we may take 
U even, X =z, W=z (mod 2), and deduce (27,). To show that it repre- 
sents all multiples of 4, take z= 2*1. Then 


a— (4m—2)22? = 4°19, Q =—32n + 28—4m + 2=2(mod4), 
Hence Q is a sum of three squares. 

THEOREM 1. For and 

f. =f: + yz both represent all integers. Expressed otherwise, X? + Y? + 


U? + (4m—1)2* and X? + W? + U*? + (4m — 2)2* both represent all posi- 
tive multiples of 4. 


9. The caseL=—8. By (24),¢—0, M is even, and N—4k—-s. Write 
M=2y. By 2MEL, »S2. We postpone to § 10 the case »—1. Here 
Then 4f—F, 


= 
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(28) F=¢+02, o=—4k—s— 2’, 

where X = 2a -+ sz, W—=2y+ vz. Thus 

(29) X=sz, W=vz (mod 2). 


Conversely, if / represents a multiple of 4, (29) hold. Evidently 
X =oz==sz (mod 2). If z is even, X and W are even. If z is odd, 
s?+ 2W? + o=0, —2v?=0 (mod 4), whence W=v= vz (mod 2). 

We have 8N—M*?=—4h= 8-56 by (26), whence o = N — 56. 
Also, N2=L. Hence 2=k=14. By 3, ¢$ represents exclusively all posi- 
tive integers except 4°(16n-++ 14). If the latter is a multiple of 4, it is 
a= 4?(64n + 56). For z= 2?, F represents a if ¢ represents 


g=—64n+ 56—c=0. 


If s=1, o and q are odd, and 4?q is represented by ¢ If s—v—0O, (28) 
gives f = u?+ 2? + 2y? + kz’, which represents all positive integers (Rama- 
nujan). Finally, let s=0, v1. Then g=14 (mod 16) only when 
; ==— 1 (mod 4), whence k = 3, 7, 11; o = 10, 26, 42, respectively. Lxcept 
in these cases 4?q is represented by 4. For z = 2?*1, 


B= a—o2? = 4°*1p, = l6n+ 14—c. 


If ¢=10, b=4(4n+1), and B is represented by ¢ If o— 26, 
b =16(n—1) + 4 is the quadruple of an odd integer; if n > 0, B is repre- 
sented by ¢. For c= 26,n—0,a— 4:56. If p> 0, take z= 271; then 
a— oz? is the product of 4-1 by 198 =6 (mod 16). Hence unless a is 56, 
it is represented by F, while 56 is not. Finally, let o—42. Then 
b=16(n—2)+4>0ifn>1. For n=0 orl and p> 0, take z= 27-1; 
then a—oz* is the product of 471 by 182 or 438, each =6 (mod 16). 
Hence unless a is 56 or 64 + 56, it is represented by F, while they are not. 


THEOREM 2. For2=k=14, s=Oorl, v—Oorl, 
a? + + 2y? + 2vyz + 


represents all positive integers except 14 when k=7 and 14 and 30 when 
k=11, with s =0, v—1 for all exceptions. Expressed otherwise, the form 
(28) with only square terms represents all positive multiples of 4 except 56 
and 120 in the cases mentioned. 
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10. The case L=8, M=2. Now 8f=G, 

(30) G=y+the2, y=2X?+W?+ h=2N—1, N=4k—s, 

where X sz, W—=4y+2. Hence 

(31) X =sz (mod 2), W=z (mod 4). 


Conversely, the latter follow if G represents a multiple of 8. Evidently 
W =z (mod 2). If z is odd, 2X¥?+1-+h=0 (mod 4), whence X?-+ N 
is even and X =s=sz (mod 2). If W—z—4r-+2, W+-z is divisible by 
4; after changing the sign of z if necessary, W—z is divisible by 4. Next, 
let z be even, z= 2¢. Then W and X are even, W—20. Then 


G = 4(o? + =0 (mod 8), o=€ (mod 2). 
Thus (31) hold. 


By (4), any multiple of 8 not represented by y is of the form 
a= 4¢-16(8n+ 7%). Take z=2% Since A112 by (26), a—hz? is the 
product of 4¢ by a positive integer = 2s + 1==1 or 3 (mod 8), and hence is 


represented by y. 
THEOREM 3. For 2=%=14,s—0 orl, 


f=wt+ + 2y? + yz + ke? 


represents all positive integers. In other words, (30) represents all positive 
multiples of 8. 


11. The case L=3. Then M—O or 1. By (24), ¢=1, M=s, 
N=4k—s. Then 12f—F, 


(32) F=¢+4h2?, 3X°+W’*, h=3k—s, 


where U = 2u, X =2a+y+s2, W=3y+ sz. Conversely, if F represents 
a multiple of 12, we may write 


U=0, X=W (mod 2), W=sz (mod 3). 


By (6), any multiple of 12 not represented by ¢ is of the form 
a==9°-9(12r-+8). By (26), hS18. Take z—3°% Then a—4hz? is 
the product of 9° by g=9-°12r+%72—4hk=20. Thus g=s (mod 3). If 
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s=1, G=—9-g is represented by ¢. Next, let s=0. If & is not divisible 
by 3, g is divisible by 3, but not by 9, whence G@ is represented by ¢. For 
k= 3, g —9(12r-+ 4) and G is represented by ¢ For k—6, g =9-12r 
and G is repreesnted by ¢ if r is not divisible by 3. When r is divisible by 8, 
take z = 2° 3°, whence a— 4hz* = 9*1-12(r—2) is represented by ¢. 


THEOREM 4. Forl1=k=6, s—0 orl, 
wu? + a2? + xy + + y? + syz + 


represents all positive integers. In other words, the form (32) represents all 
positive multiples of 12. 


12. L odd. By (24), ¢=1, N=4k—s, and M-+s is even. Write 
U=2u, X—2%+y+s, W—Ly+ Mz. 
Using (25), we get 
(38) K =4Lf = LU? + LX? + W? + 4hz’, 
(34) U=0, W=ZX (mod 2), W=Mz (mod L). 


Conversely, let K be divisible by 42. Since 4h—LZN—M’, 
W*? — M’z?=0 (mod LZ). Changing the sign of z if necessary, and noting 
that L = 3, 7, or 11 is a prime, we get (34,). Since L=—1 (mod 4), 
K =0 requires that either W, U, X are all even or that W is odd and just 
one of U and X is even. Permuting the latter if necessary, we may take U 
even, X odd. In both cases the first two congruences (34) hold. Hence f 
represents all positive integers if and only if K represents all positive mul- 
tiples of 4Z. Write 

W=X¥+2 D=4—1. 


Then we may remove the factor 4 from (33) and get 
(35) F=Lf=—¢o+h2, XE+ 
Conversely, if F' is divisible by L, 
(X + 2€)? = (mod L), X+4+2—Mz+ Lly—W, 


after changing the sign of z if necessary. Hence f represents all positive - 
integers if and only if F represents all positive multiples of L. For L=3, 


or 
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¢ is of type (8) and we may give a new proof of Theorem 4 by use of (35). 
We have 
(36) h=Lk—r, MM’), 0OSMSHBL. 


13. The case L=‘%. Since ¢ is of type (10), any multiple of 7 not 
represented by ¢ is of the form a= 49° - 49(%n + v), where v= 3, 5, or 6 is 
a non-residue of 7. By (26),h 7-21. By (36),h = where r= 0, 
1,2,or 4. Thus2=k=21. Takez—7*. Then 


a—he=G=499, 


Hence G is represented by ¢ if r 40, or if r —0 and k is prime to 7. Hence- 
forth, let r—=0, e=1, 2, or 3. Then G=49¢1y, where y= 
"n-+v—k. Next, takez—2-7* Then 


a—hz? = D = 49°18, = Tn + v — Ak. 


If v= 6, y=%n + 2 for 2, 8 = 2 for —1, 
= 7(n—1)+1 for «x =3. 


If v=5, y= +4 for 2 for «= 3, 
§= 7%7(n—1) + 4 for = 2. 


If v= 3, 8= + 2 or + 1 for «x —1 or 2. 


In all these cases, G or D is represented by ¢ if it is positive. It is negative 
only when y=6, x = 3, or v=5, If e >0 in the former 
case, z= 7*1 gives a— hz? = 49°(49-6—3), which is represented by ¢. 
But if e—0, a= 49-6 is not represented by ¢, since a— hz? = 49(6 — 32?) 
is positive only when z? = 0 or 1 and then the second factor is a non-residue 
of 7. In the second of the exceptional cases, a—=5-49°1, h—=2-49. We 
may take z—e?™, « prime to 7. According as m Ze or m<e, 


a— ha? = 49¢+1(5 — 49m-e),  49m*1(5 49e-m 


and the second factors are non-residues of 7 if positive. 
Finally, let yx = 3. Then G is represented by ¢ unless n is a posi- 
tive multiple of 7. In that case, take z—3-7*. Then 


a—hz?=49%1!m, m=—%n + 3—27 =7(n—4) +4>0, 


and the former is represented by ¢. 
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TurorEM 5. For 2=k=21;s=0,1; M—0, 1, 2, 3; M-+s even, 


ks?’ 
represents all positive integers except 42 when k = 21, s = M = 0, and 35 - 49¢ 
when k =14, s=M=0. In other words, for h in (36) and L=%, the 
form (33) represents all positive multiples of 28 except 28 - 42 and 28 - 35 - 49¢ 
in the cases mentioned. 


14. The case L=11. Then ¢ in (35) represents all positive integers 
except 121*(11n + v), where v ranges over the non-residues 2, 6, 7, 8, 10 of 11. 
Let it be a multiple a = 121¢-121 (lin-+ ) of 11. By (36), h =11k 
where 7 is 0 or one of the quadratic residues 1, 3, 4, 5, 9 of 11. By (26), 
h=2-11% Thus 3=k=22. Takez=—11% Then 


a—he? = G—=121¢9, g=121(11n + v) —11k+7=0. 


Thus G is represented by ¢ if r 0, or if 7 —0 and k is prime to 11. Hence- 
forth, let r= 0, =11x, or 2. Then G—121%1ly, y=11n + 
Next, take z= 2-11°. Then 


a— hz? = D=121%18, §=11n + v—A4Ak. 


First, let =1. For or 6, y=11n+4orlin+5. Forv=% 
or 8, 8=11n+ 3 or lln+4. Hence G or D is represented by ¢. 

Next, let-x«—2. For or 7, y=11n+1 orlin-+5. Forv—=2, 
y =.11n and G is represented by ¢ if n = 0, while if n > 0, 8 = 11(n—1)+ 
5 > 0 and D is represented. Finally, let v—8. If n—0, D=0. 

If n> 0, take z—3-11% Then ¢ represents 


a — ha? = 121¢1[11(n—1) +1]. 
THEOREM 6. For 22; s=—0,1; M=0, 1, 2, 3,4,5;M+s 


even, 


+ a? + zy + saz + 8y? ++ +s) yz + 


represents all positive integers. In other words, for h in (36), L—11, the 
form (33) represents all multiples of 44. 


15. The case L=12. By (24), ¢=0, N=4k—s, M=2, 
Also, h=3N — p’. 


(37) 12f = 12u? + 3X? + W? + hz’, X —2r-+ sz, W = by + pz. 


— 
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However this form is divisible by 4 also when z is even, and X and W are odd. 
Since the usual converse argument fails here, we apply devices appropriate 
to the necessary sub-cases. 

First, let » be even, w= 2v, whence vO or 1. In (37), W=2Y, 
Y = 3y + vz, whence 


(38) F=12f=a+hz2*, a—12u?+ 3X?+ 4Y?, 
X =sz (mod 2), Y=vz (mod 3). 


Conversely, the latter follow if F is divisible by 12. 

We have h ~12k—v7, where +r = 3s + 4y? takes the values 0, 3, 4, 7. 
By (26) and NZ2L,3=kS6. By (20), a represents all multiples of 12 
except a= 9°-9(12n-+ 8). For z= 3%, 


D=a—he=9d, d=9-12n+72—12k +7=0. 


If r=, d=7 (mod 12), and D is represented by a. 

Let r= 3. Then d is divisible by 9 only if k=1 (mod 3), whence 
k=4. If k4, d=3 (mod 12), and D is represented by a. If k—4, 
h = 45, d= 9(12n + 3), and D is represented unless 12n + 3 is divisible by 
9, whence n= 3m +2. In the latter case, take z= Thena represents 


a— = 9%1-12(3m +1). 


If >= 4, d=4 (mod 12), and D is represented by a. 

Let 70. Ifk is prime to 3, d is divisible by 12, but not by 9, when a 
represents D. If k—=3, d—=9(12n-+ 4), and D is represented. Finally, 
let k=6. Then D = 9*1-12n is represented unless n= 3m,m>0. Take 
z—=2:3°. Then a represents 


a— hz? = 9%1-12(38m —2) > 0. 
THEOREM %. For k = 3, 4, 5,6; s=0, 1; v=0, 1. 
u? + + sxz + 38y? + 2vyz + ka? 


represents all positive integers. In other words, (38) represents all positive 
multiples of 12. 


16. By (87), 


(389)  12f = y= 12u? + 1227 + W?, h= 12h — 22, 
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where W=6y-+ pz. Conversely, if F is a multiple of 12, we may write 
W=pz (mod 6). Since N= 44212, 24, k=83, 4, 5, or 6. By 
(22), y represents all multiples of 12 except a—9°-9(12n-+ 8). For 


== 3°, 


D=a—he? =9ed, d=9(12n + 8) —12k +p? =0. 


If »—1, then d=1 (mod 12), and y represents D. 
Let p= 3. If &k=—4 or 5, then d=9 (mod 12) and d is not divisible 
by 9, whence y represents D. If k —6, then d=9(12n-+1) and y repre- 
sents D. If k=3, e>0, take z= then y represents a— hz? — 9¢q, 
where g = 9(12n-+8)—3=9 (mod 12). There remains only the case 
p=k=3, e=0. Then a=—9(12n-+ 8) is not represented by F. For, 
a— hz? = 98, § = 12n + 8 — According as z is odd or even, 8=5 or 8 
(mod 12), whence 96 is not represented by y. 


THrorEM 8. For k=3, 4, 5,6; »w=—1, 3, u? + 2? + 3y? + pyz + ke? 
represents all positive integers except those of the form 9n + 6 when k=p=3. 
In other words, (39) represents all positive multiples of 12 except 12(9n + 6) 
when k = p=3. 


1%. L=12, p—s—1. Then k=4, 
h=12k—4=3-24, k=6. Hence 


(40) f=w +a? + + + yz + ke? = 4, 5, 6)... 


We first determine the integers represented by f when z is even, z = 2Z. 
Write Y=—3y+Z. Then 3f—F, 


(41) F=¢+h2Z*?, ¢=3u?+3X?+ Y?, Y=Z (mod 8). 


Conversely, F =0 implies Y?— Z?=0 (mod 3). Changing the sign of Z 
if necessary, we have Y==Z. We shall prove that F represents all positive 
multiples of 3 except 18, and, when k = 5 or 6, except 45. By (6), ¢ repre- ‘ 
sents exclusively all positive integers except 9*(3n + 2). If the latter be a 
multiple of 3, it is a=9°-9(3n+2). If e>0, take Z = then 


a—hZ?=D=9e1d, d=9?(3n+2)—h>0,. d=1 (mod 38). 


Hence ¢ represents D, and F represents a if e >0. Next, let e=0. Take 
Z4=1. Then a—hZ?=2%n+18—h=1 (mod 8) and is positive if 
n> 0O when k = 4, n >1 when k—5 or 6. 
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For the exceptional cases we employ (40) with y=0, z—1, whence 
f=w+2?+2+k. We have 


f=6 if : f=—=15 if 5/2 2 


THEOREM 9. The forms (40) represent all positwe itegers. 


This and Theorem 10 are equivalent to the fact that (37) subject to the 
condition X =z (mod 2) represents all multiples of 12 apart from the ex- 
ceptions noted. 


18. D=12, s=1, p=—3. Then h—12(k—1) 53-24, 
Thus 
(42) f=w? + 2? + az + 3y? + + kz? (k = 4, 5, 6, 7). 


We first determine the integers represented by f when z—2Z. Write 
X=2+Z, Y=y+2Z. Then f becomes 
(43) F=¢o+mZ*?, m=—4(k—1). 
By (5), ¢ represents exclusively all positive integers except a = 9°(9n + 6). 
For Z = 3*, D=a—mZ? d=9n+6—m. When k=4, 5, 6, 
d= 3, 8, 4 (mod 9) and D is represented by ¢ unless it be negative, i. e., 


unless n=O for k—=4, n=O or 1 for k=5, 6. Next, for ¢ >0, take 
Z = 31, Then 


G=a—mZ?=—919g, g=—9(9n+6)—m>0. 
When k = 5, 6, 7, g =2, 7, 3 (mod 9), and @ is represented by ¢. 
It remains to consider only the following cases: 


(1) k=4,a=—9°-6. Then* f=—a for u=—y=—0. 


(ii) 


=6or15. Neither a is represented by F, since m = 16. 
For y = 0, z =u 


Thenf—6 ifu—1,c—0; f—15 
if = 2, 


(iii) k—6,a—6 or 15. Neither a is represented by F, since m = 20. 
For =0, z= 1, f=uv? + or 15 if u—0 or 8. 


* But F does not represent any such a. For, if Z = 3me, e prime to 3, P= a— 1272 
is negative if m >e li m<e, P=I9mp, p = 9e-m- 6 — —6 (mod 9). 
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(iv) a=9n+6. For Z=—1, a—24Z?=—9(n—2) is re- 
presented by ¢ if n= 2 and n— 20 or 6 (mod 9) and hence if n+ 2 or 8 
(mod 9). Similarly, for 7 = 2, a— 24Z* = 9(n—10) is represented by ¢ 
if n=10 and n¥1 or 7 (mod 9). Hence a is represented by F' unless 
n =0, 1, or 8; whence a= 6, 15, or 78. We enquire whether 12a is repre- 
sented by (37), viz.. B= + 722%, with X —z even, where B =12u? + 3X? 
+ 

By (21), BA 72, whence 2? — 1, X —0, and X —z is odd. Thus (42) 
does not represent 6. Next, B—=12:15 for u—0, z=1, XY¥=—3, W=9, 
whence f==15 for u—0, r=y=—z=—1. Finally, B—12-78 for u=—83, 
2—=1, X =3, W= 27, whence f= 78 for u=3, r=z—1, y= 4. 


THEOREM 10. The forms (42) represent all positive integers with the 
exception of 6 when k =%. 

19. By way of summary, we note that Theorems 1-10 give 329 forms 
which represent all positive integers, and six forms which represent all but 
one or two small integers. We have not included the known forms u?+ 7 
which represent all positive integers, where 7 denotes 


+ y? + s2z2(sS7), x? + 2y? + sz? (28514), 
x? + 3y? + sz? (s= 3, 4, 5, 6). 
Including the latter 24, we obtain 
THEOREM 11. There are exactly 353 classes of positive ternary forms T 


whose minimum is 1 such that u? +-T is a quaternary form representing all 


positwe integers. 


Similar methods apply in the remaining case m = 2 (§ 7), and the num- 
ber of reduced forms 7 is again finite. The results will be published on 


another occasion. 
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Reduction Formulas for the Number of 
Representations of Integers in Certain 
Quadratic Forms. 

By E. T. BELL. 


1. Introduction. We shall illustrate by its application to three inter- 
esting formulas, summarized in § 3, an algebraic method whereby, knowing 
the number of representations of an arbitrary integer in a given quadratic 
form in r= 2 indeterminates, we can obtain a finite reduction formula for 
the number of representations in a quadratic form in rs indeterminates, s > 1. 
The formulas in any case can be obtained by strictly arithmetical methods de- 
pending on a simple use of the theory of binary quadratic forms of negative 
determinant —c; I prefer however to use an algebraic method, partly on 
account of its independent interest, but chiefly because such a method can be 
more easily extended. 

Let c, n, p> 0 be integers, of which p is prime; write (r, s) for the 
G.C.D. of r, s, and N [n =f] for the number of representations of n in the 
form f. The three formulas mentioned refer to 


+ cy’ + + ct’), 


so that N[n =f] is the total number of sets of integers x, y, 2,t=0 such 
that, for n fixed, 
n= 2? + + p(2? + ct*). 


Any integer t > 0 is of the form p*n, where (p, n) =1,a20. We shall 
call the expression of V[p'n = f],a > 0, as an algebraic function of N[n =f] 
the reduction of N[p°n =f], and similarly for any form at least one of whose 
coefficients is divisible by p. 

For c=1, 2, 3 the complete reductions of N[p"n =f] were stated by 
Liouville,* who however suppressed all details concerning the only cases in 
each instance which are not immediately obvious. His interesting and useful 
formulas therefore stand in need of proof, the more so as a direct use of them 


* Journal des Mathématiques, Vol. 10 (1865), pp. 48, 49; Vol. 11, p. 211. In §3 
here we considerably shorten Liouville’s statements of his results; also, in §§ 6-8 we 
derive formulas which unify several of his. 
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greatly abridges Pepin’s memoir * on Liouville’s numerous theorems concern- 
ing special quadratic forms which also were published without proof. 

A particular interest attaches to those reduction formulas which express 
N[p =f], a> 0, in terms of N[n =f] and algebraic functions of the divi- 
sors of the given integers n, p, c. Liouville’s theorems for c= 1, 2, 3 are of 
this type; equally simple results can be found for c—5, 7 and n properly 
restricted, also in any case where there is but one class of binary quadratic 
forms of determinant —-c, or where each order contains but one class. This 
will appear from the algebra for c—1, 2, 3, although the details are given 
only for these values of c. 


2. A difference equation. Let (m, 2) =1, n= 28m, s=0, (n, p)=1, 
and write o(j) =the sum of all the divisors of Then 


a(2%p'm) =a(p*)o(n), a=0. 


Now if N[p% =f] —y(a)¢(n), where ¢ is independent of p, a, and for p 
fixed (a) depends only on a, we need attend only to y(a) in the determina- 
tion of N[pn=—f]. The difference equation for y(a) in all cases with 
which we shall be concerned is of the type 


(1) ¥(a-+1) + = ko(p*), 


where & is independent of p, a. The solution of this equation is 


and hence, since o(p") = (p"1—1)/(p—1), we have 


a k p{r—(— a 


or, as it may be written, 


(3) y(a) — WO) — 


For m, n, s as stated, k is of the form x(s)o(m), where x(s) is a function 
of s alone. 


3. Summaries for c—1, 2, 3. Write ! 


(4) + p(2? + #), 
(5) = + 2y’? + p(2? + 277), 
(6) h =a? + 3y? + p(2? + 


* Ibid., Vol. 6 (1890), p. 1. 
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for p prime, and let m, n, s be as in § 2. 
Then, for the form f, 


(4.1) N[2%m =f] =cN[m=f], 

where c = 2, 6 according asa=1 ora >1; 

(4.2) p=—1mod4, N[pn=—f] —=N[n=f], 

(43) p=imod4, N[pn—f] =y(a), 

where y(a) is as in (2), (3), with k —160(m) or 480 (m) according as 
s=0ors> 0. 


For the form g, with p>2 (the case p—2 is treated in §5), we 
have 


(5.1) p=5,% mod 8, N[pn=—g] —N[n=—g], 
(5.2)  p=1,3 mod 8, N[p'n—g] —y(a), 
where is as before, with k = 8o0(m), 160(m), 480(m) according as 
s=0, s=1,3>1. 
For the form h, 
(6.1) p=1 mod 6, N[pn—h|] —=N[n=h], 


(62)  p==—1mod6, (p,3) —(p,n) (3,n) =1, 
=f] = N[pn—f], 620; 

(6.3) p=—1 mod 6, (n, 3p) = (m, 2)=1, n=28m, 
N[ =h] =y(a), 


with k = 80(m), 8(28*1 —3)a(m) according as s=0, s > 0. 
The enumerations (4.1) — (6.3) are evidently exhaustive for the re- 
spective forms; the only cases offering any difficulty are 


(7) (4.3), (5.2), (6.3), 


to which we return in §§ 6-8. 


4. Connection of f, g, h with principal binary quadratic forms of deter- 
minant —c. In order that n may be represented in the form z? + cy? of 
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negative determinant —c, it is necessary that (—c|n) —1, where (as al- 
ways henceforth) (r|s) is the Jacobi-Legendre symbol. 
that n be not representable in 2? + cy? it is sufficient that (—c | n) 

Consider all the representations of pn, where p, a, n are as before, in 


Hence, in order 


= — 1. 


the form 

F=2? + cy? + p(2+ cl’). 
If (—c|p)=——1, so that p is not representable in 2x*-+ cy’, then 
N[p'=F]~0, a>0, only when mod p. _ Put therefore 


pé, y= py; then 


N[ pn = 2? + + + | 
= N[p™1n = 2 + ct? + p(& + cn’) ], 


whence, repeating the reduction a sufficient number of times we have 


N[ptn = F] 


Hither this or a similar argument disposes of all cases in §3 except (7), 


and the like can obviously be extended to any form f =f, + p(f2+°':+fs), 


where 
fj = + + 1, S), 


and, with obvious modifications, to any forms f; all of the same kind in differ- 


ent sets of indeterminates. 


5. The Algebraic Method. To illustrate the simpler essentials of the 
algebraic method of finding theorems of the type § 3(7), we shall apply it to 
the form 


=a? + y? + 22? + 22?, 
verifying first that 


N[2%1m = $] = N[2%m = 2? + y? + 2 + 2), 


where (m, 2) =1 as before. Multiply within the second bracket by 2, in 
order to make the arguments 2%1m, 2%m agree, so that like powers of a para- 
meter q in the generating functions of both N’s may be compared: 


N[2%1m = = N[2%1m = 2a? + By? + 22? + 277]. 


| 
Vin — FF), 


al- 
rder 
—1. 
, in 


hen 
‘ore 


s)s 


Representations of Integers in Certain Quadratic Forms. 61 
The respective generators are 
N[2n + 2y? + 22 + 
where & refers to n ~ 0, 1, 2,° - -, and the modification of 


=X N[n— + + + 26°), 


from which are selected only those g” such that n==0 mod 2. This is accom- 
plished by changing q to —q and adding; 


93°(q) + = N[2n = 9]. 
Hence it must be shown that 
or, what is the same, 


93°(q) + 40°(q) = 265"(q"); 


which is correct, by the transformation of the second order. The verification 
is therefore complete. 

The particular point to be noted is that we select only those q” in the 
second generator for which n=0 mod p, where p (= 2 in this case) is the 
prime concerned in the given form ¢. This generalizes to, any p. 

The foregoing however is a mere verification, and we are interested in 
finding such reduction formulas ab initio. For this the order of the steps 
is reversed. We start from any known theorem concerning numbers of repre- 
sentations in binary quadratic forms, translate it into its theta equivalent, 
multiply the result throughout by theta constants to obtain a new identity 
which, when read as a relation between N’s, yields a theorem concerning 
numbers of representations. The appropriate theorems for binary forms will 
be found in Dirichlet’s Zahlentheorie, or in the papers of Pepin and others 
abstracted in Volume 3 of Dickson’s History. 


6. Form f=2?+y?+ p(2z? +77). The appropriate binary theorem is 


N[n—= 2? + y*] =4€(n), 


where (nm) = the excess of the number of 4¢ + 1 divisors of n over the num- 
ber of 4¢ + 3 divisors; é(r)é(s) = (rs) if (r, s) =1, and if p=1 mod 4 
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(p prime as always), é(p*) =a-+1. Since we need attend only to p=1 
mod 4 for f, we take only this case, and have, in the notation of § 2, 


N[pin = 2? + =a? + = 2 [pin 


as the simplest difference equation satisfied by € The translation of this 
into an identity between generators is (coefficients of q/, 7 = pn, are com- 


pared), 


ge" N [pn = + y?] + N[n = 2? + y?] ge" N[n=2* + 


in which & refers to n= 0, 1, -. We now replace this by its 6-equiva- 
lent ; the first = is the only one requiring any attention. To restate this = as 
a function of 6’s we must select from 


only those terms for which the exponent of q is a multiple of p. Let r be a 
primitive (algebraic) pth root of unity. Then 
p-1 
2, = p Ni pn + 
and hence the required 6-equivalent is 


p-1 
2 + = 2p 03? (q?), 


which is the analytic expression of our linear difference equation for N. 
Hence any consequence of this 6-identity will yield a true proposition. We 
wish here to derive from it a theorem concerning enumerations in quadratic 
forms in 4 indeterminates, and therefore we multiply the identity throughout 
by 63;7(q) or by 637(q?), choosing the latter as it gives a simpler result. We 
have then 


4?) = ++ (4?) = 2p6s*(q?). 


From the way in which the & in this was obtained, only multiples of p 
will occur as exponents of q in the power series for the first product on the 


left. Equating coefficients of q we get 
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N[pn = 2? + + p(? + 
+ N[ pn = (pa)? + (py)? + plz? + #)] 
= 2N [pn = px’ + py® + + pt?), 
+ +22 +2], = 16[2 +(—1)"Jo(m), 
the last by Jacobi’s theorem on 4 squares, where as always n = 2°m,(m,2)=1, 


s2o0. 


Hence, in the notation of § 2 we have 
y(a+1) + =16[2 + (—1)"]o (pm), 


whose solution y(a) is given by (3) with k = 16[2 +(—1)"]o(m), and this 
is equivalent to (4.3). 

In the foregoing analysis the use of pth roots of unity is not a necessity ; 
we might have used throughout the full generator in its N form. But in 
more complicated situations of a similar kind it is advantageous to have all 
generators expressed in finite theta form on account of the transformations 
thus suggested. 


Form g=2? + 2y?+ + 2t7). We start from the known N 
for the principal binary quadratic form of determinant — 2, 


N [2m = 2? + 2y*] 23 (—2 | d) =20(m), 


where (m, 2) —1, and & refers to all divisors d of m. From the properties 
of Jacobi’s symbol it follows that if r, s are relatively prime integers, then 
o(7s) =w(r)o(s). To obtain the initial difference equation for w we need 
the explicit form of w(p®), B > 0 (p prime, as always). We have 


(—2| = (—2 |p) 


and therefore 


according as (—2|p) =—1 or —1. Explicitly then 


o(p?) =3[1 + (—2| p)] (8 +1) 
+ 4[1—(—2| p)] +(—1)9], 
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so that 
40 (p ) = 28 +3 +(—1)6+(—2 | p) [28 +1—(—1)¥]. 
Henee, for a > 0, 
4w (ptt) = 2a + 5—(—1)*+(—2 | p) [2a +3 +(—1)*], 


4u(p") = 2a-+3 +(—1)*+(—2 | p) [2a +1—(—1)4], 
2a -+ 1—(—1)"+(—2 | p) [2a—1+(— 1); 


whence we have the required difference equation 
+ o(pt) = 20(p*) + (—1)* [(—2 | p) —1]. 


Referring to §3 (7) we see that it is necessary for our present purpose 
to attend only to this equation in the case (— 2| p) —1, so that 


(—2|p)=1, + = 20(p"), a>0. 
Proceeding now exactly as in § 6, we write down from this 


> N[np = 2? + 2y?] + N[n = 2? + 
N[n x? + 2y"], 


where & refers to n = 0, 1, 2,- - -, and (as in all that follows) (—2| p)—1. 
Translating as before to theta functions, with r—a primitive pth root of: 


unity, we get 
p-1 
> 93 (19q) 93 (724g?) + = (9?) 
j=0 
and therefore 
p-1 
6; 6; (q??) 03 
j= 
+ p63 (q?);(q??) 0s = 97"), 


which steps up from binary to quaternary forms upon equating coefficients of 
like powers of g. With (n, p) =1 we get 


N [prin =] + N[ptn=g] — 2N [pon =a? + 2y? + 2? + 24]. 


The last N being a known function we can proceed. For (2, m)—1, 


se 
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N [28m = 2? + y? + 2(2? + #?) |] =bdo(m), 


where b = 4, 8, 24 according as s = 0, 1, > 1, which follows at once from § 5 
and Jacobi’s theorem on 4 squares. Hence for n = 2*m, s= 0, we have 


b= —(—1)"} +4(2 1+ (— 
and therefore in the notation of § 2 (1) 


y(a+1) + =ko(p*), 
with here 
k = 2bo(m), 


whose solution is given by § 2 (3) with this value of &. This is equivalent 
to (5.2). 


8. Form + 3y*+ p(z? + 3/7). An even integer n is repre- 
sented in x? + 3y? only if m=0 mod 4. From the known theory of binary 


forms we have 
N[m = 2? + 3y?] = (—3 | d) =28(m), (m, 2) 
N [228m = + 3y?] = 68(m), s>0, 


the referring to all divisors d of m. When m==0 mod 3 the Jacobi symbol 
(—3| 3d) —0. Proceeding as in §7 we find 


48(p*) = 28 + 3+ (—1)* + (—3| p) [28+1—(—1)* J, 
which gives the difference equation 
8(p%*t) + 8(p?) = 28(p*) + (—1)*[(—3| p)—1], 
for which, by § 2 (7), we need discuss only the case p==—1 mod 6; whence 
8 (pet) + = 28(p*). 


As before, with (n, p) =1, we find from the identity between products 
of 4 theta functions derived from this, 


N[p™*1n = h] + N[ =h] =2N [pn = 2? + 38y? + 2? + 32?]. 


The last N is well known; it may be obtained if desired analytically from 
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the transformation of the third order of the theta functions. It can be 


written 


N[n = 2? + 8y? + 2? + 342] = 4[o(n) —{1—(—1)"} o(m)], 


where n = 28m, s=0, (m, 2) =1. Hence in this case 


w(a+1) =ko(p*) 


where 


= 4[1 —(— 1)* + (1 + (—1)*} (2%? —38)] o(m), 
and for this & we get (6.3) from (3). 
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Generalization of Certain Theorems of Bohl, 


[Second Paper]. 
By F. H. Murray. 


In a preceding communication with the same title * it was shown that 
a number of existence theorems of Bohl 7 can be generalized in such a manner 
as to apply to the solutions of a system of differential equations in the 


canonical form 


dx, OF dy oF 
dt OY: dt 02; 
near a point solution 7% =a; (‘—1,:~°-, m) when certain conditions are 


satisfied. In the present paper these conditions are replaced by less stringent 
ones; the methods of proof of certain existence theorems are very similar to 
those employed in the first paper, and these proofs are given here in an abbre- 
viated form. In addition, the asymptotic properties of certain trajectories 
are discussed by an extension of the methods of Bohl. 

On account of the more complicated form of certain quadratic forms 
which occur here, it has been convenient to leave undetermined certain con- 
stants which are determined explicitly in the first paper; this procedure, 
together with the algebraic method of transforming the canonical equations, 
reduces to a small amount the results common to both papers. 


1. Transformation of the canonical equations. 


Assume the function F(2,, tn; °° Yn) to be continuous, 
together with its partial derivatives of the first and second orders, within a 
region 


(D) |e | SD, ly | =D, i=—1,---,n, 
oF oF 

and suppose = 0, (t=1,---,n) if %—yx—0. The canon- 
02x; OY: 


* Amer. Journal of Math., vol. XLVII (1925), pp. 25-44. 
+P. Bohl, “Ueber die Bewegungen eines mechanischen Systems, etce.”, Crelle’s 
Journal, Band 127 (1904), pp. 179-276. 
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ical equations 
dx; OF dyi oF 


become 


(1.02) } 


and the equations of variation 


dé; n PR 
dt 2 dy 02, ér >> Oy 104 r 

(1.03) t=1,:- 
dni n  @ n ge 
2 ~ Axidy, 


possess characteristic exponents which satisfy an equation with real coefficients 


(1.04) A(a) =0. 


By hypothesis the functions F;,’, F;’’ possess continuous derivatives of the first 
order which vanish with the variables xi, y; (1, k =1,°°*, 1). 

From the classical theory of differential systems of the form (1.03) it is 
known that if a is a root of equation (1.04), the solutions can be associated 
in certain groups of the form 


= He” = ny’, 
(1.05) ett + téi’), ert (ni/’ + tni’), 


Ni (y) ett (ns tni (y- — ni 


Suppose another group of solutions represented by 


Zi’ = Gi’ = e%'ty,’, 
(1.051) 


dt 2 + 2 OY 0Yr 
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Then from the relation which is satisfied by any pair of solutions of (1.03), 
(1.06) (EP —C, 
we obtain, if a + a’ 0, 

> (éi’yi’ — = 0, 


and in general 


n n 
— nia) t [x Py, — Ya; ) 
i=1 i= 


é=1 


whence 


n 
yp, 


If a+ a’ = 0, we obtain in the same manner 
n n 
> (EP — )/ > — nia; ) t 
i=1 
>> — + > Pai) | 


eC’, 


from which further relations can be derived; but these will not be used in 
what follows. 

Suppose that no exponent is zero, and that the function F is real for 
real values of its arguments. Then it is known that to each group of v solu- 
tions of (1.03) corresponding to the exponent a there is a group of vy solutions 
with the exponent — a. 

If a is one exponent, the conjugate imaginary a is another, and a funda- 
mental set of solutions of (1.03) can be so chosen that the constants which 
appear in any group corresponding to a are the conjugate imaginaries of the 
corresponding constants appearing in the group with the exponent a. 

It will be assumed that there is only one solution in each group corre- 
sponding to a pure imaginary exponent; the solutions corresponding to a pure 

ay 


imaginary exponent ar, ( = > 0), can be normalized * in such a manner 


V—1 


* Amer. Journal., l. c., pp. 28-30. 


| 
nm, 
n, | 
«CF 
i 
| 
i 
4 
| | 
if 
| 
j 
4 
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that if (&, 9%), corresponding to ay, and (2, yi) corresponding to —a, are 
conjugate imaginary sets, then 


while if (2’;, y’s) is any other solution corresponding to —a,, 


Suppose 
and define 
Q (&, &n’) = = &) 
(1.11) r=1 Onr 
rss + + ini Am rss + dy ris 
The functions of (1.05) satisfy the equations 
(1.12) dé,’ /dt = agi’, dé’ /dt = & +08”, +, d&™/dt = + af, 
consequently 
(1.2) aé;’ Qn, n’)> a. Qn (E, 9), 
V5 


and if a+a’=0, 
— > [2s Q, (EP, 9) + Qn, ] 


.(r) ys = 0, r=1, 
i=1 


(1.3) 


n 
a > (Ey, — 0, 


r=1 
r= 2: cp. 


ay the pure imaginary exponents such that a,/W—1 


Suppose aj, a2,° °°, 
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>0, M), while ay,:,° * +, an have their real parts positive; 
suppose also any = —ar, and let the fundamental set of solutions be so num- 
bered that if the solutions (& i”, ju“) (r—=m, 
form a group corresponding to ax, the solutions (Enzi, nei) (T= 
m-+1,°°*,m-+v—1) form a group corresponding to an... Under these 
conditions make the transformation 


k 
k 


(1.4) 


the first summation extending over all the groups with the exponents a,,°- -, 
an, the second over the groups with the exponents an.1, ° °°, aon. 

If the transformation (1.4) is made in the right-hand members of (1.02), 
these become 


If for the functions Q; , Q@, we substitute their equals from (1.2), and 


rearrange the linear terms, then for the derivatives = : 4 substitute their 


equals from (1.4), linear expressions in the variables us”, vs‘ and their 
derivatives with the same coefficients are obtained on both sides of the equa- 
tions, their determinant being the determinant of the fundamental set of 
solutions and hence different from zero. Solving, we obtain 


du,™ /dt = + 6,™, /dt = + 
(1.6) /dt + + dv, /dt = 0, + + 


dtty,® /dt = + axtly,® + /dt = + + on™, | 


the functions being linear combinations of the functions 
The function Q(xy), representing the terms of the second order in the 
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Taylor expansion of F, can easily be calculated from (1.1), (1.11), (1.3), 
(1.4) ; we obtain 

Q(2,y) = + yiQu (ry) ] 
(1.7) 


kl 


rs 
M 


kirs 
where 3’ denotes a summation over the groups for which the real part of the 
exponent is different from zero; the coefficient Q Er is dif- 


ferent from zero only if a: + ansx = 0. 
Since the function F' is defined only for real values of its arguments, in 


general, it will be necessary to give ws", vs" such values that xi, yi, will be 
real; this result will be obtained if the following conditions are satisfied : 


(a) when ay is a pure imaginary, a, = pr V —1, pr > 0, Z and %, are real, 
and + Wr, = Lp Wr; 


(b) when ag is real, as = ge > 0, = Use, Ve" = Vex, Where Use and Vex 
are real; and (c) when a, is neither real nor pure imaginary, 


ar = Or + ade = Cr — dr > 0, 
= Urk + UE” = Ure + Wer 
= — = Ve, — Wek 
where Urk, Usk, Vex are real variables. 
Of the exponents a;, - * *, an suppose there are M pure imaginaries, K 


real, and ZL pairs of conjugates, each exponent counted just once for each 
group (1.6). Then these equations take one of the following forms, a, and 
azsr being conjugates: 

(a) dx; /dt = — + Xi, = + Yi, 
(1.8) dur;/dt — drUr1 +. Goss dvr; /dt drUri + 


(b) dur: /dt = Ur,i-1 + + Osi; /dt = Ur, t-1 ri + V ris 


dur; /dt = rULar,1 + Uris 


/ 
— 
an 


Murray: Generalization of Certain Theorems of Bohl. 73 


(1.8) duri/dt = Ur,i-1 + Qruri — rULer,i + Uri; 
dvyi/dt = Ur,i-1 — ri + rVLar,i + Vets 


(c) dutsr,1/adt — + + U 
AV r4r,1/at = — — + Vise, 


durr,,i/at = + Q’rUri + + U tar, is 
t= 
From the equations (a), (kb), (c) can be obtained certain fundamental 


differential inequalities similar to those employed by Bohl; 
if for each system of type (b), 


Vr Vr 
m=1 m=1 
we obtain 
Vr 
m=: 
Vr 
(1.9) +2 = hrm U em; 
m=1 
dg,(vr) /dt Gr? Ur1" 


Vr Vr 
= = Urm-1)” +2 = 1 Vim. 
m=2 


m=1 
For each system of type (c), if 
Vr Vr 
m=1 m=1 
then 


Vr 


Vr 
+ tara Lary, +z m1 + 


m=2 
Vr 
(1.91) + 9 (ttrmU Utar,mU 
m=1 
d/dt[gr(vr) + Gr (Uzsr) | — — », 
Vr 
Vr ‘ 
m=2 m=1 VLsr,m V 


* P. Bohl, “ Sur certaines équations différentielles d’un type général utilisables en 
mécanique,” Bull. Soc. Math. de France, tome 38 (1910), p. 18. 


| 
k) 
| 
i” 
| 
| 
| 
| 
i 
| 
| 
| 
| 
i 
| 
| 
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From the manner in which equations (1.6) were derived, it is seen that 
if ax, a1 are conjugate imaginaries, the functions ¢;“, ¢;‘” are conjugates, 
as are Wi’, Yi“ ; and the functions X;, Yi, Urs, Vrs are real functions of the 
present variables. From (1.08), (1.7), we obtain 


M 
(1.92) F=Fot +(%i? + Yi?) + Fs, 
i=1 


kirs 
in which the constants A,,** are real. 

Suppose (a): |% | Sa, Sa, | ure | Sa, | vrs | S a, a region such 
that the original variables 2;, y; lie within (D) when the new variables lie 
within (a). The functions Xi, Yi, Urs, Vrs possess continuous partial deriv- 
atives of the first order with respect to the new variables, which vanish with 
their arguments. If 9(%i, Yi, Urs, Vrs) is the largest of the absolute values of 
these derivatives at the point Zi, Yi, Urs, Urs), then 


| Xs | S rHi, rUre, Tyre) [ (| +1 Ge |) +B (| Ure | + | 
and similar inequalities hold, for other values of +, for the other functions. 
Instead of a single trajectory consider two, represented by the variables 
Yin Urs, Vrs, Te + U4, Yi t Urs + Vrs + U're} 


since equations (1.8) are valid for both, we obtain a similar set of equations 


(1.81) dx;’/dt — piyi’ /dt piri’ + etc., 

such that if X’;, Y’i, U’rs, V’re represent the differences between the unprimed 
functions at the two points considered, the new equations (1.81) contain the 
primed variables and functions where the equations (1.8) contain the un- 
primed variables and functions. Equations similar to (1.9), (1.91) will also 
be obtained. We shall have 


|X’ | Sg trai’, + ryi’, Urs + ture’, Ure + 
[EC |) +3 (le! | + | 
0<7r< 1, ete. 
2. Certain fundamental inequalities. 
Suppose 1,, +, lu, mi, M2,° *, Mxsz any set of positive numbers, 


and 81, S2,° *, Su, ti, to,° * *, te, 
a set of real variables. Then from the inequality 


Vm m 


| 
| 

| 

| 

| 

| 

| 
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and the introduction of the constants 
(2.1) k = hk, k S gZl, N=n—M, 


T=(1/k)V2M+N, A=—gvV2M+N 
we find 


M K 
4=1 r=1 r=K+1 


=F tel) + 
(2.2) r=1 


+ me +E Algo lee) + 1Z 


M K vr 
+] te |) +3 S meget] | 


K+L vr 


+ Mr? (| | + 2L+r,s | ) 


r=K+1 8=1 
Suppose = let 
K+2L 


G(u) = gr (Ur) 


and 
Vr 
m=2 
(2.3) 
Vr 
m=2 
K+2L K+2L 


H(u) = hr(utr), H(v) = hr(vr): 
Since the quadratic forms h,(u,), hr(vr) are positive definite, there exist 
positive constants A,, Ar such that these forms satisfy the inequalities 
hr (ur) — (ur) = 0, he(Ur) (He) = 0, 


and it will be assumed that these constants have their greatest values.* 


* It follows that the quantities X,, ,. are roots of the characteristic equations of 
the pencils of forms 
h,(u,) —A,9,(u,), h,(v,) —Ag,(%,), 
respectively. Since the transformation 


carries the second form into the first, the characteristic equations have the same roots, 


the smallest of these roots being the values of 4,, Ky 


| 
| 
| 
| 
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Suppose P, Q constants such that Q = = P, (s=1,° °° 
K +L); then 
(2.4) | SVNQ VGE(z), Dare |S VN/P VGE(z). 


If 
K 
r= +91") + Mr? Jr (Ur?) 
K+L 
then from (1.9), (1.91), if q is the smallest of the constants Axst; 
Ags, and if 7? =0, t= to, we obtain easily 
dr jdt < [—q+29(0 + VN/P) (A+ VNQ)]G(w), t—to. 
If in (a), 
q 
(2.5) g< 
VN/P)(A+ VNQ) 


then dr’/dt = 0. 


Suppose this inequality for g satisfied within (a). If the two trajec- 
tories lie in (a) for S¢St,, and 7” for t=), while G(u’) > 0, then 
dr’/dt <0, t=t,. Suppose = 0 for —¢,; at some time ” —0, 
and 7” <0, t’. Hence there exists also a time < for which 
dr’ /dt =0; but at this time 7’ < 0, hence é@r’/dt <0. From this contra- 
diction we conclude that 7” <0 when ¢t—#,, and <0, 
Hence from the definition of 1’, G(u’) > 0, tt, and if G(u’) 
t=t,, then G(u’) > 0, 

Under the assumption 7” = 0, and the preceding inequalities, 


K K+L 
d/dt G(u’) = H(u’) +2) rs’ +E Urs’ + 


r=1 8 r=K+1 8 
| 
K K+L 


ai’ | + | yi’ | + | Urs” | | |] 


> H(u’) —2g VNQ VG(w’) VN/P) VG(w’) 
= [q—29 VN/P) VNQ] G(w’). 
If the constant g satisfies (2.5), there can be found a constant o such that 


q>o> 2 VNQ (T+ VN/P); 


76 
T,8 
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and 


dG(u’)/dt > (q—o)G(w’), 
[G(w’) Jest, > E(u’) Jeato. 


From the preceding results we have the theorem: If o < q, and |, ls, 
Uy, mi, * are any positive numbers, then if a neighborhood 
(a) of the origin ts so chosen that within it the partial derivatives of Xi, Yi, 
Urs, Vre are less in absolute value than either of the quantities 


whence’ 


q o 
VN/P)(A+VNQ) 2VNQ (C+ VN/P) 
where A, T, P, Q, are the constants previously defined, this neighborhood has 


the following properties: 1. if the coordinates of two trajectories lie in (a) 
for StSt, and tf when t=—t, G(u’) €0, 


M K 
K+L 


Gr (0 rar) | G(u’) < 0, 


then 1°<.0 fort, < and 
[G(u’) Jeet, (t-te) ] rato 


2. If the coordinates defining two trajectroies lie in (a) for tZ2t=t, 
(t, >t), and tf for t= G(v’) ~0, and 


M K 
4= r=1 
K+L 


+.B + ge(w'sor)] — G(0’) S0, 


then s’ <0 for tp >t =t,, and 
[G(v’) > (tort) TG rato 


The last part of the theorem can be demonstrated in the same manner. 


3. An application of the energy integral. 


M 
In what follows it will be assumed that the form > + (%;?+ 9,7) 
i=1 


is definite ; suppose the constant A so large that for all values of the variables 
concerned, 


— A[G(u) + G(v)] S D = A[G(u) + G(v) ]. 


kirs 


v7 
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Assume p, v, p positive constants which satisfy the inequality 


put + p/P) <1. 


Since the first and second partial derivatives of F; vanish with its arguments, 
we can assume that in (E), a sub-region of (a), 


(3.1) | [> (zi? Yi”) (Urs” + Ure’) |. 
T;8 
Then suppose 6 > 0 such that the region defined by 


(S) G(u) Sw, G(v) Sv 
lies within (E), and let (R’”) be the region 


M 
(R”) + Ge) S 


From the energy integral, if a trajectory lies in (R’), (S) at ¢—¢,, in (R”), 
(S) at ¢=—?,, 


M M 
+ | t=t, = + [ +> + + Arg is + F; | t=te 
i=1 i=1 kirs 


[ D’ Ars ter Vis + t=t, 


kirs 
S pd? + 4Avd? + p[d? + pd? + 1/P (28° + 2v8") 

(3.2) + < ®. 

Let (H) be a region lying within (E), and such that in (H) the partial 
derivatives of X;, Yi, Urs, Vrs are less in absolute value than 6 where 

6< ____. 
2V (V2M +2VNr/P ) 

Suppose d so small that the region 


M 
> + 9?) S@, Sra’, G(v) 

lies within (H), and construct the regions 


(A) 492) <a, 


M 
(a) > (Zi? + S pd’, 


G(u) Svd2, G(v) Svd?. 


| 
i=1 | 
| 
(B) 
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It results that tf (Zi, Yi, Urs, Urs) lies in (a, B) for t= t., m (A, B) for 
t= t,, and if G(u) =rd*? for t= t,, then dG(u)/dt > 0, 

For when ¢ = ¢,, from (3.2) 
dG(u)/dt = H(u) —20V NQG(u) 


[Val VE + VIN/P VG(u) + 


> [qv —20V NOQv (V2M +2VNv/P)]@ > 0. 


4. Haistence theorems. 


With the definitions of the domains (A), (a), (B) of the preceding 
paragraph, we obtain at once, by the reasoning of the first paper ($4, pp. 
40-42) the following theorems: 1. To every set of values (Zi, Yi), (Vrs) of 
the region 


M 
+97) G(v) Sva’, 
i=1 


can be adjoined a set of values (urs) of the region G(u) S vd? such that the 
trajectory having these initial values at t =t, remains in the region (A, B) 
for t> to. 


2. To every set of values (Xi, Yi), (Urs) of the region 
M 
= (Zi? + Yi") = pa’, G(u) = 
i=1 


can be adjoined a set of values (vrs) of the region G(v) S vd? such that the 
trajectory having these initial values at t = t, lies in the region (A, B) for 


t< to. 
Let (H’) be the region (H) for the particular set of values 1, 1, =: > 


= ly = 1, my = = mM <1; from § 2 the quantity 7” formed 


for two trajectories which lie in (H’) for t= ¢,, 
9 
(xi? + yi”) + m?G(v’) — G(u’) 
i= 


must be positive or zero, otherwise G(w’) would increase indefinitely with ¢. 


Consequently 
M 
(4.1) G(u’) (ai? + yi?) + mG (v’), to; 


and tf two trajectories lie in (H’) for t=t,, and tf at any instant t’ > to, 
vi’ = yi’ = Vrs’ = 0, then urs’ = 0, and the two trajectories coincide. 


ts, 
| 
| 
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In the same manner we obtain the result that if two trajectories lie in 
(H’) for tSto, and tf at any instant U < ty, vi’ = yi! = Urs’ = 0, then 
Vrs’ = 0, and the two trajectories coincide. 

From the last two theorems it follows that the set of values (urs) cor- 
responding to the given set of values (Zi, Yi, Vrs) in the first existence theorem 
above is uniquely determined; these values, functions of the remaining vari- 
ables, will be denoted by [w,s]: similarly, the values (vrs) whose existence is 
asserted in the second existence theorem, will be denoted by [v,s]. The 
continuity of these functions follows at once from (4.1) and the corresponding 


inequality 


M 
(4.2) (ai? + yi?) + m’G(wv’), tS. 
i=1 


If a set of values (Zi, Yi) satisfying (a), § 3, is fixed, the functions 
Urs — [ure], — [vre], = 1 a Vrs r=] K 21 
are continuous in Uys, Vrs Within and on the boundary of 
If no set of values of u,;s, Vrs exists for which these functions vanish simul- 
taneously, from a theorem of Bohl a set of values must exist on the boundary 
for which one of the systems of equations must hold: 


Urs — [Urs] =—ANUrs, Vrs — [Urs] Ure, N > O, 


(1 +X’) ttre = [tre], (1+ A’) [tre], 


Consequently, 
= G@([u]). or = 
and therefore either 
= G([u]) or (1 = G([v]). 


But this is impossible, since the values [urs], [vrs] satisfy the inequalities 
(B). Consequently values must exist for which 


Urs = [ure], Urs = [Vrs ] 


from which we obtain the theorem: to every set of values of (a) can be ad- 
joined a set of values of (B) such that the trajectory possessing these initial 
values at any instant remains in the region (A, B) for all values of t. 


| 
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If two trajectories remain in (A, B) for all values of ¢, we shall have 


M 
» (vi? + yi”) + m’°G(v’) 2 G(w’), 


M 
(ai? + yi”) + 2 G(r’), 
i=1 


from which 


Hence if at any instant x;’ = yi’ = 0, then Urs’ = vrs’ = 0, and the two tra- 
jectories coincide. 
5. An mequality. 


Let 1, m, h be positive, and let a, B satisfy the inequalities a = 1, mq,4, 
gr*4 S B. Assume k& a number satisfying 


(a ha 


2B(M +N) 


If 
M 

= et [22S (ay? + yi?) + m2G(u’) —G(v’)], 
i=1 


where 2;’, Yi’ Urs’, Urs’ are the differences for two trajectories in the region 
(C) where the partial derivatives of X;’, Yi’, Urs’, Vrs’ are less in absolute 
value than k&, then from (1.81), 
M 
dw/dt = 2hw + + 227? Xi + yi’Vi’) 


i=1 


r,s 

— A (v’) + 2D 
Now 
[Xi], | |, | Tre’ 


? | 


(| | |) + | + | 1) 


Sk/a V2(M+N) V DP (xi? + yi?) + mG (u’) + G(r’) 
Also, 


+m 


Urs’ | + S gr84| |] 


SB FN) VES FEW) 


| 


82 Murray: Generalization of Certain Theorems of Bohl. 


and consequently 


do /dt = 2hw + + —2 (2m + en) 
M 

{SP(ai? + yi?) + m’*G(u’) + G(r’) }] 
i=1 


+2N)]o 
! a 


+ g — (M+ + qm?G(u’) 
a 
2k SBk 
> + (MN), — (M+). 


From which 
d/dt [e**tw|] = 0. 

If w» is positive for any value of ¢ > to, w will therefore increase until 
one of the trajectories leaves the region (C); if o 0 for ¢=—f, but the two 
trajectories do not coincide, G(v’) > 0, hence dw/dt > 0 and » takes on posi- 
tive values. Hence if both trajectories remain in (C) for ¢ > fo, and w is 
bounded, then » < 0 unless the two trajectories coincide. 


6. Asymptotic character of certain trajectories. 

Suppose h a positive constant less than gq, g,, and suppose A, B so small 
that the region 
(J) |, | re | SB 
lies within (a) and in (J) the partial derivatives of X;, Yi, Urs, Vrs are less 
in absolute value than y, a constant to be fixed later. Then we have the 
theorem: The constants y, H, H can be so chosen that if H < A, H <B, 
and (K) is the region 
(K) |,| 9: | SA—HA, | | ors | SB—AH, 


then tf a trajectory o lies in (K) fort = ty, suppose (vrs°) its set of v-values 
at tt: if an arbitrary set of v-values is chosen such that 


(Urs = Ure’)? (Gr 


then there is one and only one trajectory whose v-values coincide with these 
for t= to, which remains in (J) for t > to, and approaches o asymptotically 
as t—> , in such a manner that the differences of the coordinates, multiplied 
by e*t, remain bounded for t > to. 
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Let the coordinates of o« be represented by the variables Zi, Yi, Urs, Tre, 
and suppose e so small that the region 


|,| |S A +2(e*—1)H, | ure], | vre |S B+ 2(e*—1)H, 


lies in (a). Choose a constant @ < so small that for 4; —6@StS tho, 
—1)H, | are |; 
Then if 

(6.01) | | | Ore | Seed, 


| Gre | SB—H + (e*—1)H. 


we shall have 


|, | |S A+ 2(e%—1)H, 
Tre + | S B+ 2(e%—1)H. 


| tire + 


Making the substitution 


(6.1) w= en, 

Ure = tre + Ure = Tre + Wee 
we obtain from (1.8) the equations for the new variables &, 7i, Zre, Wre- 
These become 


7? 


(a) dé&/dt=hé: — pini + Xi’, dyi/dt = pis thy t+ Yi’, 


(b) dzr,/dt = (h+ dwr;/dt = (h Gr) Wri + 


dwyi/dt = + (h— Gr) wri + r=1:-:-K. 


(6.2)  dan/dt = qr)én On’,  dwr/dt = — qr) + Ver’, 
dz,,i/dt = 2r,i-1 + (h + Gr) Zri — 
dwyi/dt = + (R— Gr) Wri + + Fed, 
(c) deters /dt = qr’trr + Gr) + 
= — + — Gr) + 
= + t+ (h + Gr) + 
dwysr,i/dt = — Yr’ Wri + (h— Gr) + 
2's 
where 
= X Yi» Ure, Ure) — Xi (Zi, Yi, Urs, Tre) |; 
Ure’ = Yi, Urer Ure) Ji, dirs, Tre) |, ete. 


the variables 7;. 9;. Urs. Urs being given by the substitution (6.1). 


7=t—t, > —8, 
| Gam 2+ 
J 
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The right-hand members of (6.2) are bounded and satisfy a Lipschitz con- 
dition for t= to — 0, if the dependent variables are in the region specified 
above, hence the existence of solutions in a certain neighborhood of ¢ =, is 


assured. 
Suppose the forms Gn(z), Gn(w) defined by 


Ga(z) = (Gr + Gi (w) = (dr —h) 
18 
If ¢ = to, and the inequalities 
(6.3) = (é? +77) +Gi(z) SH’, Gi(w)S H?, 
i 


are satisfied, then if H,(z) is obtained from H(z) by replacing g, by gr +h, 
while Hi(w) is obtained from H(w) by replacing gr by gr —h, we obtain 


as in § 2, 


M 


= (0 + + 
i=1 


M 


If Pr, On, Pr On are constants such that 
Ai(z) Hi(w) =v Gi(w), 


Qn = (Gr + = Pr, = (qr 
then 


5 LE tnt) + Gala) ] + 


— 2% [V2M > + V Gi(z)] X 


| + 42) ta = +4 uw) G,(w) | 


If 2p is the smaller of 2h, ji and if H< H/VP,, H < H/VPy, H < H while 


(6? + + Gi(z) = H?, Ga(w) H? 


( 
if 
( ( 
(( 
be 
su 
in 
ex] 
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then 
(6.4) d/dt[ ~ (é:2 + 9:2) + Ga(z)] = — V2M + VN/Pi + VN/Pi] 


(V2M + 


We shall assume 7 so small that the expression on the right multiplying H? 
is positive; then the left-hand member is also positive. 
Also, 


Ga(w) S —H(w) + (qr — | 


= (w) xX VGi(w) 


Now suppose Gi (w) = H?; then if + 97) + Gi(z) S H?, 


(6.4) S t+ 29 VON (V2M + VN/Pa + VN/Pi) 


it will be assumed that 7 is so small that the term on the right multiplying 
H? is negative. From the choice of H it follows that when the inequalities 
(6.3) are satisfied, the first inequalities for the present variables are satisfied 
also. It will be assumed that A, B are so small that (J) belongs to the region 
(C) of §5. 

On‘ the basis of these inequalities, we can show * that to every set of 
w-values satisfying Ga(w) S H? at least one set of values (é:, nt, 2rs)o Can 
be chosen satisfying 


(6.5) (&? + 4”) + Ge(z) S 


a 


such that (wrs)o, (£i, nis 2rs)o are the values for ¢ = 7? of a solution remaining 
in the region (6.5) for t > to. 

Now suppose a set of v-values (V;s) chosen, such that Ga(V —v°) S H?, 
and let a solution of equations (6.2) be determined in accordance with the 
existence theorem above, such that at wre —=—Urs° + Vrs, (S=1 


* Bohl, 1. ¢., p. 230. The proof of the proposition stated above is so similar that 


it is unnecessary to repeat it here. 
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‘vr, P= K+2L), while the solution remains in the region desig- 
nated for ¢ > to; for ¢ >¢, the variables &, yi, 2rs, Wrs satisfy 


|, S A, 
| dirs + | Urs + wWrs | = 


while for ¢ = ft, Ure° + Wrs = Vrs. The trajectory defined by 


Urs = Urs + Vrs = Ure + Wee, 


satisfies the first part of the theorem; it is evident that the differences 
multiplied by e”*, remain bounded. 

The second part of the theorem follows easily from § 5. In these inequal- 
ities suppose 1 m=—1, and suppose the quantity w» constructed for two 
solutions remaining in (J) for ¢= ¢, in such a manner that the coordinate 
differences, multiplied by e’*, remain bounded for ¢ > fo, the constant h being 
such as to satisfy the inequalities of the present paragraph. Then 


M 


is bounded, and consequently = 0 fort If the differences v;s’ vanish 
at ¢t = to, it follows that the remaining differences x;’, yi’, Urs’, vanish also, 
and the two trajectories coincide. Hence when the v-values are fixed, only 
one trajectory having the properties stated can exist. 
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Applications of the Determinant and Permanent 
Tensors to Determinants of General class 
and Allied Tensor Functions. 


By M. CRAMLET. 


1. Introduction. The geometric scheme of representing the elements 
of a detetrminant as an arrangement of rows and columns or of layers, sec- 
tions and sheets * becomes cumbersome in the general study of determinants 
of elements with several indices. The natural analytical method in these cases 
is to confine attention to the general term of the determinant. Various meth- 
ods are used to determine the signs of this term and to exclude the terms that 
should not appear. The method most like that used here is the scheme of 
Cayley + of using a square determinant in this capacity. Recently {§ papers 
have appeared dealing with square determinants by a new method. It is 
along these lines that the work of this paper will proceed. Tensor notation 
will be used throughout but that is merely for considerations of symmetry. 
The results are valid for any tensor character that might be assigned. It is 
important to note, however, that with the transformation character that we 
shall assign our elements, the laws of transformation of the functions of the 
elements can be read from the indices. The modifications for different trans- 
formations on the elements will be at once apparent. To illustrate, suppose 
Jrst 18 One Of n® elements that appears in our work, and that it is a coefficient 
in the trilinear form gp0,7°y%z7, then since the ground form is invariant the 
indices in grst are all covariant. The transformation character of the deter- 
minant or whatever function of these elements that is discussed in our work 
can be read off at once. On the other hand suppose the g;s: are the n* coefti- 
cients in a system of m quadratics n, then the first index 
on the g is not of tensor (or transformation) character. It might be well in 
this case to indicate it as grsr. The modifications for this situation can be 
made in any equations in which the first index is not multiplied by inner 
multiplication by a tensor. 


* Hedrick, Annals of Mathematics, Vol. 1 (1899), page 49. 

t Trans. of Cambridge Phil. Soc., London, Vol. 8, Part 1, 1844 (2-43), S. 11, p. 85-8. 
(Reference from Lecat’s Lecgons sur la Théorie des Déterminants & n dimensions.) 

t Murnaghan, Math. Monthly, XXXII, 1925, p. 233. 

{ The writer: Annals of Mathematics, Vol. 27, No. 4, September, 1926. This will 
be referred io hereafter as paper I. 
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2. Definitions and Conventions. We define a tensor* by the equation, 


On by bn 


Dp 


where the n! terms formed by permuting the a’s and multiplying each by an 
arbitrary scalar, are added. The n! scalors are distinguished by the m! per- 
mutations p. When the scalors are +1 according as the permutation p is 
even or odd we have 
Qn 4 
(2.2) 


This tensor vanishes unless the a’s and b’s are the same set of distinct integers 
when it equals + 1 according as the substitution to bring them to the same 
order is even or odd. This will be called the determinant tensor.+ When 


the scalars are all equal to 1 we have 


On (2.3) 


which will be called the permanent tensor. It vanishes unless the a’s and b’s 
are identical sets. For the non-vanishing cases the values are 1 if the sets 
are composed of distinct integers. If there are q groups with pi: °° py 
repeated or equal indices p,+---pq >m the value of the tensors is 


An inner product of y and a set of elements such as 


&n 


will be called a tensor function of the elements. The y’s are of the tensor 
character indicated yet they are absolute invariants. If a 6 is.used in place 
of the y we have a determinant function, and likewise, using a 7, a perma- 
nent function. For functions involving elements with more than two indices 
we may use combinations of these three. The degree of a tensor function will 
be equal to the number of elements occurring in a term. This term will 
replace the less significant word, order. The set of m indices in the same 
position on the elements of a term will be called a range. Greek letters, only, 
when repeated will indicate a summation. Repeated letters are called umbral. 
For the non-repeated we venture the name virtual. The virtual ranges on the 
elements will be called fixed as is customary. The class is one greater than 


* Properties of this tensor were given in a paper read before the Society April 3, 
1926. 

7 Murnaghan has suggested the name “ generalized Kronecker symbol.” Either 
(2.2) or (2.3) for n= 1 becomes the Kronecker symbol. 


— 


Tensors to Determinants of Generall Class and Allied Tensor Functions. 89 


the number of umbral ranges, i. e. one greater than the number of tensor 
multipliers when each multiplier has one virtual range. If there is but one 
fixed range the class is equal to the number of ranges. 

These definitions and certain general properties will be illustrated by the 
function, 


This is a class 3 (or 3 range) determinant function of degree n. This ex- 
pression has n*" components corresponding to the integers assigned the virtual 
ranges. The value of the component of such an expression when the virtual 
ranges take the sequence 1- - - n is ordinarily called a determinant or in case 
7's replace 8’s, a permanent with appropriate names for combinations. We 
shall not be careful, however, to make this distinction and will speak of 2.5 


for example, thinking of any component, as a determinant. 


3. Elementary properties of the fundamental tensors. The elementary 
properties of a tensor function are apparent in the notation. 


(1) The virtual range on the multiplying tensor z or 6 is permanent 
or alternant. 

(2) Any substitution on an umbral range of elements only is equivalent 
to the same substitution on the virtual indices of the multiplying tensors 
carrying the same umbral indices. 

(3) If the tensor function has a y for a multiplying tensor the fixed 
range is in general neither permanent nor alternant, otherwise the fixed range 
is permanent or alternant according as the number of &s having one range 
multiplied by inner multiplication (i. e. summed) is even or odd. The proof 
of the property is as follows: Consider the right member of (2.5) (to illus- 
trate) when the 7’s are rearranged by a substitution S. They may be ar- 
ranged in the original order by a permutation S-! on the g’s. But this re- 
arranges the remaining umbral ranges by a substitution S-!. The order of 
these, however, may be restored by renaming the umbral ranges. This will 
leave the umbral indices on the multiplying tensors only, permuted by sub- 
stitution S. If the number of 8’s is even these indices may be arranged in 
the original order and the determinant has been brought from one form to 
another with a different arrangement of indices in the fixed range without a 
change in value. If the number of §’s is odd there will be a change in sign 
if S is odd and the range is alternating. An illustration of this is found in 
the two range determinant 


I crn) (an) Irie, I (3.1) 
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The parentheses on the range letters indicate that the range is alternating, 
while the brackets will indicate that the range is permanent. The indices in 


a fixed range need not be distinct. 
(4) Alternating ranges may be interchanged. For example 


(tr) tral tn) (an) 
for the number of changes of sign to bring (sn) to (tn) is balanced by the 
number to bring (tn) to (Sn), (Sn) and (tn) being permutations of the same 
integers. 
4, Permanents, Transformations of Coordinates, Inner Products. The 
following combinatory properties of the permanent tensor will be found use- 


ful: * 


Pi pn 41 
or more generally: 


The proofs follow readily from (2.1) and (2.3). 

Analogous to the method of finding the equation of transformation of a 
determinant + the equation of transformation of a permanent may be found. 
For covariant elements ars the permanent is 


Any interchange of the virtual ranges of the expression leaves its value un- 
changed, for example, interchange (s,s.), follow this by interchanging (p:p2) 
and then note that 7s?a--- on — 77P2-+- Pn and the expression has been changed 

to original value without a change of value. By the definition of the tensor 
(4.3) may be written 


Consider a typical term of this expression under & in which the p’s have been 
permuted by a substitution S. Before allowing the p’s to be summed, per- 
mute the elements by S. Then summing with respect to the p’s this term is 
Ar,o,°** rag, Where the o’s are a permutation S of the s,°--s,. The entire 
expression is then 


gn 


* For the determinant tensor see reference (3), Article I. 
7 See Paper I. 
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In a similar manner we arrive at the fact that any permanent range may 
be fixed and the others umbral if there are an even number of 8’s. This may 
be illustrated by 


(on) (Tn) (un) (fin) 
(8n) (tn) (tn) (4 ) (4 11) 


Let = (8%) be substituted for the first tensor. Consider a typical term in 
which the permutation of the o’s is 8S. Permute the a’s by S. Then per- 
mute all umbral letters on the elements by S-!. If there is an even number 
of 8’s the total effect is to permute the fixed range to order S(rn’) = (pn) say. 
Summing all such terms we have 


(pn) (Tn) (un) ( f 
9 
(rn) (tn) (un) (mn) (4.12) 


as equivalent to (4.11). The method of changing from a fixed alternating 
range to a fixed permanent range will be illustrated in (5.1) and (5.2). 

The transformation of a permanent. We may represent the transformed 
permanent by 


’y,s, Which we write for brevity 
, = --/[pn] , qr’ len] 
[Tn] [8n] ~~ [rn] (a Prin) 1 Tira] [sn] (a’ Pron) 


Ox, dag 


0x 
= 1, Nn! [rn] 02’ p, [sn] ( anBn ) 
In this expression the [a,] appearing on the partial derivatives is a range of 


a permanent and therefore permanent in character, so we can continue the 
equation 


[on ] [an] 


If the summation in this equation is extended only to distinct sets the nu- 


merical factor is removed. A bar under the umbral letters will indicate this 
convention. Thus 


This equation is identical in form to equation (6.4) of paper I. There the 


J’s represent the Jacobians or determinants of the transformation, while here 
the K’s represent the permanents of the transformation. In (6.4) paper I 


e 
e 
) 
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the inner product of two alternating Jacobians reduces to the square of the 
Jacobian. This simplification of the inner product of permanents is not 
possible since the absolute values of non-vanishing permanent components 
differ. Unlike determinants the sets of indices in brackets need not be equal. 


Inner Products of Determinants and Permanents. There is a very close 
similarity in the inner products of these two functions. 


Let Dw (en) = ( densn) 
(pn) : 
(an 
and E yy) (dn) ~ (€anan) 
then the inner product of these will = (dens) 


=n!>d ) (densne, 
p 


Now consider a typical term where the permutation on the p’s is 8S. Permute 
the d’s by S. This will leave the s’s permuted by S. Sum the p’s getting 


to within a sign 
n! where S(sn) — (on). 


The sum for all permutations S under & will be 


(qa 
(dean) 


If we set 
fa? = (4.6) 


the inner product 


(8) =n! 8" (for) Prem (4.7) 


(On) an Can 


Since the (yn) ranges are alternating this gives the ordinary multiplication 


theorem for determinants 


= Fis) with (4.6). (4.8) 
For permanents, every step in the proof will hold down to (4.6) and (4.7) 


but the argument by which (4.8) was derived depended upon alternating 
indices and is not valid for permanents. We will here state the results, if P, 
Q and RF are permanents with elements 


Yor and where 757 = p*qap (4.9) 


then 


n 
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5. Alternative forms for determinants. Decomposition. There are 
methods of writing a determinant alternative with (2.5) for example 


or 
eee ore 5 
J tral (sn) (tn) = 5 it Jon8nTn (5.2) 


Equation (5.1) is obvious when it is noticed that [pn] is a permanent range, 
there being in this case an even number of inner multiplications with 6’s hence 
all of the n! term in the sum for [pn] are equal tO 9p..3¢9n)¢taye Hquation 
(5.2) may be proved directly or verified in this manner. By making a sub- 
stitution in the (s,) and arguing as we have above, (s,) may be seen to be 
alternating. By changing s to p in the right member of (5.2) and multi- 
plying by §e1--- on the member is multiplied by m! since the inner product 
Leese 
of two alternating ranges is mn! times the ordinary product. But this gives 
thie ale ividi ‘4 
the right member of (5.1) which equals »! 9;,,3¢¢n)ctn) + By dividing by n! 
equation 5.2 is obtained. In these equations the first range is permanent. 
Corresponding equations are obtained when the other ranges are fixed. 


5. Cayley originally considered higher determinants as sums of two 
range or ordinary determinants, and methods of decomposition are current in 
the literature.* A single illustration will serve to show how this is done with 
the present notation. In the notation of (2.5) the fixed range is virtual, so 
the number of ranges is one greater than the number of tensor multipliers 
of the elements. When the expression is written 


Sn 


the expression in the brackets is a class 2 determinant since but one 8 is used 
in expressing it. The first range in this case must be considered as being 
the bipartite range 110;, -** Ton. When o is summed there will be n! such 
determinants. Hither § may be taken out of the brackets. The extension to 
the types of determinants classified at the close of article 4 is apparent. 


6. Multiplication of Determinants. Demonstration of the applicability 
of the method to more general functions than heretofore contemplated. The 
muliplication of two cubic determinants of order 2 illustrates the process. 
However, for file multiplication the process will be carried out in its greatest 
generality. 


*L. H. Rice, American Journal of Math., 1918, p. 242. 
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Element Multiplication. For the product of the two cubic determinants 
g and h we have 


h = 8%%2 §7172 §21b2 


8182 tite 
8182 tite 110171 1202T2 
where KIBY Graz (6.2) 


The class of & is hinds by counting the 6’s and adding 1. In general then 
for classes p and gq the product & will be of class p—1-+q—1+1= 
p+q—1. In the present case r, and a are permanent indices (property 3, 
article 3) so can be taken identical. The point so apparent here is that in 
this method not a property of the 8s is utilized, we might have 8's, z’s, y’s 
or in fact it is immaterial what they are. The class of & can be increased 
to p+ q. There will by this method always be two fixed ranges, in this case 
the ranges r and a. The range 7 can be changed to p,°-+ pn and the deter- 
minant multiplied by 1/n! 7 Spe where p is summed. If r had been alter- 
nating we should multiply by 1/n!dm---o otherwise the product would 
vanish. Now the class will be equal to the number of ranges. Remarks such 
as this one on raising the class are so trivial by the analytic method that one 
hesitates to make them but from the geometric treatment they seem of more 
importance. 

File Multiplication. The analytical analogue of the geometric nomen- 
clature, “file multiplication” is “inner multiplication.” Beginning with 
equation (6.1) and setting ¢=c—p and summing (the number of choices 
of umbral indices is the product of the numbers of alternating ranges, each 
leading to a multiplication theorem) we get 


BiB2 
8182 Bi Be Y2 | h h 


8183 


8182 * ry, 1101 
where Kos = Gray heey (6.3) 
8:82 [x 1202 + 1101 


from the interchange of umbral indices B,; and £2 and use of alternating 
property of ot Thus 


=k (bybe) 

[rire] (8183) 
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where as before we exhibit the permanent and alternating ranges by brackets 


and parentheses. 


Multiplication of a p range and a q range tensor function. 
Tay’ * *Tqay are arbitrary sets of quantities which may be permanent 


tensors, determinant tensors y’s or otherwise. 
plied are 


— On « Tn U1 

A Pay Sy + Sn ty tn Uy 
and 


The determinants to be multi- 


where A contains at least one alternating range indicated by one 6 and B 


contains at least one 6 (but an even number of alternating ranges by property 


3, art. 3) while the remaining ranges are either permanent or alternating as 
indicated by the x’s and 8’s._ The total number of &s in B will be designated 
by K (1). The number of z’s is immaterial. 
contains the indices which must be a permutation of 1- - so the 


product of these 6’s in the product AB will be 


2 My Ln 
uv 


The last 6 in both A and B 


When this sum is multiplied with the product of the h’s the sum of h’s is 


uv 


In (1) the sign of a term is + 1 according as the class of the permutation of 


subscripts of v’s is even or odd—thus in (2) the sign of a term will be 


according as the permutation with respect to indices of any letter say a, is 


even or odd. Consider a term of this sum and let S represent the permuta- 


tion of the v’s in it. The v’s may be brought to former order by a permuta- 


tion S-! of the h’s which has the effect of permuting the other indices by S-}. 


All of the umbral indices may be brought to the original order but this will 


This is of no 


leave the lower indices on the z’s and 8’s permuted to order S. 


consequence on the z’s but to bring the indices of §’s to the original orders will 


require a permutation S-! on each of K —1 8s so if K —1 is even the sign 


of the term is unchanged, otherwise it is changed when S-! is an odd permu- 


tation. The a’s are in the order S-! to which they were brought by the per- 


mutation on the h’s. In case the sign is unchanged it is of class S which 
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Pp 


is also of class S-1 which indicates that if a,° * ay» is the new order of the a’s 


the sign of the term is 6%-:-% . By substituting 


the product may be written 


En Ni +++ Nn T1 TnOn.+-+ Tn 


The class of C is found by the number of 7”’s, x’s, &’s equal to p—1 
+ q—1, the original number, — 2 lost by inner multiplication + 1 regained 
by the artifice + 1 by which the class exceeds the umbral ranges = p + g —2. 


The Associated Contravariant Tensor. The following method of obtain- 
ing an associated contravariant tensor shows that the existence of an associated 
contravariant tensor depends upon the property of an alternating tensor. Let 
the three indices of g;s: be covariant. The only possible determinant in this 
case, other than a permanent, is given by (2.5). The equation of transforma- 
tion may be written 


Ox 
1? “Dy 
G [Tn] (8n) (tn) G (n) (tn) ) (6.5) 
vr, 
So OG tn) 2 ( OXp, 
09’ por 0x’, 09 dg 
OU’ Odapy 0%, Oxp 


The first partial of the determinant with respect to the element has the co- 
variant character of the determinant and an additional contravariant character. 


If the determinant had no permanent range, products such as (7) would 
not appear in the equations of transformation. There would then he an even 
number of alternating ranges, hence an even power of the Jacobian of the trans- 
formation would appear in (6.5) and in (6.6) so if the first partial is divided 
by the determinant the resulting tensor will be of contravariant character 
wholly. Only an element with an even number of indices, when all repre- 
sent tensor character, can have an associated tensor of opposite rank. To the 
writer’s knowledge it is only such elements that are considered as representa- 
tives of a line element in differential geometry. 

A “ Bibliographie des déterminants 4 plus de deux dimensions” to April 
1924 is given by M. Lecat, Bruxelles. 
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Transformations Leaving Invariant Certain 
Partial Differential Equations of Physics. 


By R. D. CARMICHAEL. 


We shall write the general wave equation in the form 
(1) 
Any solution of this equation will be called a wave-function. 


For the case when n is 4, the transformation %=y, 
t,= V—1ct reduces equation (1) to the classical form 


x? dy? 02? Ot? 


of the wave equation in theoretical physics. This is the justification for call- 
ing (1) the general wave equation. 

The principal object of this paper is to determine the generators of the 
group of point transformations (the coordinates of the points being complex 
numbers) under each transformation of which a wave-function is a relative 
invariant. The group contains a certain linear group as a subgroup (see § 6). 
It is shown (see § 11) that this subgroup and two particular Cremona trans- 
formations generate the whole group when n > 2. When the equation is taken 
in the form (1), one of these Cremona transformations is real and the other 
is imaginary; but they are both real in the case of the classical form of the 
equation in mathematical physics. 

It turns out that the group involved is identical with a certain group of 
conformal transformations (see §3) which has been the object of several 
investigations. It was necessary, however, to redetermine the generators of 
this conformal group, since certain imaginary transformations are quite as 
important for our problem as the real ones which have been the object of 
previous investigations. 

Incidental to the principal object, we have determined the corresponding 
groups for two other partial differential equations of physics (see §2). (An 
additional equation is treated in the following paper by Mr. Goff.) 

In each of the three problems, the case n = 1 is trivial. The case n = 2 
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(see §%) is exceptional in character and is shown to lead to a well-known 
group of transformations on two variables. For this reason its exceptional 
features are not treated in detail. 

Incidental to our main purpose we obtain new proofs of four important 
theorems stated below as theorems IV, V, VI, VII. The principal result is 
contained in theorems VIII and IX near the end of the paper. 


§1. The partial differential equations defining the transformations of 
the wave equation into itself. Let the function v(y1, y2,° + *, Yn) be a solu- 
tion of the wave equation 


n 
(2) — =), 

k=1 
and let us make the transformation 
(3) Yo» * Yn) wF0 
where w is a function of 2, 22,° + +, Z, to be chosen conveniently and where 
* Yn are taken as functionally independent functions of %2,° °°, 


Zn Which are to be so determined in connection with w that w shall satisfy 
equation (1) whenever v satisfies equation (2). Every wave-function v is 
thus to be transformed into a wave-function u by means of (3). 

Since 1, yi, yiys (14 J), Yi? — Yn’, are solutions of (2), it follows from 
(3) that w and the y; must be such that w, wyi, wyiy; (147), W(yi? — yn’) 
are wave-functions. These facts imply that w, y1, ye, °° *, Yn must satisfy 
the following system of partial differential equations: 


n 07 
f (); 
k=1 
(4) n OY: 
i Y j 
— =0, 
> (0y:/dax)? = t=—1, 2,°° +, n—1. 
k=1 k=1 


From relation (3) we have 


{ 
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Ow py dy; 
i=1 1 OY: dy; OY: ” 


From these equations it is easy to see that wu, as defined by (3), is a solution 
of (1) provided that equations (4) are satisfied. 
Let us use p for the common value of the n sums 


n 
(dy: 2, °° +, a. 


If we employ 8;; to stand for 1 or 0 according as j is or is not equal to i, 
then the conditions in the last two lines of (4) may be put in the form 


n 


dy; — 


(6) Ox, 02, ijP J 1, 


Since the Jacobian (functional determinant) 

OY; 
| 
( ) 027. 


is different from zero (owing to the assumed functional independence of the 
y's) is follows from (6) that p can not be identically zero. 

The results so far obtained may be put into the form of the following 
theorem : 


THEOREM I. Let w, 91, yo,* Yn be a set of functions of 2%, %2,° 

a, such that ys, Yo, * * Yn are functionally independent and w340. Then 

a necessary and sufficient condition that every solution of (2) shall be trans- 

formed through (3) into a solution of (1) is that the functions w, y1, Yo; 

‘5 Yn Shall satisfy the following conditions :—a function p exists, not iden- 
tically zero, such that 


(8) 


and the equations 
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Ow dy": 

(9) = ( 02,2 ) 
n 0?w 

10 = = 0, 

2 0%," 


are verified. 


The totality of transformations (3), subject to the conditions just named, 


obviously form a group. 


Corellary. A necessary and sufficient condition equivalent to that in 
theorem I is that each of the following functions shall be a wave-function: 
1, 2,° °°, 2,°°°, n—1, J—1+1,°°°, 2; 
w (Yi? — Yn?), 1, 


§ 2. Certain related groups of transformations. The method used in 
the preceding section may be employed in the proof of some results related 
to theorem I. 


THEOREM II. Let w, yi, Yn be a set of functions of 2%, %2,° 
Such that * Yn are functionally independent and Then 
a necessary and sufficient condition that every solution of the equation 


(11) +v=0 
shall be transformed through the relation 
into a solution of the equation 
(13) 
OK 
is that the functions w, Y2,° Yn Shall satisfy the following conditions: 


a function p exists, not identically zero, such that 


Oy; 
k=1 02%, Ox; 


(14) 


and the equations 


(15) ft 


2 w— 
k=1 02%. 02}, 027." 


> 2 ow ) = 0 


n-2 
*In §5, these conditions are replaced by the simpler one, w =p See also 


the latter part of § 5. 
7 See the preceding footnote. 
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n 
(16) > 


—> 
k=1 
are satisfied. 

The necessity of these conditions is easily proved by aid of the following 
solutions of (11): sin y; and cos y; (or ev'V-1 and e-u'V-1) taken in conjunc- 
tion, and sin (yiVs) sin (yjV¥1—s) where 1-47 and s is a constant differ- 
ent from 0 and 1. In order to show that the conditions are sufficient one 
may proceed as in the proof of the corresponding part of theorem I. 

The partial differential equations of the group of transformations 2; 
into y; and w into v such that uw, as defined by the transformation vu = 
V(Y1, * Yn), Shall be a solution of the equation 


(17) > = 0 

k=1 
whenever (41, Y2,° °°, Yn) is a solution of the equation 
(18) (dv/ayx)? = 0 

k=1 


may be determined in a similar manner. The following sufficient conditions 
are readily obtained: 


> (dw/dr,)? =0; 


k=1 
n a 0 j 
tur], 2,°°-, 23 
k=1, Ox, 
OY; 


i564); 
k=1 


k=1 


> (dy; = =z (dyn/0r.)? 40, 
k=1 


+=1,2,°°°, 


That these conditions are also necessary may be shown by means of the follow- 
ing solutions of (18): v=1, v=f(yit y;V¥V—1), v= f(y: 
(‘4 7), f(t) being an analytic function of 

On multiplying the second of the foregoing equations by @y;/dx1 and 
summing as to t from 1 to n, we find (by aid of equation (21) in the next 
section) that dw/dr, =0, 11,2, - -, , so that w is a constant. Without 
loss of generality we may take it equal to unity; and this we do. 

Then we have the following result: 
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THEOREM III. A necessary and sufficient condition that every solution 
of (18) shall be transformed by the relation u=v(y1, Yo, * Yn) into a 
solution of (17), the functions Yn Of V2,* * *, Un being func- 
tionally independent, is that a funtcion p shall exist, not identically zero, 


such that 


OY: Oy; 
§3. The transformation group of theorem III. The transformation 
group of theorem III, transforming 2, %m into Y1, Yo, * *, Yn Con- 


tains the corresponding transformation groups of theorems I and II. We 
shall therefore treat first the transformation group of theorem III. 
Let J denote the functional determinant 


OY; | 


| 


Then by the multiplication theorem for determinants, we have 


Oy; | 
0x; 0x; 


—— | OY; n 
k 


the last member having been obtained by use of (19). For J we have there- 
fore the value 
(20) J =m 


Since J £0, it follows that p can not be identically zero. 

If we multiply both members of equation (19) by dy;/dx; and sum as to 
j from 1 to n, we have 

Le ( = Ox} 0x1 ) 02x, 3 3 

On holding / fixed and allowing & and 17 to vary, we have a system of linear 
equations in the n quantities inclosed in parenthesis. The determinant of 
the coefficients of this system is J and is therefore different from zero. On 
solving the foregoing system for the indicated quantities, we have 


Thus equations (19) imply equations (21); it may be shown similarly that 
equations (21) imply equations (19). 
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If we take y:, y2, °° *, Yn to be the coordinates of a Euclidean space, 


then the element ds of arc length is given by the equation 
ds? = dy,? + ++ ++ dyn’. 


OYm 
But dym Axx, m = 4, 2, 
k=1 


Substituting these values of the dym into the preceding equation and simpli- 
fying by the aid of (19), we have 


(22) ds? = dy? + dyo* dyn? = p(da,? + dian’). 


Hence the mapping effected by our transformation is conformal in character. 
Conversely, if a transformation of the y’s into the 2’s is such that equation 
(22) is satisfied, then it may readily be shown that equations (19) are satis- 
fied. Hence the transformation group involved in theorem III is the group 
of conformal point transformations in a Euclidean space of n dimensions, the 
coordinates of the points being complex numbers. 

$4. Determination of p for the conformal transformations. If we con- 
sider the equations 

Ti = Li (Yr, Yo» * * * Yn) 


of the inverse transformation, then on differentiating with respect to aj we 
have 

0x; 

02; 

Multiplying by 0ym/0x;, summing as to j from 1 to n and reducing by aid of 
(19) we have 


04 m 02; 
(23) —" =p —, i, m—1,2,°°-, 
0x; m 


These equations are useful in considering the relation of a transformation to 
its inverse. 

From equations (19) or (21) we may obtain by elimination a system of 
equations involving a single one of the functions y; alone. For this purpose 
we have from (19) by differentiation the equations 


yj Oy; Op 


k=1 02}; k=1 Ox}; 027.021 0x) 
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On multiplying by dy;/éa,, summing as to 1 from 1 to n, and simplifying by 
aid of (21), we have 


4 _ Op _ Oy; 


p k=1 02% i=1 021 02,021 0x1 
The second term of this equation may be simplified by means of the relation 


l=1 021 9x 02%, 02% 


obtained by replacing & by 7 in (21) and differentiating both members of the 
resulting equation with respect to a;,; thus we have 


Op dys 
(24) 0x" 0x1 021 


From (24) may be derived some consequences which will be of use later. 
If we sum as to / from 1 to n we have 
Op 
For the case when n = 2 this shows that y; and y2 are both wave-functions. 
On multiplying both members of (24) by 0y;/0%m and summing as to j 


we have 


Op 
02m 021? Ste 021 


(26) 


j=l 
On summing both members of this equation as to / we have 


2—n Op 
(27) ; 2 02m 


We may now determine * explicitly the possible values of p. For this - : 
purpose it is convenient to write 


(28) p= 


* This is based on the work of R. Beez, Zeitschrift fiir Mathematik und Physik, 
20 (1875), 253-270; it is necessary to repeat the argument for the sake of some addi- 
tional formulae not given by Beez. Some simplifications in the computation are also 
effected. Moreover, in our problem, we have to examine a case not treated by Beez. 


n 0 OY; | 
1 
1 a 
| 
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If then we substitute into (19) and (21) the values of dy;/da, and Oy;/0rx 
as given in (23) we have 


Ox; Oa; i Ox; 


om 


Differentiate both members of the last equation with respect to y; and in the 
result replace 02;/dy;, and 02;/dy, by their values as given by .(23) and sim- 
plify; thus we have 


= 


or 


0x; 
(30) )+ 


From the evident relation 
we have by means of (23) the relation 


1); 


o da; dy: 


2 Oo dx; 2 


Combining this with (30) we have the relation 


0 ( An; 
31 
( ) 0x; ) 
Differentiating with respect to 2m we have 


0? ( ) “1 ( 0x4 


rae... 
5 dy, J 0x; OL 


; 

+ (ae Dan \ 
+ $3; 


For the first and third terms in the brace { } substitute their values as given 
by (31); on simplifying we have the following equation: 


| 
J Ox, 
o dx; Oo Oo 
is 
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6? 0x; 


On interchanging m and j and comparing the resulting equation with (32) 


itself we see that 


02; 02% dy 02m Sim mlz ) 


When n > 2 it is possible to choose 1, j, m all different; if we do so we have 


, n>2, 4, j, m all different. 


Multiplying by 07m/dy: and summing as to / we see that 


Po 
(34) a 0, 14], n> 2. 
From this it follows, when n > 2, that 0c/0x; is independent of each of the 


x’s except 7. Hence we have 
(35) o = X,(%1) + X2(t2) + Xn(An), n> 2, 


where X;(2;) is a function of 2; alone. 
In (33) take 7 =7 and m1; then we have 


02m 2 dc Oc 0 0a 
For each m there exists a y: such that 0%m/dy: is not identically zero. Taking 
m —=1 we see that 60 /dx;7 is independent of ¢ if 1 belongs to the set 2, 3, 

-, n. Then taking m= 2 we see that is independent of if 
belongs to the set 1, 3, 4,---,m. Hence when n > 2 it is independent of i 
fori—1, 2,:-+,n. But from (35) it follows that 0’o/dz;? is a function of 
a; alone. It must therefore be a constant. Calling this constant 2¢ we may 


therefore write 


(37) n> 2. 
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Before undertaking to evaluate o completely it is desirable to obtain still 
another simple restrictive condition on it. In (32) take j=1, m1, > 2, 
and simplify by aid of (34) ; thus we have 


( 02; )- dc 
In (32) take m= i, j At, n > 2, and simplify by use of (34) ; thus we have 
1 04; @o 1 do do 


In the last equation change j to m, subtract the resulting equation member 
by member from the preceding one, and simplify; then by the aid of (34) 


we have 
2 00 dim 
o Oy: OXm Oy: 
a0 n 0 0x; 
Rae 1 te & @ 02%, 
02m 


the last line being obtained by aid of (31). Hence, since 0%m/dy: is cer- 
tainly not zero for every J, we have by aid of (37) 


(38) 


= 4¢o 2. 
k=1 OL m? 


For determining o there are two cases to be examined. Suppose first 


that c is different from zero. Then from (35), (37), and (38), it follows 
readily that o has the form 


(39) o=c > (%—ax)’*, n>2, 
k=1 


Where aj, ao," * *, a» are constants. Later we shall show that this value of o 


if 

| 

) 

| | 

n 

[ja 
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Suppose * next that c—0. Then from (37) and (35) it follows that 


o is a linear function of 21, %2,°°+*,%n. If we write 
(40) o = fit, + Bot +° + B, n > 2, 


we see from (38) that we must have 


(41) Bi? + Bo? +° Bn? = 0. 


The transformations which arise for such values of o will be considered later, 
In particular we must consider the case when o is a constant, that is, when 


Bo, * * Bn are all zero. 


§ 5. The determination of w. By means of conditions (8) and (9) in 
theorem I or the equivalent conditions (14) and (15) in theorem II it is 
possible to find the value of w explicitly in terms of p. For this purpose we 
multiply both members of (9) by dyi/0x1, sum as to 7 from 1 to n, simplify 
by aid of (21) and reduce; thus we have 


1 0 1 n 0 n 72. i 
wW 021 2p i=1 0x1 k=1 


On simplifying this equation by aid of (27) we have 


Hence except for a trivial constant factor (taken equal to unity) we see from 
this equation and (28) that 


n-2 2-n 


(42) w=p* =o?” ., 


Thus w is completely determined in terms of p except for a constant factor. 
When n =2 we have w—1. 

We shall now show that this value of w may be used to replace equation 
(9) in theorem I and equation (15) in theorem II. For this purpose it is 
sufficient to show that (9) is implied by (8) and (42). In this proof we 
may use equations (21) to (27) since they are all implied by (8). Differen- 
tiating both members of (42) logarithmically we obtain the equation preceding 
(42) ; thence by the aid of (27) we pass to the first equation of this section. 
On multiplying both members of that equation by 0y;/dx; and summing as to 


* The case c = 0 is not treated by Beez in the article already cited. 


| 
1 2—n 4p m—2 
log = — 2p 021 (log p), 1=1,2, 
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iat 1, we obtain an equation which reduces readily to the form of (9). Hence, 


in theorems I and II equations (9) and (15) may each be replaced by (42). 

Equation (42) gives the value of w for the transformation groups in- 
volved in theorems I and II. In the case of theorem III we saw that the 
corresponding quantity w must be a constant, whence it was taken to be unity. 
When n = 2 it is unity even in the cases of theorems I and IJ. When w is 
unity the wave function is an absolute invariant in the usual sense of the term 
in the theory of invariants. The same is true for solutions of equation (13). 
Also we have absolute invariants in all cases when p=1. In general, the 
solutions of equations (1) and (13) are relative invariants, the multiplying 
factor w being a power of the Jacobian J of the transformation, since as we 
see from (20) and .(42) we have 


3 
tw 


(43) w= J 


By aid of the values of o given in (39) and (40) it is seen that w, as 
defined by (42), satisfies equation (10). From this it follows that condi- 
tions (8) and (9) in theorem I are superfluous, except as they are involved 
in finding the value of w in terms of p. Hence the groups involved in theo- 
rems I and III are the same. With the given possible values of o, and hence 
of p, and by aid of (42), it is seen from equation (16) that p must be unity 
in the case of the group in theorem II. 


ae §6. The case when p is a constant C. When p is a constant C, whence 
ao is a constant, the second member of equation (31) is zero. Hence with the 
aid of (23) we see that 


From this it follows that each function y; is linear in the variables 2,, 2», 
*, 2%. Then we must have y; of the form 


Necessary and sufficient conditions that equations (19) are satisfied are then 
the following : 


When these conditions are satisfied we have in (44) transformations belonging 


n | 
s 
y 4 
i 
4 
3 
| 
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to the groups involved in theorems I and III. In order that the transforma- 
tions shall belong to the group involved in theorem II it is necessary to add 
the further condition that C shall be unity, as one sees by the aid of equation 
(16). In view of the result at the end of § 5 it is then seen that the trans- 
formations (44), with the named conditions, constitute the whole of the group 
defined in theorem III. 


§%. The case whenn=2. When n= 2 it follows from equation (42) 
that the function w of theorem I is unity. Hence, when n= 2, the trans- 
formation groups of theorems I and III coincide and each of them is the 
group of conformal transformations in the Euclidean plane. This group is 
sufficiently well known to justify us in not treating it further. The results 
given in the next section, however, are valid for all values of n. 


§ 8. On the class of possible solutions when p is given. From (23), 
(19) and (21) we see that 


(1/p )2 8:;/p. 
Hence the value of p for a given transformation is the reciprocal of its value 


for the inverse transformation. 
Again, if we have a second transformation of the y’s into z2’s we have an 


equation of the form 
Then we have 


2; i 02; OYm 


Hence if two transformations have the characteristic functions p and p then 
the product of the two have the characteristic function pp. 


n n 0z 
02; 02; 
= bij pp- 
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ma- Now let S and 7 be two transformations with the same characteristic 
add function p. Then p™! is the characteristic function of the inverse S of S. 
ion Hence the transformation S17 has the characteristic function 1. If we write 
ns- 

yup S1T =I 


then I is a transformation with the characteristic function 1. But from the 
last equation it follows that T — SJ. Hence if two transformations have the 


9 
¥ same characteristic function p then either one of them is equal to the other 
ha multiplied by a transformation with the characteristic function 1. In § 6 we 
ie have already determined all such transformations and have seen that they are 
Its linear. It follows that if we have a particular transformation associated with 
a given characteristic function p we can determine at once all the transforma- 
tions having the same characteristic function. Therefore it will be sufficient 
), henceforth to determine only a single transformation having a particular char- 
4 acteristic function p. 
§9. The case when o ts a quadratic function and n> 2. The necessary 
q form of o in this case is given by equation (39). It remains to construct a 
particular transformation associated with this value of o. Now we have 


n 0 2 n n 


= (1 1/o* 
( )22 ke Oa, (= + ( ) m= 02m 


km 
n 2 
= (1/s) ¥ ) = 4¢/s, 


the last two members being gotten by aid of (23), (21) and (38). On com- 
paring with the paragraph containing equation (39) we see that constants B 
must exist such that 


= » (yx — Bx)?. 
k=1 


Comparing with (39) we see that we must have F 


(2% — ax)? (yx — Bx)? = 1/c?. 


k=1 k=1 i 


q 

a 

4 
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This suggests the inversion transformation with respect to an n-dimensional 
sphere and hence the solution 
ai 
It is readily verified that these values of the y’s satisfy equation (19) for the 
given value of o. Hence they afford a particular transformation associated 
with the given value of c. The remaining transformations associated with 


t= 1, 2, 


are then determined by means of the results in § 8. 

On taking appropriate translations in the a-space and the y-space re- 
spectively and adopting a convenient magnification in one or both spaces we 
may reduce the preceding transformation to the case in which c 1 and the 
a’s and f’s are all zero. This change can be effected by means of linear 
transformations belonging to our group. Hence, except for suitable linear 
transformations, the preceding transformation may be written in the form 


(47) Yi = where s? = 


This is a real Cremona transformation of period two. It leads us to the follow- 
ing known theorems * (in the i..st case by the aid of the value of w given 
in §5): 


THEOREM IV. Jf w—f(%, %2,°**, Ln) is any solution of the wave 
equation (1), then 
1 


is also a solution. 


THEOREM V. If u=f (2, %,° *,%n) ts any solution of equation (17) 
then 
is also a solution. 


§10. The case when o is a linear function and n > 2. In this case o 
has the form given in equation (40). From (46) and the results at the end 
of § 4, the latter carried over to the inverse transformation, it follows that we 
have 


(48) 1/o = yiY1 + ynYn + y- 


* For the case n = 3 see Kelvin, Liowville’s Jowrnal, 12 (1847): 259. For general 
n see Bocher, Bull. Amer. Math. Soc., 9 (1903): 459. For n=4 see Bateman, Proc. 
Lond. Math. Soc., (2) 7 (1908): 75. The methods of the present paper are different 
from those of these authors. 


ve 


| 
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where 71; y2,° * *» yn are constants such that 
From (40) and (48) it follows that 
(50) (Bit + Bote ++ + B) 
By aid of (23) equation (31) may be put in the form 


Gl 


Ox; 02; 02x; 02; 02; 


In the present case, by aid of (40) and (48), this equation reduces to 


0? 0 0 
+ bun = 0. 


02402; 0x; 


On making the transformation 
(51) yr = 
the equation takes the form 


02 02x; 


y= 0. 
On taking 1 and j different in this equation one sees that ¢, has the form 


where 7';(2,) is a function of x; alone. Moreover, on taking 7 7 we see 
that each T1;(xi) has its second derivative with respect to 2; equal to — y:. 
Hence ¢; can be written in the form 


th = — (21? + +° + En”) + + + + 


where the ¢’s and 7’s are constants. Then from (51) it follows that we must 
have 


2, 


where the constants 8; and y; satisfy conditions (41) and (49). Since it is 


8 
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necessary and sufficient that the y; shall satisfy equations (19), we now impose 
the conditions stated in these equations. Thus we have 


Ox, 


(vite + Bryi—eix) + — jx) 


= > + ying + yini — Yi » + Byi] — [> Breix + Byi), 


k=1 
the last member being gotten by aid of equation (41). Hence, in order that 
the transformation defined in (52) shall be effective, it is necessary and suffi- 
cient that the constants involved shall satisfy the following conditions addi- 
tional to that expressed in equation (41) : 


(53) + ying + = = Bree + Byi = 0, 


It thus appears that condition (49) must be implied by (41) and (53), since 
there was no need of employing (49) in the derivation just given of the 
necessary and sufficient conditions that an effective transformation shall be 
afforded by equations (52). 

We have seen that the y’s must satisfy (49). From (40), (50), (52) we 
see that the following conditions also must be satisfied: 


(54) + yB = 1, + yBi = 0, 1=1, 2, 
k=1 k=1 


It is clear, however, that these conditions on the constants are not independent 
of those previously derived. 

Let us apply the linear transformations of § 6 to simplify transformation 
(52). By a translation of the z-space we may reduce B to zero. Similarly 
by a translation of the y-space we may reduce y to zero. Now at least two of 
the quantities B,, -, Bn are different from zero, since the sum of their 
squares is zero and we are supposing that o is not a constant. Without loss 
of generality we suppose that Bn.s40, Bx40. If we now replace a; by 
V—1 /Bn, 1, we have a linear transformation belonging to 
those in §6. In terms of these variables o has the form 


Bits’ + Bow! + Bn-10'n-1 V— - 


where 


Bi? + Bo? + = 1. 


n n 
n n 


t 


Certain Partial Differential Equations of Physics. 115 


Now if we revolve the (n —1)-space 2’n_, 0 in such a way as to leave the 
(n—1)-space 2,’ 0 invariant, we can bring the former into coincidence 
with the (n—1)-space Biya’ +° Let 2”, 
av’'n-1) Un’’ =n’ be the new coordinates after this revolution; then o has the 
form o = V—1 aa”. The transformations employed all belong to 
the group of § 6. Hence, by means of transformations belonging to the linear 
group of § 6, we can reduce o to a normal form in which the coefficients B are 


Bri=1, 

Let us now take the original value of o in this normal form. In view of 
the results in § 8 it is sufficient to obtain a single transformation of the type 
of (52). In order to find a convenient one we may also similarly normalize 
the y-coordinates so that y= y: = y2="** = yn-2 = 0, yn-a = 1, yn = V—1. 
Then it is easy to obtain a particular solution of equations (53). Thus we 
have for a particular form of transformation (52) the following: 


‘a Vi 


i= 
+ V—1 & 
1— 3? 
1+ s? 
where s? 2,2 -+ -+- 2,7. For this transformation and the resulting 


theorems in the special case when n = 4 see Bateman, I. c., p. 77. 
Then we have the following theorems: 


VI. If to, Xn) ts any solution of the wave 
equation (1) then 
U = + V—1 (Yr, Yor Yn) 
is also a solution, the functions y; being defined by equations (55). 


VII. If Un) ts any solution of equation 
(17) then u=f(y1, °° *, Yn) also a solution, the functions y; being 
defined by equations (55). 


The transformation defined by equations (55) is an imaginary Cremona 
transformation of period two. A real Cremona transformation of period two 
in 2n variables may be obtained from it by separating the given variables into 


| 
Byi], 
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real and imaginary parts and equating the real and imaginary parts of the two 
members of each of the resulting equations. 


§11. Summary of results. Some results proved but not explicitly stated 
in the foregoing sections may be summarized in the form of the following 


theorems: 


THEOREM VIII. The transformation group involved in Theorem I and 
that involved in theorem III are the same. The group is that of the conformal 
transformations in space of n dimensions characterized at the end of §3, 
imaginary transformations being allowed. The transformation group of theo- 
rem II is that subgroup of the foregoing in which p=1. (It is completely 
characterized in § 6.) 


THEOREM IX. The conformal transformation group of theorem I or III 
(see end of § 3), when n > 2, is generated by the linear transformations of § 6 
and the two particular transformations (47) and (55). 


By means of the transformation indicated in the second paragraph of the 
introduction, it is seen that the two particular transformations (47) and (55) 
are real when they are taken in the form applicable to the wave equation as it 
is generally written by the mathematical physicists. 


Transformations Leaving Invariant the Heat 
Equation of Physics.* 


v0 
By JouN A. GOFF. 
d 
IS We shall write the general heat equation in the form 
dv 
d (1) kaa 
For the case when n = 3, the transformation y, cpt/k, y1= 2, Y2o=Y, 
3, y; = reduces equation (1) to the classical form 
)- 
dv fu Fo Fe 
To of the heat equation of theoretical physics. 
6 §1. The Partial Differential Equations Defining the Transformation 
of the Heat Equation into Itself. Let the function v(Yo, * * Yn) 
e be a solution of the heat equation (1) and consider the transformation, 
) (2) U= W°V(Yo, °° Yn)s w 
where w is a funcion of 2%, +, Also the functions yo, 1, °° *, Yn 
are to be functionally independent functions of Xo, 2%, ° * Yn so determined, 
in connection with w, that whenever v satisfies equation (1), wu will satisfy 
the following equation: 
du 
We shall prove the following theorem: 
THEOREM I. Let w, yo, Yn be a set of functions of %1,° Xn 
such that yo, Yn are functionally independent and Then a 
necessary and sufficient condition that every solution of (1) shall be trans- : 


formed through (2) into a solution of (3) is that the functions w, yo, Y1; 
"++ 5 Yn shall satisfy the following conditions: a function p exists, not iden- 
tically zero, such that 


Ow Ow 
4 
( ) 02 k=1 02%," 


* Written with the aid of the foregoing paper by Professor Carmichael. Refer- 
ences to Carmichael’s paper are indicated by prefixing the letter C; thus, (C30) refers 
to Carmichael’s equation (30). i 
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OY; 

5 

(6) dij P> J — 2, nN, 
(7) = 02,” 2 0x, 7 i, 


k=1 
Since 1, yj, (1547), and yj? — yn? (1, J = 1, 2,°° +, m) are solu- 
tions of (1), it follows that the functions w and y; must be such that w, wy;, 
wyiys (1547), and w(y;? — yn?) for 1, —1, 2,- are solutions of (38). 
These facts imply equations (4), (5), and (6). Substituting the value of 
from (2) into (3) and simplifying by aid of (1), (4), (5), (6), we have 


dv 


n dYo 2 

The particular solutions of (1), namely, y,—4y;* and yoy; — %y;* 
(j=1, 2,---+, ), show that (7) must be satisfied and that (8) must be 
satisfied for 71, 2,---+, m. That (8) must also be satisfied for 7 —0 
then follows immediately. The conditions of the theorem, with the exception 
of the condition that p be not identically zero, are thus shown to be necessary. 

Formulae (C19) and (C21) can be used since both of these are implied 
by (6). Let us multiply (8) by dy;/dy:, sum as to j from 1 to n, and sim- 
plify by means of (C21) getting 


Then (7) becomes 


dYo 
(10) 


The Jacobian (functional determinant) is 


j, 1,2,° °°, 0, 


which becomes, by virtue of (9) and (10), 


by; | 

J = | 


the Heat Equation of Physics. 


Oy 
Ox; 


J =p j, 1, 0, 


so that 


N+2 
(11) J =p? by (C20) 
Consequently it is necessary and sufficient that p be not identically zero in 
order that the functions yo, y1,° * *, Yn be functionally independent functions 
of 2%, %1,° °°, Xn. Moreover conditions (4) to (8) can easily be seen to be 
sufficient conditions. Hence the theorem. 


§2. Determination of w, Yo; 41,° °°; Yn. From equations (9) and (10) 
we get 


so that p is a function of x, alone. If we consider the equations 
= Vi (Yo, °° Yn)s t= 1, 
of the inverse transformation, then on differentiating with respect to x; we 


get 
n 02x; 


k=0 OY, 02; 
Then by virtue of (9) we have 
0x; 
Hence we can use equations (C31) and (C23) together with (12) to show that 
(13) 0, j, k, 1, 2, % 


The functions y; are then of the form 


where the a’s are functions of 2, alone. 
In view of (13), equation (5) becomes 


Oyj I(log w) dys 
Now (6) implies that 
n @ Ou; 
(C21) dip, k, t= 1, 2,°° 


Differentiate both sides of this equation with respect to 7, and then by means 
of (15) obtain 
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@(log w) Oyj dys » @(log w) dy; Oy; 
2 j=1 Oxy 2 02 jr (0m 


When this is simplified by means of (C21) and (12), the result is 


é?(log w) 


—1%; = —1% 1, 2,° °°, 
02,021 4 Okl p’/p 4 


(16) 


where p’ denotes the derivative of p(x.) with respect to its argument. For 
convenience we have denoted the ratio p’/p by &. 
From (15) with the aid of (13), (14) and (16), we get 


0?y; 


= $f-a;1 2, 


But from (14) alone we can obtain directly, 


where the prime denotes the derivative with respect to xz. Consequently 
a’ 31 p’/p, whence 
(17) aj = Cy pl, 0, 


where the quantities Cj; are constants. These results show that the y; can 
be written in the form, 


n 
(18) yi = + yi} 
m= 
where, by virtue of (6) and (C21) the constants satisfy the conditions 
k=1 


and the y; are functions of 2, alone. 
From (16) we see that log w is of the form 


n 


(20) log w = — (tm — Bm)? 


where the functions B are functions of 2, alone. Moreover, since w must 
satisfy (4), log w must satisfy the relation 


From this, with the aid of (20), we obtain the an equation : 


k, t==1, 2, 


the Heat Equation of Physics. 


(22) (em —Bn) — (n/4) 
+ (1/1632) Bn)? =0. 


This equation implies that 


(23) 20’ and = (n/4) 
and, if that 
(24) B’m = 0, % 


Now form 0z;/dx, from equation (18) and compare with @y;/0x, ob- 
tained from (15) with the aid of (18) and (20). This leads to the following 
differential equation : 


(25) Vi + + ~ CimBm = 0, j= 1, 2, pe 
m= 


Remembering that = p’/p and supposing for the moment that p is not a 
constant, we then see, by aid of (24), that 


(26) CimBm + C;/p%, 2,:°-,n 


where the C; are constants. In this case the Bm are also constants. Conse- 
quently (18) becomes 


when p is not a constant. 


§3. The Case when £=0. In this case p=k?=constant. From 
(25) we see that each y; is a constant which we may choose in such a way as 
to enable us to write 


(28) = k ~ Cim (Xm Bm) j=l, 2, 
m= 
where the constants Cjm must satisfy conditions (19); namely, 
n 
(19) = CimCim = 
m=1 


From equation (10) we see that y, is of the form 


(29) Yo = k? (to — Bo). 
Furthermore from (2) and the last of (23) (which becomes ~’ —0) we see 
that w is a constant and without loss of generality we may take 


(30) w= 1. 


| 
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If pt = 1, we may write (28) and (29) as follows 


9 


Cjm(%m — Bm); j=0, * N, 
m=0 


where Com = Cmo = Then 


t,j=0, 1, a, ° 0; 


— bu 


for if 70 and 1-40, it reduces to 0; if 70 and 1—0, it reduces to 1; 
while if 1540 and 7=£0 it becomes simply 8; by (6). These results show 
that the transformation described by equations (28), (29), (19), and (30) 


is conformal when p? = 1. 
§4. The Case when £40. By using equations (23) we get the fol- 
lowing particular solutions: 


(31) = — 2/2, B=log 2/2); 


whence by aid of (20) and (27) we have 


2 2 2 
(32) w = exp + 22 + tn \ 
0 


By (10) we have 
(33) j= 1/2. 
The constants which have been omitted in obtaining these particular solutions 
can be removed by a transformation described under § 3. Thus we can state 
the following theorem : * 
THEOREM II. If u—f(%, %, %n) is any solution of the heat 


equation (3), then 


420 J 
f (—1/2, * Ln/Z) 


is also a solution of (3). 


§ 5. Summary. The results show that the transformation of theorem 
II together with the transformations described by equations (28), (29), (19), 
and (30) generate all the transformations of the form of equation (2) which 
leave the heat equation invariant. 


* For the case n=3 see J. Brill, Messenger of Mathematics, 2 (1891), p-. 137; 
for the case m= 1 see P. Appell, Journal de Mathématiques, (4) 8 (1892), p. 187. 
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The Application of Fractional Operators to 
Functional Equations." 
By H. T. Davis. 


1. Fractional Operations. In a recent paper ¢ the author showed how 
it was possible to express certain types of integral equations as differential 
equations of fractional order and by operating on the latter to obtain solu- 
tions of the original equations. Since that time several examples of such 
equations have come to his attention, but the methods of solution employed 
in these cases made no use of the idea of fractional operations. For example, 
John R. Carson in a series of papers { has made an analysis of the calculus 
which Oliver Heaviside developed and used so effectively in the circuit prob- 
lems of electrical theory. Carson shows very elegantly how the Heaviside 
calculus is formally equivalent to the process of solving an integral equation 


of the form 


but he does not point out the interesting fact that the most striking feature 
of the Heaviside calculus is also a special application of the fractional operators 
of Riemann.|| In particular, the series expansion of the Heaviside operator 
for the general case of the non-inductive cable which Carson characterizes as 
“hopelessly complicated to either interpret or compute ” is reduced with con- 
siderable simplicity by the methods of this paper.{ 


* Presented to the Indiana Section of the Mathematical Association of America, 
May 8, 1926. 

t “ Fractional Operations as Applied to a Class of Volterra Integral Equations,” 
American Journal of Mathematics, Vol. 46 (1924), pp. 95-109. 

t+“ Theory of the Transient Oscillations of Electrical Networks and Transmission 
Systems,” Trans. A. I. H. E. (1919); “ The Heaviside Operational Calculus,” The Bell 
System Tech. Journal, (1922) ; “ Electric Circuit Theory and the Operational Calculus,” 
Ibid., Vol. 4 (1925), pp. 685-761; Vol. 5 (1926), pp. 50-95; Bull. Amer. Math. Soc., 
Vol. 32 (1926), pp. 43-63. 

§ Electromagnetic Theory, Vol. II. 

|| A similar observation has recently been made by P. Levy: “Le calcul symbolique 
d’Heaviside,” Bulletin des sciences Mathématiques, Vol. 50, 2nd series (June, 1926), 
pp. 174-192. 

{ Bell System Tech. Journal, Vol. 5 (1926), p. 59. 123 
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It seems important, therefore, that some effort should be made to give 
wider introduction to these fractional operators and to suggest that a standard 
symbolism be adopted for them. The great elegance that can be secured by 
proper use of fractional operators and the power which they have in the solu- 
tion of complicated functional equations should more than justify a more 
general recognition and use. The aim of the present paper is to present the 
case of the “ fractional calculus.” 

The idea of fractional differentiation was clearly in the mind of Leibnitz 
who said in a letter to G. F. A. de Hospital in 1695 that he was considering 
the meaning to be attached to d?/dz*. The first real use of fractional opera- 
tions was apparently made by Abel who, as early as 1823 used both the sym- 
bols d4/da4 y(x) =1/Vx fiy(a) dx and y(x) = Vx in solving 
the problem of the tautochrone.* 

Although there is mention of fractional operators in the works of Euler, 
Laplace and Fourier, it was J. Liouville who first made an extensive study 
of the new calculus. Altogether he devoted eight memoirs totaling about three 
hundred pages to the subject making application to problems in geometry and 
mechanics.¢ The effectiveness of this work was somewhat limited, however, 
by the fact that a specialized definition of fractional operations was employed. 
Most of the modern work in this field goes back to Riemann who, in a paper 
developed during his student days but posthumously published, gave general 
definitions of fractional differentiation and integration.{ 

Slightly generalized these definitions may be stated as follows: 


By the fractional integration of a function u(x) we shall mean the 


operation 


u(t) dt, v>0, (1) 


and by fractional differentiation, 


(a—t)” 


ome 
dam1 u(t) dt, 1, 2, 


<1, (2) 


(x) 


* “Solution de quelques problémes a l’aide d’integrales définies.” Werke (1881), 
Vol. 1, pp. 10-27. 

Journal de Vécole polytechnique, cah. 21 (1832), pp. 1-69; 71-162; 163-186; cah. 
24 (1835), pp. 17-54; cah. 25 (1836), pp. 56-84; Journal fiir Mathematik, Vol. 11 
(1834), pp. 1-19; Vol. 12 (1834), pp. 273-287; Vol. 13 (1835), pp. 219-232. 

¢“ Versuch einer allgemeinen Auffassung der Integration und Differentiation.” 
Gesammelte Werke, Leipzig (1892), pp. 353-366. 
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It is important to note that the constant “c”.4s a vital part of the opera- 
tional symbol. Thus Riemann let c= 0 in his original definition and Liou- 


ville let c—=— ©. 
Various symbols have been suggested for these operators. We have 


already mentioned d’/dx” and d”/da’, but they are somewhat cumbersome 
and do not convey information regarding the lower limit of the integral. 
This symbolism was effectively used by Liouville who, as has been stated, 
employed a fixed limit. Emil Post in discussing the solution of the problem 
dy /dat — y/x =0 in the Mathematical Monthly several years ago * suggested 
the use of the symbol {d-’/da-’},*. Riemann used 0,” and 0,” for his opera- 
tors and P. Levy in a recent paper + employed the symbols J,”’[¢()] for 
integration and D,”[¢(«)] for differentiation where the lower limit of the 
integral was zero. In a contribution to the subject by H. Laurent in 1883 


x 
we find the operation expressed by at .[ J. Hadamard, who made funda- 


mental use of Riemann’s operator and gave some generalization to it in his 
work on analytic functions, employed D,”’¢(2).§ The mathematical encyclo- 
pedia, to which the reader is referred for a longer bibliography on the subject, 
uses D,”” f(a) as a symbol for the general operator.|| Since all of the informa- 
tion is conveyed in a simple manner by the symbols -D,”u(x) and -Dz”u(z) 
and since they are only slight variations of those which appear to have had 
the largest usage, the author has adopted them in his work. It is occasionally 
convenient to use the easily understood abbreviations u(x) and wu” (2) 
when there is no question as to the lower limit of integration. 


2. Some Properties of the Operators. For convenience in reference a 
few of the properties of fractional operators will be stated. Where proofs 
have been omitted they will be found elsewhere. 

The significance of the definitions given in the preceding section attaches 
to the fact that 

+ v(x) |] = u(x) + De“ v(x), 
cDr"k u(x) = k u(x), (3) 


* Mathematical Monthly, Vol. 26 (1919), p. 37; see also Vol. 25 (1918), p. 172, 
problem 433. 

+ Bulletin des Sciences Mathématiques, Vol. 47 (1923), pp. 307 and 343. 

+ Nouvelles Annales, series 3, Vol. 3 (1883), p. 240. 

§ Journal de mathématiques, Vol. 57 (1892), p. 154. 

|| Encyklopidie der Mathematischen Wissenschaften, II, A, 1, pp. 116-119; II, A, 
1l, pp. 770-772. Also Encyclopédie des Sciences Mathématiques, T. 11, V. 5, fas. 1, 8. 
{ Davis, loc. cit. 
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The analogue of Taylor’s theorem is found in the following development 
of the integration operator: 


This is easily proved by replacing u(t) in (1) by its expansion about the 
point z. We thus get, 


dt, 


De” u(r) = 


which easily reduces to (4). 
Another useful expansion is that analogous to the formula of Leibnitz 
which may be stated as follows: 


(2) V(2)] =U (2) (2) Vr” (2) + 
(0) (2) — (5 


In the application of these formulas it is often useful to note the fol- 
lowing derivative : 


ds, we then have 


v(t) = c) +f" Dag (2, s) u(s) ds 
+f # — u’(s) ds. (6) 


This formula can be established by the method of “ partie finie” due 
simultaneously to R. D’Adhemar and J. Hadamard in 1904.* 


Let us write 
_ (7* 8) u(s) 
since there is no longer a singularity in the integrand for the upper limit of 
the integral, we can apply the ordinary rule for the differentiation of a func- 
tion defined by an integral. We then have 
{ 0/dx g(x,8) u(s) 
— + 9(2,s)u(s) (x—s)” \ ds 
+ 
*R. D’Adhemar, Thesis, Journal de Math. (1904); see also Hwercices et lecons 
d’analyse, Paris (1908), p. 150; J. Hadamard, Congres de Mathém. (1904), pp. 265-271. 


g(z,x—e) 


€ 


it 
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Since 0/da (x —s)-” =— (x—)-”, we can substitute this value 
in the integrand of the preceding integral and integrate by parts. We then 


get, 


_ 0/0x 8) u(s) + 0/08 [g(a,s) u(s) ] g(z,c) u(c) 


dx (~x—s)” 
and in the limit we have formula (6). 
It will be useful in another section of the paper to know that the well- 
known operator: * 


(d/dx + k)" f(x) d"/da" f(x)}, ma positive integer, 


is also valid for fractional values of n. This can be proved by means of for- 
mula (5). 
We first write, 


ete dn/dan (ot f(x)} (oto f(2)}, 


where m is an integer and p» a fraction between 0 and 1. Then in (5), setting 
U(x) equal to e** and V(x) equal to f(x), we get, 


{el f(x) } = 
eka qm/ dam [ (a2) — + = 1) 


eke dm /dam [eke (1 — phat + kee? (2x) ], 


where we have used the abbreviation z= d-1/dz1._ This expression then re- 
duces to 


ete dm/dam [ek (14 (2)] = 
(d/de + k)™ (1 f(x) = (d/de + k)™# f(a). 
3. The Solution of Fractional Equations. One of the simplest cases 
of fractional equations is that of the linear equation with constant coefficients 


involving derivatives (or integrals) of orders which differ by multiples of 1/n, 
n an integer, from one another. Thus let us consider the equation, 


DL (u) = Aou + + +--+ Ag f(z), (7%) 


* J. Edwards, Differential Calculus, 3rd ed., London (1904), p. 70. 
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where we have designated by u the derivative -D,”u(#). The function f(z) 
will be characterized later. 

If we replace u(x) by d/dx f?u(s) ds and the fractional derivatives by 
their formal values we see that (7) is equivalent to the following functional 
equation : 


G (a, s) u(s) ds = f(z), 


dx ‘ (a—s) 1-1/n 
where 
A 
G(x, s) Ay (a S) + r(1—1/n) (x s) + + T(1/n) 


Upon integrating both sides we arrive at the integral equation: 


The existence theorem is well known and has been stated under a variety 
of conditions for f(x).* It is sufficient to assume that f(z) is analytic in the 
neighborhood of ac. Since the first equation is derived from the latter 
by differentiation, it is clear that a solution of the integral equation will also 
be a solution of its derivative. 

Having thus assured ourselves of the existence of at least one solution of 
equation (7%) we can proceed to find this solution by the following formal 
method. This method will be found to apply equally well to an equation 
involving fractional integrals differing by multiples of 1/n from one another 
and since the details are essentially the same no separate discussion will be 
made of this case. 

Representing (7) symbolically by P, 


P=A, + A,zi/n A,2?/n “+ 


we can rationalize this algebraic equation in the following manner. Multi- 
plying P successively by z!/", 2?/", etc. we get the system of equations, 
P = Ao + 4+ 4+ An 


gi-1/nP — ZA, zAozi/n + 2A,22/" + 


*T. Lalesco, Journal de Mathématique, Vol. 73 (1908), pp. 125-202. 
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Eliminating 21/", z?/",- - - from the right hand side of these equations 


we obtain 
Ay — P, A2, An-1 
ZAn-1 gi/nP, Ao, A, baw 
zAn-2 z2/nP, Ap, Aus = 0. 
ZA, gi-i/np, 2Ao, 2Asz, Ao 
From this it follows that 
Pa(s) =A(2), 
where Ai, Ao, 
A (z) g2/n, 2A n-19 Ao, A (S) 
gi-1/n, zAo, 2A3, A, 


and A(z) is the same determinant with the first column replaced by Ao, 2An-1 


Since A(z) is a polynomial in z of degree less than or equal to n—1, 
it is clear that A is a rationalizing factor and equation (7) can be formally 


replaced by the ordinary linear equation 
A(d/dz) u= A(d/dz) f(z). (9) 
To this equation there corresponds an integral of the form 
u(x) = W(a, t) A(d/dt) f(t) dt, (10) 


where W(z,t) is the Cauchy function familiar from the method of variation 
of parameters, provided W(x,t) A(d/dt) f(t) does not become infinite to a 
higher order than p< 1 at 

To show that (10) is the unique solution of (7) we shall need the fol- 


lowing theorem : 


THEOREM 1. The homogeneous equation L(u) =0 has no solution ex- 
cept the trivial one u=0, provided cA. 


9 
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Proof: Equation (7) is equivalent to the equation 
A,(a@— s)-2/™ 
Ta—i/n) * Tad—2/n) 


> -(1-1/n) 
\ u(s) ds =0, (11) 


which may be written from (6) in the form 


L(u) = Aou(c) I(x, 8) w’(s) ds + u(c) I(2,¢) =0, 


where I(z,s) is the integrand of (11). Now unless u(c) =0 we should have 
a contradiction since lim Z(u) would then be infinity instead of zero. Simi- 

larly we may prove that u(c) = u’(c) =- (c) =0, which would not 
be the case unless ¢ were a singular point of the equation A(d/dx) =0. But 
since c= oo is the only singular point we conclude that Z(w) —0 can have 
no solution except the trivial one. 

As illustration we shall consider two special equations which are of fre- 
quent occurrence in application. 

I. As a special case of (7) we have 


u+ u=f(z). (12) 
Rationalizing we obtain the ordinary equation 
d/dxu—*u = f(x) —Af(z), (13) 


whose solution is 
u(x) =e { f(t)}o™ dt. (14) 


In order to see that this solution behaves properly at the point «—c, 
write (12) in the expanded form 


u(c) u’(t) 


‘ If f(x) is bounded at the point x—c, then it is clear that u(c) must 
be zero if we are to avoid a contradiction. Since (14) is the only solution 
of (13) which satisfies this condition it must be the solution of (12). 


Example: Solve .DeAu—u=—1. 


We get as the formal equivalent of this equation the ordinary equation 


u’ —u—=— (1+ 1/Vze), 
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whose solution is 


u(x) =e" f 2 (—-1— 1/V zt) dt, 


This function can be shown by direct substitution to satisfy the original 


equation. 


II. Consider also 


cDztu + u =f (2), (15) 
whose rational equivalent is the equation 
d/dx u—u=dAf’ (x) — f(z). (16) 


The solution of this equation is 
u(z) = Cea + cDaf(t)} dt. (17) 
In order to determine C consider equation (15) in the form 
— u(t) 
1/Vr dt Au(z) == f(z). 


It is clear that for ec we have Au(c) —f(c). Then from (17) we 
get u(c) =C e1/”*¢ —1/)f(c) and (17) can be written in the form 


u(x) =1/d + faa f’(t) Dt f(t)} dt (18) 
Example: Solve 


Detutu=—i. 
This rationalizes to the equation 
1 
V ra 


whose solution is 
1] 
u(x) =e? —e* f et dt. 
0 


This solution can be expanded in two ways: 
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The second, although everywhere divergent, is the asymptotic expansion 
of the first. 


4. Generalization of Abel’s Integral Equation. Among special types of 
integral equations we find the following: 


kat = a Va—t u(t) dt, (19) 


which was first solved by N. Hirakawa* and subsequently discussed by T. | 
Hayashi ¢ on the basis of an elegant method due to P. Humbert. 
The generalization 


= (27 —t)™ u(t) dt, 
m=0,1,2,°°: 


which includes both Hirakawa’s and Abel’s equations § as special cases is easily 
solved by fractional operators. This can be written 


f(x) +v +1) De 


Operating on both sides with the symbol -D,™*”*! we have 


(x) 
dm2 | 


Recalling elementary properties of the gamma function and making use 
of (6) we can write (21) as 
1 {(—1)m1 T(m+v+1) f(c) 
T'(m +v+1) T(1—v) T'(v) (x—c)™™! 


m T(m+v) _ f’(c) 
I'(v) (7—c)™” 


u(r) = 


T'(v) -(1-v) f(m+2) 


* The Téhoku Math. Journal. Vol. 8 (1915), pp. 38-41. 

+ The Téhoku Math. Journal, Vol. 10 (1916), pp. 56-59. 
t Proc. Edinburgh Math. Soc., Vol. 32 (1914), pp. 19-29; Vol. 33 (1915), pp. 35-41. 
§ Abel’s equation is the special case m = 0, 0. 


) 


m+v+1 


22) 
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But it is clear that if the value of u(x) given by equation (22) is sub- 
stituted in (20), the integral will not be convergent unless 


Under these conditions u(x) defined by (22) will be the desired solution. 
These results can be summarized as follows: 


THEOREM 2. If in equation (20), f(x) eaists together with its first 
m+ 2 derivatives in the interval cS %, and tf f(c) =" 
f™(c) =0, f° (c) 0, then the solution of equation (20) exists in the 
open interval c< 22%. The solution is given by formula (22). 


We can apply these results to equation (19), assuming a to be a positive 
integer, and obtain the solution with the greatest ease. Since we have 
T'(a—1) 


*/2, it follows that, for a > 1, 


Det a(a— 1) =a(a— 1)k 


T(a+ 1) a-8/2 
and fora=—1, 
k 1 2k 1 
u(z) = — 
T(%)T(4) V2 Vea 


5. Lalesco’s Integral Equation. Lalesco’s integral equation is a gener- 


alization of the problem that we have just solved.* 
Thus let us consider 


f(z) u(t) dt, (23) 


where we assume that 


G(x, t) = go(t) + gi(t) (e@—t) ++ gn(t) 
It is then apparent that (23) can be written in the form 
f(z) =T(1—v) (gou) +1(2—v) De (24) 
We can now operate on both sides of equation (24) with the operator 
eD,"*1-¥ and obtain the following ordinary differential equation of nth order: 
T(1—v) d*/dx" (gu) + T(2—v) (gu) 
(n-+1—v) guru = (25) 


*T. Lalesco, “Sur l’équation de Volterra,” Journal de Mathématique, Vol. 73 
(1908), pp. 125-202. In particular §9 of part I. 
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Since equation (25) is a linear non-homogeneous one, its solution will be 
given by the method of variation of parameters in the form 


— > (2) + f t) f(t) dt, (26) 


where w,(2),* °°, Un(x) are n linearly independent solutions of the homo- 
geneous equation, W(a,t) is the Cauchy function, and the c; are arbitrary 


constants.* 
Returning to equation (24) let us operate on both sides with .D,-’. We 


shall then have 

f(x) =T(1—v) +1(2—v) Dat + 
T'(3—v) Da? (g2u) +° +0 (gnu). 

Then if (x) = 0, we must have 

go(c) u(c) =0. 

Operating again with .D,! = d/dz we get in similar fashion 
[go(x) w(x) + (1—v) gi(c) u(c) = 0. 

Repeating this n times we have as boundary conditions the system: 

go(c) u(c) =0, 

[9o(x) u(x) + (1—v) gi(c) u(c) 


[9o(4) u(x) + (1—v) + (1 —v) (2 —v) g2(c) u(c) =0, 
(27) 


+ (1—v) [gu] 
(n—1—v) gu (c) u(c) =0. 


The determinant of this linear system in u(c), u’(c),° (c) is 
seen to be equal to [go(c)]". Hence if go(c) £0, we shall have u(c) = 


* An important example is found in Whittaker’s treatment of Abel’s integral 
equation: “On the Numerical Solution of Integral Equations,” Proc. Royal Soc. of 
London, Vol. 94(A) (1917-18), pp. 367-383. In Whittaker’s paper the g;(t) are con- 
stants, If the roots of F(z) = (1—a)a,z"-1 + (l—a) 
+(1—a) (2—a)... (m—a)a, =0, are A, B,..., N, the Cauchy function will be 

eA(a-t) eB(a-t) eN (a2-t) 
By simple considerations Whittaker’s form of solution is then obtained from (26). 
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be u’(c) =" * *=u%l(c) =0, and it follows that in (26) we get after n 
successive differentiations : 
n 
(c) =0, (j= 0, n—1). 
0- Since the determinant of this system is the Wronskian, W = (uw, w’2,°°', i 
un»), which is different from zero for x = c¢, provided g,(c) ~ 0, it follows 
that c, = Co * *=Cn=0. In this case the solution of (23) ts unique. 
We next turn to the problem of solving (23) when go(z) 


= gr-1(%) =0, gr(c) #0. In this case we replace equation (24) by the 
following : 
+1—v) [g-(z) u(x) ] 
(r+ 2—v) [gris (x) u(x) 
Do” [gn(x) u(x) ], 


and equations (27) by the system: 


gr(c) u(c) = 0, 
[gr(x) w(t) + (r +1—v) w(c) =0, 


[gr(a) w(x) ] (r+ 1—v) [gras] 
gn-1(¢) u(c) =0. 
By reasoning similar to that in the first case it is clear that equation (23) 
will have a unique solution if g,(c) 0. 
It is next of interest to consider the case where 


Go(%) = (t—6)"Go(z), gu = 
Gn-1(%) = (2). 


It is further assumed that Go(c), Gi(c),° ++, Gnr(c), gn(c) are all 
different from zero. 

Then system (27) is satisfied for any value of u(c), u’(c),: 
whatsoever. The solution of (23) is therefore the general one (26) with n 
arbitrary constants. 

Two examples follow: 


Ajtu(t) 


Baample 1. f(2)— at + dt. 


This can be reduced to | 
A, (au) + Ao T(%)u = f(z). 


| 
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From A,d/dz (xii) we get where cm 
— (A, +44:2)/A:. Hence we have 


u(x) + 1/A, (x/t)° f(t) dt, 
provided the integral exists and c > —1. 


Example 2. Consider the homogeneous equation 
2 
A; f u(t) dt+ Az f t (x—t)1" u(t) dt+ A; 
0 
(x—t)2" u(t) dt =0, 
0 


which reduces to 
(x) + a,ru’(z) + a.u(x) = 0 
where = Aj, = 4A, + A2(1— 1), = 2A, + A2(1 
+ A;(2—v)(1—»). 
The solution will be 
u(x) = + 
where 7, and 7, are roots of the equation 
aor(r—1) +ar+a,—0. 
It is easily shown that this function satisfies the original equation pro- 
vided r,; and r, are greater than — 1. 
The case of equal roots presents more difficulties. Then the solution will 
be u(x) = 2" log x where r is a root of the two equations 
ar(r—1) +ar+a.—0, 
aor + +a,—0. 
To verify that w(x) is indeed a solution consider the integral 


@ 
I(z) = f log t — t)# dt. 
0 
Making the transformation t= zs, we shall have 
I(2) = f (log « + log s) a (1—s)4“ds = 
7 0 


TH+) 
J) 
— {loge B(A-+ 1, w+ 1) +B(A+1, p-+1) 
[y(a +1) 


{ log 
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Substituting log in the original equation and noting that 
A,B(r + 3, i-~#} + A,B(r + 2, 2—v) + A;B(r + 1, 3—rv) = 0, we get 
for the left hand side the expression 

A,B(r + 8, 1—v)[y(r + 8) —¥(r+4—y)] 

+ +2, 2—v)[¥(r +2) —y(r+4—»)] 
A,B(r +1, 3—v)[¥(r+1) 

which we must show is identically zero. 

If we make use of the fact that y(r +2) =y(r+1) +1/(r+1), 
y(r+3) =y¥(r +1) +1/(r+1) +1/(r +2), it will be seen that the 
above expression reduces to 


Ai(r +2) + Ai(r+1) + 42(v—1). 


This is readily seen to be equivalent to aor + a,(r—1) + a, which is 


equal to zero. 


6. Sbrana’s Integral Equation. Another type of equation which yields 
easily to the calculus of fractional operators is that of F. Sbrana who, how- 
ever, used other methods in his discussion.* 

Sbrana’s equation is of the form 


f(z) = f u(t) K(«—t) dt, (28) 
h K 2 ecu 
where * 0<a<l. 


If we let z—=v—t, and u—2—+t, we shall have 


(v—t)* 


If this is substituted in equation (28) and the order of integration in- 
verted by Dirichlet’s formula, we shall have 


é (c—v)? dv dt = f(z). 


This equation can be written in the form 


*“ Sopra certi equazioni integrali di Volterra, risolubili con procedimenti finiti,” 
Atti dei Lincei, Vol. 32 (1) (1923), pp. 156-157. 
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By successive integrations we obtain the solution 


As a special case the kernel can be a Bessel’s function since 


e7t(z-2u) 
Jo(z) =1 — du.* 


It is easily proved that similar treatment can be made if the kernel is of 


the general form 


K 2 een d Un-1 een-1 (tn-1-Un-2) 
( 0 (2 Un-1)% (Un-1 — Un-2)* 


%. Applications in Electrical Theory—Heaviside’s Calculus. The fol- 
lowing applications to the flow of electricity in a cable are considered both by 


Heaviside and Carson although the methods which they employ do not ex- 
plicitly make use of Riemann operators. 


Consider the case of the flow of electricity in a non-inductive cable with 
distributed resistance R and capacity C per unit length subject to an impressed 
voltage V(t) at the point 0. 

The differential equations of the cable are + 


RI = — V, 
(0/at) V =— (8/dz) I, 


where x is the distance measured along the cable from a fixed point («= 0), 
I is the current at the point z, and V the corresponding potential. 

Using a method that is familiar in the theory of linear partial differential 
equations,{ we can replace 0/0t by z and by elimination of J obtain the equation 


2RCV = #V/d2?. 
Symbolically we shall then have 
V (a, t) =e V,(t) + V,(t), 


*Riemann-Weber, Partiellen Differentialgleichungen, (Braunschweig), (1900), 
Vol. 1, § 73; (1919), § 75. 

t Starling, Electricity and Magnetism (1916), p. 449; Jeans, Mathematical Theory 
of Electricity and Magnetism (1915), p. 333. 
t Forsyth, Differential Equations, 5th ed. (1921), pp. 522-23. 


| 
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where a= Vz RO and V;(t) and V.(t) are arbitrary functions. 
By differentiation we get 


I(2,t) =V2zC/R V;— Vo). 
‘Assuming that the cable is infinitely long so that the reflected wave is 
absent we may set V,= 0 and have as the symbolical solution of our problem 
I=V2C/R e* V,(%), 


where V,(¢) is to be determined from the boundary condition. Since this 
asserts that V(x,t) | c- = V(t), it is clear that Vi(¢) == Vo(t). Hence the 
complete solution of the problem is found in the expansion: 


I(z,t) = V2z0/R V,(t), 


2 Az)* 


. 
— VC/R {2Ax + Vo(é), 
where for brevity we have written A = VRC. Therefore we get 


If Vo(t) =1, that is to say if a unit E. M. F. is impressed at c= 0, 
we get the well-known solution,* 
RC2? 
I(z,t) = * 
If Vo(t) =sin ot, that is to say if an alternating E. M. F. is impressed 
on the circuit, we have at x =0, 


Kon d t SiN ws 
1(0,t) = 


ds. 


C 0s wS 
Rr 0 (t—s)3 


In establishing the convergence of the series representing I(a,¢) in the 
general case for real values of ¢ greater than 0, it will be sufficient to assume 


* Jeans, loc. cit., p. 334. 


f 

y 
h 
i 

dt 3! dt? 5! d# 
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that Vo(t) is analytic in the neighborhood of #0. Let us first consider 


co 
I(z,t) as a function of x. If we replace Vo(s) by a series > ad» s” and make 


the transformation y = s/(t—s), we have 


T(n+1) Vr 


where ¢(¢) is an analytic function. 
Hence for values of 540, and for n sufficiently large, we have the 


inequality 


‘ 
| | $(t)} | <k nl, (29) 


where & is a constant independent of n. From this it follows, making use of 


Stirling’s approximation for n!, that 


lim 


qd” 1/n 


Similarly the nth root of the general term of the second series leads to 
the same limit and the convergence of I(z,¢) regarded as a function of z is 
established. Regarding it as a function of t, we may use 


| Az | 


as a dominating series for the first series and 


n! kn 


1 
n=1 (2n—1)! 


as a dominating series for the second. 
The case of a cable with distributed resistance R and capacity C with 
distributed leakage G per unit length leads similarly to the differential equa- 


tions: 


(Cz + G) V =— I. 


=0 
| 
— 
2n-1 
| 
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‘: If the cable is infinitely long we shall then have 
ce 
V = Vo(t) 
*(RB °(RB 
(RB)*/? (RB)*/? 
where we have used the abbreviation B= Cz + G. 
Making use of the expansion of the operator (d/dx +k)” given in the 
second section and writing A = G/C, this result may be put into the form, 
(RC)? 
30 
3 3/2 
Vr Vi 
[ The solution for the case where unit E. M. F. is impressed at 0 can 


be obtained with considerable ease by this formula. We first recall that: 


t ers 0 ers 
dg == f ds — f ds, 
0 Vt —S§ fee) Vit Vt —S§ 


0 ers 
Substituting this value for the integral in (30) and taking the indicated 
derivatives we get 


V = 4 E _a(RC)4 


+ 2-2-2 ] Vr 

Making the transformation t —s = é, (31) reduces to 

= VRC 
V = e-2#VRG 
4 36 
(RC)*c* (RC) Jer dé 


1(48)° 
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That this, together with the value of J(z,¢) derived from it, forms a 
solution of the original differential equations can be verified by direct sub- 
stitution. The boundary condition is obviously satisfied. 

The general case of the inductive cable can be treated similarly by opera- 
tional methods, but since it does not involve the use of fractional operators 
the discussion will be omitted. 

A third problem employing fractional derivatives is that of determining 
the voltage at a cable terminal when an E. M. F., f(¢),.is impressed on a long 
cable of distributed resistance R and capacity C per unit length through a 
condenser of capacity C,. We have just seen that the current entering a cable 
whose terminal voltage is V is VC/R 2 V and the current flowing into the 
condenser will be Coz {f(¢) —V} since f(t) —V is the voltage across the 
condenser. 

Equating these two values and dividing by z we obtain the following 
equation for the determination of V: 


(C/R C.) V=C,f(t). 


The solution is easily found since this is a fractional equation of order 
—4. From equation (18) we obtain 


0 RC? 


If f(x) —1, this reduces to 


R Co? Tv 0 V s 
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Classification of Quadrics in Hyperbolic Space. 


By JAMES PIERPONT. 


§1. Introduction. 


In this paper,* I wish to present an extremely elementary, though some- 
what lengthy method, of classifying quadric surfaces in H;. Instead of 
using the usual projective methods of Cayley and Klein, I make use of the 
more direct methods of Riemann as I have developed in three earlier papers.t 
A brief resumé of this method now follows. 

Let & 7, ¢ be three real numbers, we say they define a point, whose 
codrdinates are these numbers. How the codrdinates are chosen is not speci- 


fied. The metric is defined by 


Boe 


where +e, 
p= +1, do= + Vde + dy? + dt. 


ds 


H-straights ¢ are defined by 
5 fds = 0. 


To perform the variation it is convenient to introduce a new set of codrdinates 


= 


4 4 42 


Let us also set : 
[ab } = ab, + {ab} = ab, asbs — 


We find the codrdinates 2, - - 2, satisfy the relation — — R’. 


*Cf. P. Barbarin, “Etudes de Géométrie Analytique Non-Euclidienne,” Mém. 
Couronnés, Acad. Roy. Bruwelles, Vol. 60 (1900-1) ; J. L. Coolidge, Trans. Amer. Math. 
Soc., Vol. 4 (1903), p. 161; J. T. Bromwich, Jbid., Vol. 6 (1905), p. 275. 

t The Amer. Math. Monthly, Vol. 30, p. 425; Vol. 31, p. 26 (1923-4) ; Proceedings 
of the Mathematical Congress, Toronto (1924); Bull. Amer. Math. Soc., 1926. The last 
two have not yet appeared. 

+ For H-read hyperbolic, for e- read euclidean. 
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A relation [ar] =0 defines a H-plane; H-straights are the intersection 
of H-planes and conversely. 

In order to have a model of our H-space we will regard &, », ¢ as rectan- 
gular coérdinates of a point in e-geometry as in ordinary analytic geometry, 
Our H-space in this model is represented by the interior points of the e-sphere 

= (0. H-planes in the model are e-spheres cutting A 0 orthogonally, and 
H-straights are e-circles also cutting 4 = 0 orthogonally. 

Besides the points so far considered it is necessary to consider the points 
(é, », £) for which A=0. We call these ideal and define their 2 codrdi- 
nates by 


hence for ideal points {2} 0. We must also introduce certain points called 
imaginary ; their codrdinates have the form 


where ¢, ¢, are real numbers such that {¢?} — Hence imaaginary 
points also satisfy the relation {z?} —— R?. It is understood that when we 
speak of points we mean in general non-ideal points. 

If a is a point real or imaginary, the locus of the points z such that 
{ax} = 0, is called the absolute polar of a. We say x is conjugate toa. Ifa 
is imaginary its polar is real; if a is real its polar is imaginary. All H- 
straights perpendicular to a real plane pass through its pole. 

The four planes 7; x,—0O form a tetrahedron, the tetra- 
hedron. Each face 2,==0 is the polar of the opposite vertex which we 
call Ax. 

Let a:b denote the H-straight joining the points a, b. If x is a point 
on this straight we denote the length of the segment a, x in the sense a to # 
by (a, x) ete. Then 

sinh (2,b)/R sinh (a,x) /R 


= may + nby om] --- 4, 


The cross ratio of four points (a, b, 2, 2’) is n/m: n’/m’; it is harmonic 
when ——1. It is important to observe that a bisects the segment z, 2’, 
when (a, b, x, x’) is harmonic and 6 is conjugate to a. 


ion 
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The linear transformation 


== t,jm1---4, ai; real 
j 


defines an H-displacement or a motion if ds and = {a}. In 
particular any e-rotation about A, in our model is an H-displacement. The 


transformation 


= 2’, 2,’ cosh 6—z’, sinh 
=— x,’ sinh 6 + 2,’ cosh 0 


defines a rotation about the edge A, - Az of the r tetrahedron. 
Finally the angle @ between two curves meeting at a point (£7) is 
given by 


_[# #] 


16R* dé d& dn dy , dt dk 


Hence this angle in H-measure is the same as the angle in e-measure between 


the corresponding curves in the model. 


§ 2. Classification. 


Any quadratic relation between 2,‘ - - 2, defines a quadric surface. 
Here we consider only surfaces whose equations are homogeneous, Viz. : 


The a’s are real numbers and the determinant | ai; | 40. We call the quadric 
+ x2? + 2,42 =0 the absolute, it remains invariant under a dis- 


placement. 
If a= (am a) is a point we call 


the polar of a relative to F. Thus the plane {ax} —0 which we called the 
absolute polar of a, is the polar of a relative to the absolute. 

Suppose now the point g is such that its polar relative to F coincides 
with its absolute polar. If g is real, every chord of F through g is bisected 
by g; that is g is a center of F. If g is an imaginary point, its polar is a 
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plane of symmetry, i. e. all chords perpendicular to this plane are bisected 
by it. Such g points must satisfy the equations: 


(d11 91 + + Gisgs + drags = 0 
+ (@22 go + + = 0 

+ (dss —k) + dssgs = 0 (1) 
+ + + (as + gs = 0. 


We call these the g-equations. For these to hold it is necessary that & satisfies 


the equation 


which we call the H-equation; its roots we denote by a, b, ¢, d. 
If we subject the surface F to a displacement the new surface will have 
an equation 
F = a’ = 0 


whose H-equation will have the same roots as (2). By this displacement the 
surface F’, the point g and its polar will behave as a rigid body. If g is real, 
imaginary or ideal, it remains such after displacement. This- being so, we 
have only to consider the various possibilities that may arise and in each such 
case move the surface relative to its center, plane or axis of symmetry in order 
that its equation takes on a simple form.* The following theorems readily 
proved except the last, simplify the work: 


THEOREM 1. To two different roots of the H-equation do not correspond 
the same g. 


THEOREM 2. If ga, go correspond to two different roots a, b of the H- 
equation then {ga, go} =0. 


THEOREM 3. Moreover both ga, gv cannot be ideal. 


THEOREM 4. If gq is real, no gy can be ideal. 


* Dr. P. D. Schwartz in a paper given at the May meeting of the American Mathe- 
matical Society has shown how conics in H, may be classified by this method. Cf. 
Bull. Amer. Math. Soc., Vol. 32 (1926), p. 315. 


TurorEM 5. The H-equation has at least two real roots.* 
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I find now the following types of quadric surfaces, the numbers on the 


right are the types given by Bromwich. 


No g necessarily ideal. Roots of H all real. 


», arbitrary. 


Roots unequal. (k —a)(k—b)(k—c)(k—d) =0 
(1) ga real. ba,” + cx? + — = 0 
One pair of equal roots (k — a) (k—b)(k—c)* =0 
(2) ge real. ax,? + bx,” + = 0 
(3) gc imag., ga real. + 23”) + — az,’ = 0 
Two pairs of equal roots (k —a)*(k—b)?=0 
(4) ga real. b(x,? + 2,7) + — 2,7) =0 

or b(2,?-+ 2,”) + a(x3?— 2,7) =0 
(5) ga, go imag. + + 2737) = 0 
Triple Root (k—a)*(k—b) =0. 
(6) gp real. a(x,* + 22” + x7) — bry? =0 
ga real. ba,” + + — = 0 
(8) 9a, go imag. + bx.? + + — 0 
Quadruple Root (k—a)*=0 
(9) ga real + + — x,?) = 0 


(10) ga necessarily imag. + x2? + + = 0 


(11) ga imag. or ideal. 


Id 


Ic2 
IIIb 


Te 
IIIc 


IVb 


*T. J. Bromwich, Quadratic Forms. Nr. 3 of the Cambridge Tracts in Mathe- 
matics, Cambridge, 1906, p. 69. The proof as given by Bromwich rests on the theory 
of elementary divisors. A very simple proof of a generalization of this theorem which 


I have found will be given in another paper. 


I 
Tb2 
Ib1 
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A pair of conjugate imaginary roots. 
Two real unequal roots (k —a) (k —b)(k —a—iB) (k —a + i8)=0, B0 


(12) ga, g» imag. + + + Va? + 0 
(18) a=b. ga imag. + 227) + + + IIb 


One g necessarily ideal. Roots all real. 


One pair of equal roots (k —a)(k—b)(k—c)? =0 

(14) ge ideal, ga, go imag. + + + — 0 TITa 
Two pairs of equal roots (k —a)?(k—b)? = 

(15) ga imag., go ideal. a(x,? + + + = 0 
Triple root. (k—a)*(k—b) =0 


(16) ga ideal, gy imag. 


bay? + axe? + + — — IVa 


We note that Bromwich has no case corresponding to (5). 
We now illustrate on a few cases how the reduction is effected. 


Type 6. We move F so that ga coincides with A,=— (0, 0, 0, R). The 
g-equations give 
= Ago = O43 = 0,7 — (3) 


We rotate F about A, till g, coincides with A, = (iF, 0,0,0). Then the new 
coefficients satisfy (3) and in addition the g-equations give 


0, = Ag, = 0. F now has the form 


bay? + — + = 0. 
It will not confuse the reader if we omit from the ai; here and in the following 
the accents. We now rotate F about the edge A,- A, of the + tetrahedron 
through the angle 6 determined by 


tan 6 = 
— Are 


when @22 > d33, in the latter case thru 645°. Then F has the form 


+ + — =0, 


= 
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whose H-equation is (k —a)(k—b) (k 


(k = 0. By hypothesis 
H=(k—a)*(k—b)?; hence az. = a, a33 = b or doz = 0, =a. We thus 
get type 6. 


Type We move F' so that ga becomes A; =(iR,0,0,0). The g-equa- 
tions give d21 = =0, We rotate F about A, till coin- 
cides with A, = (0, 0, 0). Then ay. = = = 0. It is 
easy to show that | | < | @s3 + |. We can therefore rotate F about the 
edge A, Az so that as, 0 if it is not already 0. F now has the form 


whose H-equation is (k — a) (k —b) (k — (k + = 0. Hence either 


The g-equations in case 1° show that g, may be taken real, which is contrary 
to hypothesis. Hence only 2° exists and this is type 7. 


Type 10. We move F so that ga coincides with A, = (1R, 0, 0,0). The 
g-equations give 4, = 4, = 0. The plane z,—0 cuts F in 
the conic ® whose g-equations are obtained by deleting the first row and col- 
umn of (1). The H-equation of this system is Hy = (k—a)*=0. 


Suppose now that the g-equations of ® admitted a solution for k =a, 
to which corresponds a real non-ideal point y= (ge, gs, gs), then g’= 
(0, 92, Js, Js) 18 a solution of the g-equations of F for k =a, and F belongs 
to type 9. We will suppose therefore that the g-equations of ® admit only an 
imaginary y, then g’ is imaginary. We rotate F about A, till g’ coincides 
with A,= (0, 0, 0), then the g-equations give = a4. = 0, 
(go = 0 and F has the form 


az,? Aa,” + Ba? = 0 
whose H-equation is 
(k —a)*{k® — (a— + a} — (k —a)*=0. 


Comparing coefficients gives a— B= 2a, y—aB—a?; hence a=a-+ey, 
8 = ey —a, & =1, and F has the form 


a(a,? + + (a + ey) (a— ey) = 
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This surface cuts the edge 7, = 0, 7, = 0 in two points given by 
(a + ey) 23” (a ey) + = 0, 
one of which is (0, 0, 2ef, 2h). If the other is P, we rotate F about A,: A, 


so that P coincides with A,=— (0, 0, 0, R) when a=ey. F now has the 


type (10) 
As further illustration of our method let us investigate the cases which 


arise when 
H = (k—a)*=0 and g ts necessarily ideal. (C) 


We begin by rotating F about A, so that ga coincides with L = (0, 0, 2R, 2R), 
The g-equations give noW — 413, —=—A23, O31 = A—A33, = 
as3— 2a. The corresponding H-equation is 


(k— a)? Qi2 = (k—a)*=0. 
G12 22 — k 
Comparing coefficients gives a1: 22 = @, = 9, Gis = B, Gos = y, hence F 
has the form 
a(x," + 22") + + (y — 2a) + — 
+ 2y — + 2(a— = 0 (4) 
This cuts the edge x, =x, 0 in ZL and in another point determined by 
nts” + (yn — + 2(a— = 0. 
Suppose (C,1) that PAM = (0, 0, —2R, 2R). We rotate F about A,- A, 
so that P coincides with A, = (0, 0,0, R). These codrdinates set in (4) give 
n = 2a, and F reduces to the form 
A(x? + 227) + + 2B — + Py — — = 0 
whose first two g-equations for k =a give 
B(9s—9s) =9, —gs) = 0. (5) 


Suppose (C,1, a) that B or y is #0, then g;=g,. The last two g-equations 
give 
BG: + 9. (6) 


Thus the g-equations for k —a admit the solution 


R 


Ja = » — 0, 0) 


i 
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This is contrary to our hypothesis (C’) that ga is necessarily ideal; hence case 
(C, 1, a) does not exist. 


Suppose (C, 1, b) that B=y=0. The g-equations for k—a admit 
the solution ga = (1R, 0, 0, 0) = A; which is also contrary to (C). 


Suppose (C, 2) that PM. In this case a rotation about A,-° A, does 
not move M and the foregoing does not apply. We may proceed as follows. 
Setting the codrdinates of M in (6) gives y =a and F has the form 


+ 22") + ary? — ar,” + 28( 2173 — 4124) + — 0 


whose g-equations for k =a give again (5), (6). Hence if B or y is ~ 0 the 
g-equations are satisfied by ga = (py, — 8, 0, 0), p as above, i. e. by a non- 
ideal point contrary to (C). Finally if 8 —y —0, F reduces to type 9 whose 
g-equations are satisfied by any g, contrary to (C’). Thus case (C) does not 


exist. 
In a similar manner all possible cases may be treated. 
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